Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,851 A→T A182A (GCA→GCT mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,8510AT92.5% 106.7 / ‑0.8 40A182A (GCA→GCTmogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (16/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAAC  >  NC_000913/9607‑10093
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
accAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAACCGGCACCAGAAGTGGTTGCTGCATTc                                                                                                                                                                                                                                              >  2:163635/3‑251 (MQ=255)
   cgttgcgccatatacgCCTGCATCTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGTGTGATTCGCAAACAGGCGCTGAGCTTTAACTTACCCGGTCAGCCGAAGTCTTTTAAAGAGACGCTGGAAGGTGTGAAGGACGCGTCGGGTAACGGTGTTGTACACGTTATTTTTGCCACCGTACCGTACTTTATTCAGTTGCCGTAAGGGCCATGCGTAGAACCGGCACCAGAAGTGGTTGCTGCATTcc                                                                                                                                                                                                                                             <  2:69881/234‑2 (MQ=255)
    gATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGAc                                                                                                                                                                                                                                          <  2:395726/251‑1 (MQ=255)
     aTGCGCCAGATAAGCCTTCATGTTGTATCAACGGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGGTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGAcc                                                                                                                                                                                                                                         <  2:342278/251‑1 (MQ=255)
       gcgcCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTAATGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGGGGTGGCAGCATTCAGCCCg                                                                                                                                                                                                                                        >  2:186499/1‑250 (MQ=255)
                   ccTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAAAAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTACAGAGACGCTGGACGGGGTGACGGACGCTGCGGGTAACGTTGTGCTACACGGTATTTTTGCCAGCGTACCGTCCTGCATTCAGTTGCTGTAAGGGCCATACGTTGAACCGGCACAAGAGGTGGTTGCTGCATTCAGACCACAGAGTGCAAg                                                                                                                                                                                                                             >  1:451568/1‑249 (MQ=255)
                         ttttGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACg                                                                                                                                                                                                                     >  1:456589/1‑251 (MQ=255)
                            tGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGAAAGACGCGACGtt                                                                                                                                                                                                                   >  2:191727/1‑250 (MQ=255)
                                    cTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATa                                                                                                                                                                                                           >  2:246031/1‑250 (MQ=255)
                                                 gcgTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAGCTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGATGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCg                                                                                                                                                                                             <  1:387530/251‑1 (MQ=255)
                                                 gcgTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCg                                                                                                                                                                                             <  1:401209/251‑1 (MQ=255)
                                                     cAGGTTGGTGTGATTCGCAAACAGGCGCTGAGCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGGGACGTTAGCGAATAAAAAAATCCCCCCGAGCg                                                                                                                                                                                         <  1:437963/251‑1 (MQ=255)
                                                       ggttggtgTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGAGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAAGCGACGTTAGCGAATAAAAACATCCCCCCGAGCggg                                                                                                                                                                                       >  2:150834/8‑251 (MQ=255)
                                                         gtggtgTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggggg                                                                                                                                                                                     <  2:332432/246‑1 (MQ=255)
                                                           ggtgTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGGAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgggggg                                                                                                                                                                                    <  2:81152/247‑1 (MQ=255)
                                                             tgTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCTTACTGCATTCAGTTGCTGGAAGGGCCATAAGTTGAAAAGGCCCCAGAGGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGATAGCGGATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                                                                                                                 >  2:406384/2‑251 (MQ=255)
                                                                                  gATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGAt                                                                                                                                                              <  1:206645/249‑1 (MQ=255)
                                                                                    tCCTTAACTTACCCGGACAGCCGAAGTCTATTAAAGAGACGCTGGCTGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGGACTGCATTCAGTTGCTGGCAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGTAGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGCGGTATCAAAACGATTAGTGGGATTCAc                                                                                                                                                          >  2:350554/1‑251 (MQ=255)
                                                                                        tAACTTACCCGTTCTTCCGAAGTCTATTAAAGAGACGCTTGAATGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTTTTTTTTCCAGCGTACCGTACTGCATTTATTTGCTGGAATGGCCATACGTTGAAACGGCACCAGATGTGGTTGCTGCATTCAGTCCTAAGAGTGCAAGTCGCGACGTTAGCGAATAAAAAATTCCCCCCGAGCGGGGTGATCTCAAAACAATTAGTGGGATTCACCa                                                                                                                                                        <  1:459686/249‑1 (MQ=255)
                                                                                            ttACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGc                                                                                                                                                  >  1:477654/1‑251 (MQ=255)
                                                                                                  ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGAc                                                                                                                                                                                                                                          >  2:244565/1‑157 (MQ=255)
                                                                                                  ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGAc                                                                                                                                                                                                                                          <  1:244565/157‑1 (MQ=255)
                                                                                                   gTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                                                                                                                           <  1:355330/251‑1 (MQ=255)
                                                                                                                 aTTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACt                                                                                                                               <  1:163635/249‑1 (MQ=255)
                                                                                                                           cGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                                                                    >  1:79416/1‑250 (MQ=255)
                                                                                                                                  aGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                                                                                                            >  1:391422/1‑251 (MQ=255)
                                                                                                                                         aaGGACGCTGCGTGTACCGTTGTGGTTCTCGGTATTTTTTCCAGCGTACCGTTCTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCATAAGTGGTAGCTGCCTTCAGACCGAAGAGTGCAAGACTCGACGTTAGCTAATAAAAAAATCCCCCCGAGCGGGGTGATCTCAAAACAATTAGTGGGATTCACCAATCTGCAGAACGGTGCGACCAAACTGCTCGTTCATTACTTCACCCATCgtc                                                                                                     <  2:469835/251‑3 (MQ=255)
                                                                                                                                             ccGCTGCGGGTTACGTTGTGGTACACGGTATTTTTGCCCGCGTACCGTACTGTATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATGCAGCCCTAAGAGTGCAAGACTCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGTCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCGCCagat                                                                                                 <  2:451568/250‑1 (MQ=255)
                                                                                                                                               gCTGCGGGTAACGGTTTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAACGGCAGAACGGTGCGACCAAACTGCCCGTTCAGTACTTCACCCATCGCCagatag                                                                                               <  2:153427/251‑1 (MQ=255)
                                                                                                                                                         aCGTTGTGGTACACGGTATGTTTGCCAGCGTACCGTACTGCATTCAGTTGGTGGAAGGGCCATCCGGTGAAACGGCACCAGAAGTGGTTGCTTCATTCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAACCAATTAGTGGGATTCACCAATCGGCCGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACTCATCGCCTGATAGATTGCACTgg                                                                                     <  1:470116/251‑1 (MQ=255)
                                                                                                                                                                   ccactgtagttttGCCAGCGTACCGTACTGCATTCAGTTGCTGAAATGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCTAGCGGGGGTATCTCAAAACAATTAGTGGGATTCACCAATCGTCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGa                                                                            <  1:38289/241‑1 (MQ=255)
                                                                                                                                                                                         ccgtgctgcagtccattgattgAAGGGCCATACGTTGAAACTGCACCATAAGTGGTTTCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAATCAATTAGTGGGATTCACCACTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCCCCCCAGATCAGCCCAATCCAGCCGGCAAAg                                                     <  2:156575/230‑1 (MQ=255)
                                                                                                                                                                                               tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGCTCAGCCCAATCCATCCGGCAACGTGgatg                                               >  1:38041/1‑251 (MQ=255)
                                                                                                                                                                                                          gCTGGAAGGGCCATACGTTGAAACTGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAGCGCCAGATAGATTGCACTGGCACCACAGAGCAGACCAATCCAGCCGGCAAAGTGTATGATTGCGGCGt                                     <  1:422348/250‑1 (MQ=255)
                                                                                                                                                                                                                     cATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGt                         >  1:313888/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGtt                        <  1:191727/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGCACCAGAAGTGGTTGCTGCATGCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACATATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCca             <  2:291351/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCACCGGAAGTGGTTGCTGCATTCAGACCGCAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACAACTCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGcc              >  2:214286/1‑249 (MQ=255)
                                                                                                                                                                                                                                     accagAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGTATAAAACAATCCCCCCGAGCGGGGGGATCTAAAAACAATTAGTGGGATTCACAAATCGCCAGGACGTTGCGACCAAACTGCTCGTTCAGGACTTCACCCATCTCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGCCAAAGTGGATGATTGCGGCGTACCCGGCACTGTTACCGATCGCcagc           >  1:38951/5‑249 (MQ=255)
                                                                                                                                                                                                                                              ggTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAAc  >  2:423885/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
AACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAAC  >  NC_000913/9607‑10093

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: