Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,812 C→A 8.4% H275N (CAT→AAT)  lgoD → L‑galactonate oxidoreductase; L‑gulonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,596,8120CA 8.4% 2040.3 / 34.5 649H275N (CAT→AAT) lgoDL‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (375/198);  new base A (2/74);  total (377/272)
Fisher's exact test for biased strand distribution p-value = 1.83e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 659 total aligned reads displayed.

AGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAG  >  NC_000913/4596663‑4596961
                                                                                                                                                     |                                                                                                                                                     
 gTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                                                      >  1:7743278/1‑150 (MQ=255)
  tttGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCAt                                                                                                                                                     >  1:8468780/1‑150 (MQ=255)
  tttGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCAt                                                                                                                                                     <  1:7344090/150‑1 (MQ=255)
   ttGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATa                                                                                                                                                    <  2:7154051/150‑1 (MQ=255)
      gtggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaa                                                                                                                                                 >  2:5666059/1‑150 (MQ=255)
        ggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaaa                                                                                                                                               <  2:2008493/150‑1 (MQ=255)
          ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                             <  1:7567195/150‑1 (MQ=255)
          ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                             <  1:8754428/150‑1 (MQ=255)
          ttCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGa                                                                                                                                             >  1:9791238/1‑150 (MQ=255)
            cGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGaaa                                                                                                                                           <  2:5951450/150‑1 (MQ=255)
                      gattgTGATCGACGCTACTGGTAATAAACATGCGATTATTAACACCGTGATTTTGTTTCGTCACGGCGGCACGGTGTTTTTTGTCGGCCTGTTTAAAGGTGTGTTGCTTTTCTCCGTTCCGTATTTCCATAAAAAAGAAACGACgatgat                                                                                                                                 <  2:8007321/146‑1 (MQ=255)
                       aaaGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGatgatg                                                                                                                                >  1:7978680/1‑150 (MQ=255)
                       aaaGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGatgatg                                                                                                                                >  1:1759950/1‑150 (MQ=255)
                       aaaGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGatgatg                                                                                                                                >  2:3476850/1‑150 (MQ=255)
                         aGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGAGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATggg                                                                                                                              >  1:5029373/1‑150 (MQ=255)
                           tGATCGACGCGACTGGTAATCAACATGCGATTAATAACACGGTTATTTTGATTCGTCACGGCGGCAGGGTGGTATTTGTGGGCGTTTTTAATGTTGAGTTGCAGTTCTCCGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCa                                                                                                                            <  2:5074425/150‑1 (MQ=255)
                              tCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGcc                                                                                                                         >  1:3102702/1‑150 (MQ=255)
                                gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                                                       >  1:10890253/1‑150 (MQ=255)
                                  cgcgACAGGTAATCAACTTGCGATGAATAACACGGTTATTTTTTTTGGTCACGGCGGCAGGGTGTTTTTTGTCGGCCTTTTTAATGGTGAGTTGCTGTTCTCGGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCaa                                                                                                                     <  2:499267/150‑1 (MQ=255)
                                   gcgACAGTTAATCAACATGCGTTGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTTTTTGTCGGCCTGTTTAATGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                                                    <  2:7228789/150‑1 (MQ=255)
                                        aGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccac                                                                                                               >  1:3580930/1‑150 (MQ=255)
                                         ggTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacg                                                                                                              >  2:7412051/1‑150 (MQ=255)
                                          gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                             >  1:1501959/1‑150 (MQ=255)
                                          gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGAGCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                                                             >  2:9439428/1‑150 (MQ=255)
                                           tAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgcc                                                                                                            >  2:6945163/1‑150 (MQ=255)
                                            aaacAACTTGCGTTGAATACCACCGTGAATTTGTTTCGTCACGGCGGCAGGGTGTTTTTTGTCGGCCTTTTTAAAGGTGTTTTGCTGTTCTCCGATCCGTAATTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCg                                                                                                           <  2:2565275/147‑1 (MQ=255)
                                            aaTCAACATGCGATTAATAACACCGTTAATTTTATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTTTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAAGTCCATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCg                                                                                                           <  2:8948935/150‑1 (MQ=255)
                                            aaTCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCg                                                                                                           >  2:2168298/1‑150 (MQ=255)
                                              tCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGa                                                                                                         >  2:8572697/1‑150 (MQ=255)
                                                   aTGCGATTAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGGGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGAtt                                                                                                    <  2:9807202/150‑1 (MQ=255)
                                                    tGCTATTAATAACACCGTTAATTTGATTCGTCACGGCGGCAGGGTGTTTTTTGTCGGCCTTTTTAACGTTGTGTTGCTGTTCTCCGATCCGTAATTCAATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGAttt                                                                                                   <  2:2470116/150‑1 (MQ=255)
                                                      ctaTGAATAACACCGTTAATTTTATTCGTCAGGGCGGCAGGGTGGTATTTGTCGGCCTTTTTAAAGTTGAGTTGCAGTTCTCCTATCCGTATTTCCAAAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                                 <  2:6735937/148‑1 (MQ=255)
                                                      cGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                                                 >  1:5253250/1‑150 (MQ=255)
                                                       gATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGc                                                                                                >  1:11807473/1‑150 (MQ=255)
                                                        aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                               >  2:954121/1‑150 (MQ=255)
                                                        aTGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCt                                                                                               >  1:4439848/1‑150 (MQ=255)
                                                           aaTAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaaa                                                                                            >  1:2442097/1‑150 (MQ=255)
                                                            aTAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAg                                                                                           >  1:4600759/1‑150 (MQ=255)
                                                             tAACCCCGTTATTTTGATTCGCCCGGGCGGCCGGTTTTTTTTTGTGGGCCTTTTTAATGTTGTGTTGCTTTTCTCCTAACCTTAATCCAATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgt                                                                                          <  1:8162977/150‑1 (MQ=255)
                                                               acacCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcg                                                                                        >  2:1477876/1‑150 (MQ=255)
                                                                 acCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAgtcggt                                                                                      >  2:6831698/1‑150 (MQ=255)
                                                                    gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                                                   >  2:4331096/1‑150 (MQ=255)
                                                                     tGATTTTGATTCGTCACGGCGGCACGGTGTTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGAc                                                                                  <  2:12141026/150‑1 (MQ=255)
                                                                         tttGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGAt                                                                              >  2:11877400/1‑150 (MQ=255)
                                                                            gATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGc                                                                           >  1:4199079/1‑150 (MQ=255)
                                                                             aTTCTTCACGGCGGAACTGTGTTATTTTTCGCCATTTTTAAAGGTTATTTCCAGTTCACCCAACCGTAATACAATAAAAAAAAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCg                                                                          <  2:4380792/150‑1 (MQ=255)
                                                                              ttCGGCCCGGCGGCACGGTGTTATTTGTCGGCCTGTTTAACGGTGAGTTGCAGTTCTCCGAGCCGTAATTCAATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCgg                                                                         <  2:7777566/150‑1 (MQ=255)
                                                                                cGTCCCGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                                                       <  2:10214252/150‑1 (MQ=255)
                                                                                cGTCAGGGCGGCACGGTGGTATTTGTCGGCCTTTTTAATGTTTTGTTGCTGTTCTCCGATCCGTATTTCCTTAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                                                       <  2:3992202/150‑1 (MQ=255)
                                                                                cGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGTGTTGCAGTTCTCCGATCCGGAATCCAATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                                                       <  2:8153297/150‑1 (MQ=255)
                                                                                     cggcggCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGaaa                                                                  >  2:3677433/1‑150 (MQ=255)
                                                                                       ggggccgggtgttATTTGTCGGCCTTTTTATCGGTGTGTTCCAGTTCCCCTACCCGGAATTCCATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGaaaaa                                                                <  1:1362269/143‑1 (MQ=255)
                                                                                              ggtggtTTTTGTGGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGACCCGCAATCCCAAAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgc                                                         <  2:11633982/150‑1 (MQ=255)
                                                                                                 tgtttttGTCGGCATTTTAAATGGTGTGTGGCAGTTCCCGGATCCGTAATTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                      <  2:1571096/146‑1 (MQ=255)
                                                                                                 ggtATTTGTCGGCCTGTTTAACGGTGAGTGGCTGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGa                                                      <  2:827303/150‑1 (MQ=255)
                                                                                                  gtATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGAc                                                     >  1:12211167/1‑150 (MQ=255)
                                                                                                   tATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACa                                                    >  1:6877806/1‑150 (MQ=255)
                                                                                                    ttttGTCGGCATTTTTAACGTTTATTTGAAGTTCTCCTATCGGGAATTCCATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                                   <  1:2765279/149‑1 (MQ=255)
                                                                                                    aTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                                   >  2:7612782/1‑150 (MQ=255)
                                                                                                    aTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                                   >  2:5238608/1‑150 (MQ=255)
                                                                                                    aTTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACat                                                   >  1:6105095/1‑150 (MQ=255)
                                                                                                     tttGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatg                                                  >  2:11400560/1‑150 (MQ=255)
                                                                                                      ttGTCGGCCTTTTTAAAGGTGAGTTGCAGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                                 <  2:1845525/150‑1 (MQ=255)
                                                                                                      gtgtgGGCCTTTTTAATGGTGCGTGGCTGTGCTCCGATCCGGAATTCCACAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                                 <  1:10543295/149‑1 (MQ=255)
                                                                                                      gtGTCGGCCTGTTTACCGGTGCGTTGCAGTTCTCCGAGCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatga                                                 <  2:1712354/149‑1 (MQ=255)
                                                                                                        gTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatgatg                                               >  2:10735214/1‑150 (MQ=255)
                                                                                                           ggCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTa                                            >  2:4738536/1‑150 (MQ=255)
                                                                                                            gcatttttAATGGTGTGTTGCAGTTCTCCGATCCGTAATGCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTaa                                           <  1:372562/145‑1 (MQ=255)
                                                                                                             ccTGTTTAACGTTTCTTTGCCGTCCCCCGAGCCGTAATCCCATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc                                          <  2:5448798/150‑1 (MQ=255)
                                                                                                                 tttAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCAt                                      >  1:7089493/1‑150 (MQ=255)
                                                                                                                  ttAACGGTGTGTTGCATTTCTCCGATCCGGAATTCCAAAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc                                     <  2:4523126/150‑1 (MQ=255)
                                                                                                                  ttAACGGTGAGTTGCTGTTCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc                                     <  2:3823537/150‑1 (MQ=255)
                                                                                                                  ttAAAGGTGAGTTGCATTTCTCCGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc                                     <  2:6051219/150‑1 (MQ=255)
                                                                                                                   tAATGGTGAGTTGCAGTTCCCCTATCCGGAATTCAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg                                    <  1:2563865/150‑1 (MQ=255)
                                                                                                                   tAACGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg                                    <  2:9194969/150‑1 (MQ=255)
                                                                                                                            gTTGCTGTTCTCCGATCGGGAATTCCATAAAAAAGACACGACGATGACGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGtt                           <  1:5813163/150‑1 (MQ=255)
                                                                                                                             ttGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                          >  1:11194810/1‑150 (MQ=255)
                                                                                                                             ttGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTc                          >  2:11155478/1‑150 (MQ=255)
                                                                                                                              tGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAAAGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCg                         >  2:10632906/1‑150 (MQ=255)
                                                                                                                                cAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGcc                       <  2:10995646/150‑1 (MQ=255)
                                                                                                                                 aGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCa                      <  1:949422/150‑1 (MQ=255)
                                                                                                                                   ttCTCCGATCCGTAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACg                    <  1:12036959/150‑1 (MQ=255)
                                                                                                                                    tctcCGATCCGGAATTCCAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGc                   <  1:6348073/150‑1 (MQ=255)
                                                                                                                                    tctcCGAGCCGCAATACAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGc                   <  1:5513144/150‑1 (MQ=255)
                                                                                                                                     ctcCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTCTCCGTTCGCCACGCt                  >  1:3988562/1‑150 (MQ=255)
                                                                                                                                       cctaTCCGTAATTCCATAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTgg                <  1:9230030/147‑1 (MQ=255)
                                                                                                                                        cGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc               >  2:4814100/1‑150 (MQ=255)
                                                                                                                                        cGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc               >  1:1000523/1‑150 (MQ=255)
                                                                                                                                        cGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc               >  2:10010492/1‑150 (MQ=255)
                                                                                                                                        cGAACCGTAATTCAAAAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGc               <  2:4161004/150‑1 (MQ=255)
                                                                                                                                         gttCCGTAATTCCATACAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCa              <  1:10202539/148‑1 (MQ=255)
                                                                                                                                         gATCCGTATTTACATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCa              <  2:6364907/150‑1 (MQ=255)
                                                                                                                                         gATCCCGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCa              >  1:11393456/1‑150 (MQ=255)
                                                                                                                                          aTCCGGAATTCCATAAAAAAGAAACGACGATGATGGCAGCCGCAACGCCACGCCGGAAGTTTTTGCTAAAGTCGGTCGACTGATGGCGGTAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGa            >  2:3021096/1‑150 (MQ=255)
                                                                                                                                           tCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGa            >  1:8878739/1‑150 (MQ=255)
                                                                                                                                            ccGTAATTCAAAAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGaa           <  1:1889382/150‑1 (MQ=255)
                                                                                                                                                aaTTCCAAAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCt       <  1:4738536/150‑1 (MQ=255)
                                                                                                                                                  ttCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTAc     >  1:2021630/1‑150 (MQ=255)
                                                                                                                                                  ttCCATAAAAAAGAAACAACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTAc     >  1:6456507/1‑150 (MQ=255)
                                                                                                                                                   tCCATACAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACg    <  1:2329874/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAG  >  NC_000913/4596663‑4596961

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