Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,541,452 G→A 6.8% Q39* (CAG→TAG)  panE ← 2‑dehydropantoate 2‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,541,4520GA 6.8% 169.2 / 2.1 59Q39* (CAG→TAG) panE2‑dehydropantoate 2‑reductase
Reads supporting (aligned to +/- strand):  ref base G (23/32);  new base A (2/2);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCAATCGCCGGGGCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCCGTTATAATGAGCCCG  >  NC_002947/1541304‑1541597
                                                                                                                                                       |                                                                                                                                                 
gCAATCGCCGGGGCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGc                                                                                                                                                  >  1:1586664/1‑150 (MQ=255)
      gCCGGGGCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGt                                                                                                                                            >  2:149066/1‑150 (MQ=255)
       ccGGGGCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGtt                                                                                                                                           <  1:861727/150‑1 (MQ=255)
         gtggCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGc                                                                                                                                         <  1:328142/148‑1 (MQ=255)
         ggggCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGc                                                                                                                                         >  2:70578/1‑150 (MQ=255)
         ggggCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGc                                                                                                                                         <  2:1406398/150‑1 (MQ=255)
           ggCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGccc                                                                                                                                       <  1:70578/150‑1 (MQ=255)
                     tAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCCgg                                                                                                                             >  2:188041/1‑150 (MQ=255)
                     tAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAgg                                                                                                                             >  2:1453631/1‑150 (MQ=255)
                          cTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAggcaagg                                                                                                                        >  2:584094/1‑145 (MQ=255)
                          cTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAg                                                                                                                        >  1:260558/1‑150 (MQ=255)
                           ttGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTCAGACGTTTGCTCTTCCGATCTCCGCTACGC‑‑GTTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAgg                                                                                                                       <  1:1197811/150‑1 (MQ=255)
                           tgcaaGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAgg                                                                                                                       <  1:1227076/146‑1 (MQ=255)
                                ggCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGcac                                                                                                                  >  1:1823611/1‑150 (MQ=255)
                                   caccaGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGc                                                                                                               >  1:258231/1‑150 (MQ=255)
                                        gcagcCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGc                                                                                                          <  1:1228278/150‑1 (MQ=255)
                                        gcagcCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGc                                                                                                          <  1:567803/150‑1 (MQ=255)
                                            ccGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTgc                                                                                                      <  1:1898259/150‑1 (MQ=255)
                                             cGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTTGTCCTGCTCGACCAGGGATCGGTCAGCGGGCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGCTCAAGGGGACCGCCTTGCCTgcg                                                                                                     >  1:119626/1‑150 (MQ=255)
                                              ggTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTgcgc                                                                                                    <  1:669493/150‑1 (MQ=255)
                                               gTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGcgcg                                                                                                   >  2:753859/1‑150 (MQ=255)
                                                tGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGt                                                                                                  <  1:163387/150‑1 (MQ=255)
                                                     cGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTACAGCCGTTGCCCGTCACGCAGGATCAgg                                                                                                                       >  2:1352723/1‑124 (MQ=255)
                                                     cGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTACAGCCGTTGCCCGTCACGCAGGATCAgg                                                                                                                       <  1:1352723/124‑1 (MQ=255)
                                                             gcgcTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTc                                                                                                                                    <  1:982306/103‑1 (MQ=255)
                                                             gcgcTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTc                                                                                                                                    >  2:982306/1‑103 (MQ=255)
                                                               gcTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGc                                                                                                                                         >  2:1470820/1‑96 (MQ=255)
                                                               gcTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGc                                                                                                                                         <  1:1470820/96‑1 (MQ=255)
                                                               gcTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGcc                                                                                   >  1:950265/1‑150 (MQ=255)
                                                                cTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGCTTTGGCCCTGCACGACCAGCCTCAGGCCCCCGTCCTCCTGATAGTCCGGCAGCCGTTTCCCGTCAGGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCCGGccc                                                                                  >  2:709755/1‑150 (MQ=255)
                                                                    tgCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAg                                                                              >  1:1148545/1‑150 (MQ=255)
                                                                    tgCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAg                                                                              >  2:147821/1‑150 (MQ=255)
                                                                         tt‑‑‑tCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGc                                                                         <  2:1650873/150‑1 (MQ=255)
                                                                                 gCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAgg                                                                    <  1:1077993/150‑1 (MQ=255)
                                                                                 gCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAgg                                                                    <  2:1304393/150‑1 (MQ=255)
                                                                                    ggTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGACGTTGCCCGACACGCAGGATCAGGCGCACAGCCGTGCCTGCTCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTa                                                                 <  2:857970/150‑1 (MQ=255)
                                                                                      tATGGCGTAGCGCCGGGTTTGGTCCTGCTGGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCTTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCACGGCTAcc                                                               >  2:296538/1‑150 (MQ=255)
                                                                                       aTGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCg                                                              <  1:147821/150‑1 (MQ=255)
                                                                                             tAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCg                                                        >  2:82131/1‑150 (MQ=255)
                                                                                              ggcgcCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCga                                                       <  1:1932878/149‑1 (MQ=255)
                                                                                                cgcCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTGCCGGCGCCgaga                                                     <  2:143725/150‑1 (MQ=255)
                                                                                                  ccGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTCGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAat                                                   >  2:752840/1‑150 (MQ=255)
                                                                                                        ttGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCa                                             <  1:1559812/150‑1 (MQ=255)
                                                                                                        ttGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTCGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCa                                             <  1:752840/150‑1 (MQ=255)
                                                                                                            tCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGtg                                         <  2:258231/150‑1 (MQ=255)
                                                                                                                gCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAg                                                                     <  2:1102949/118‑1 (MQ=255)
                                                                                                                gCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAg                                                                     >  1:1102949/1‑118 (MQ=255)
                                                                                                                gCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTACAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGtgctgc                                     >  1:1824687/1‑150 (MQ=255)
                                                                                                                  tCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTc                                   <  2:1148545/150‑1 (MQ=255)
                                                                                                                         gggTCAGGCCACCGGCCTGCTGATAGGCCTACAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCg                            <  2:1824687/150‑1 (MQ=255)
                                                                                                                           gTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGtt                          <  2:1586664/150‑1 (MQ=255)
                                                                                                                              aGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGaa                                                    >  2:1618971/1‑121 (MQ=255)
                                                                                                                              aGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGaa                                                    <  1:1618971/121‑1 (MQ=255)
                                                                                                                                 ccACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCaa                    <  1:753859/150‑1 (MQ=255)
                                                                                                                                       gCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCCGtt              <  2:260558/150‑1 (MQ=255)
                                                                                                                                            ctgATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCCGTTATAAt         <  1:1453631/150‑1 (MQ=255)
                                                                                                                                              gATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGTCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCAGTTATAATGa       <  2:784314/150‑1 (MQ=255)
                                                                                                                                              gATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCCGTTATAATGa       <  1:885991/150‑1 (MQ=255)
                                                                                                                                                   gCCTGCAGCCGTTGCCCGTCACGCAGGGTCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGCCGCCGAGAATATGCCAGGTGCTGCCCATATGCGTTCCGGAACCGGTTATAATGAGCCCg  >  2:372333/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GCAATCGCCGGGGCAGCGTCATAGGCCTTGCAGGCCACCAGCAGCCGGTGGATCGGGCCTTGCGCTTGTGCCGTT‑‑‑TCGGCCGGTATGGCGTAGCGCCGGGTTTGGTCCTGCTCGACCAGGGTCAGGCCACCGGCCTGCTGATAGGCCTGCAGCCGTTGCCCGTCACGCAGGATCAGGCGCACAGCCTTGCCTGCGCGTGCCAGCCGGCAGGCCCACAGACTGCCAAGGCTACCGGCGCCGAGAATATGCCAGGTGCTGCTCATATGCGTTTCGCAACCCGTTATAATGAGCCCG  >  NC_002947/1541304‑1541597

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