Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,354 A→G 5.8% E272G (GAA→GGA)  waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,3540AG 5.8% 313.7 / 7.4 121E272G (GAA→GGA) waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base A (52/59);  new base G (4/3);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 122 total aligned reads displayed.

CGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATC  >  NC_000913/3809207‑3809500
                                                                                                                                                   |                                                                                                                                                  
cGCTGCGCCGCCAGTGGGCTCCACACAGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGaa                                                                                                                                                   <  2:530193/149‑1 (MQ=255)
 gCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGATTGTGGTGATCGCAGCACATCAGGCATTTTTACAGCAATTCCCGAATTAATTGCTCATCCTGGTACCCCGTCAGCAGGAAc                                                                                                                                                  >  2:1521904/1‑149 (MQ=255)
  cTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATACGGCAcc                                                                                                                                                 >  2:1456491/1‑148 (MQ=255)
       ccgccAGTGGGCACCACACCGCCCGGTATGGGTTGCCACCGGCACTCACGAAGGCGAAGCGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGTATTTATTGCTCCTCCTGGTACCCCGTCATCCGGACCGCTTcc                                                                                                                                            >  1:1052975/1‑149 (MQ=255)
        cgccAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTccc                                                                                                                                           >  1:545499/1‑149 (MQ=255)
           gatgtggCTCCACACCGCCCGGTATGGATTGCCACAAGCACTCACGAAGGCTAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACGCCGTCATCCGGAACGCTTCCCGGa                                                                                                                                        <  1:538138/144‑1 (MQ=255)
            aGTGGGCACCACACCGCCCGGTATGGATTGCCACCCGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCGCATCCTGGTACCCCGTCAACCGGAACGCTTCCCGGAt                                                                                                                                       >  2:173112/1‑149 (MQ=255)
            aGTGGGCACCACACCGCCCGGGATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGGGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGAt                                                                                                                                       <  1:1416311/149‑1 (MQ=255)
             gTGGGCACCACACCTCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATg                                                                                                                                      <  1:43756/149‑1 (MQ=255)
             gTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATg                                                                                                                                      <  2:697767/149‑1 (MQ=255)
              tGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCACTAATTTATTGCTCATCCTGGTACCCCGTCATCCGGCACGCTTCCCGGATGc                                                                                                                                     >  2:1395795/1‑149 (MQ=255)
               gggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCACATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCg                                                                                                                                    >  2:672776/1‑149 (MQ=255)
                ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                   <  1:844045/149‑1 (MQ=255)
                ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGGCATCCGGGAAGCGTTCCGGATgacg                                                                                                                                   <  1:559103/149‑4 (MQ=255)
                 gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCACACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATgacgg                                                                                                                                  <  2:835555/149‑5 (MQ=255)
                 gCACCACACCGCCCGCTAGGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAACAGGCATTGTTACAGCAATTCCCGAATTTAGTGCTCATCCTGGTACCCCTTCATCCGGAACGCTTCCCGGAtagatc                                                                                                                                  >  2:538138/1‑144 (MQ=255)
                  caccacACCGCCCGGGATGGAAGGCCACCACCACTCACGACGGCGAAGAGAGTGGGGGGAGCGCCGCACATCAGTCATTGTTACAGCTCTTCCCGAAGTTATTGCTCATCCGGGGACCCCGTCATCCGGAACGCTTCCCGTATGCGAtt                                                                                                                                 <  2:1050118/149‑1 (MQ=255)
                   acgacacCGCCCGGTATGGATTGCCACCAGCACACACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATATATTGCTCATCCTGGTACCCCGACATCCGGAACGCTTCCCGGATGCGATTa                                                                                                                                >  1:138172/4‑149 (MQ=255)
                    ccacACCGGCCGGTATGGATTGCCACCAGCACTCAAGTAGGCGAAGAGTGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTaa                                                                                                                               <  2:1490214/149‑1 (MQ=255)
                     cacaCCGCCCGGTATGGATTGCCACCTGCACTCACGAAGGAGAAGAGAGTGTGGTGCTCGCCGAACATCAGGCATTGTTTCAGCACTTCCCGTATTTATTGCTCAGCCTGGTACCCCGTCATCCGGAACGCTTCCGGGCTGCGAtgacc                                                                                                                              >  2:895696/1‑145 (MQ=255)
                      acacCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAcc                                                                                                                             <  1:941451/149‑1 (MQ=255)
                      acacCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAcc                                                                                                                             <  2:1269214/149‑1 (MQ=255)
                        acCTCCCGGTATGGCTTTCCACCAGCACTCAGGATTGCGAAGATAGTGTGTTGATCGCCTCACATCAGGCATTGTTACAGCAATTCCAGAATTTATTGCTCATCCTGTTACCCCTTCATCCGTAACTCTTCACGGATGCGATTAACCtt                                                                                                                           <  2:1195208/149‑1 (MQ=255)
                           gCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTTGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTc                                                                                                                        <  1:868521/149‑1 (MQ=255)
                            cccGGTATGGATTGCCACGAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTcc                                                                                                                       <  1:1429999/149‑1 (MQ=255)
                              cGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGc                                                                                                                     >  2:378632/1‑149 (MQ=255)
                                gTATGGATTGCCACCAGCACTCACGCAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCa                                                                                                                   <  1:1507559/149‑1 (MQ=255)
                                gTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCa                                                                                                                   <  2:1448382/149‑1 (MQ=255)
                                   tGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGc                                                                                                                <  2:1505049/149‑1 (MQ=255)
                                    ggATTGCCACTAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCt                                                                                                               >  2:55327/1‑149 (MQ=255)
                                     gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTAGCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                              <  2:1199479/149‑1 (MQ=255)
                                     gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                              <  2:1145214/149‑1 (MQ=255)
                                       ttGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCACGGATGCGATTAACCTTGTCCGCCAGGCTGGa                                                                                                            >  1:292658/1‑149 (MQ=255)
                                          ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTa                                                                                                         >  2:40085/1‑149 (MQ=255)
                                          ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTa                                                                                                         >  1:45208/1‑149 (MQ=255)
                                          ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGGTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTa                                                                                                         >  2:285638/1‑149 (MQ=255)
                                            accaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAg                                                                                                       >  2:489950/1‑149 (MQ=255)
                                            accaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCACGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTacg                                                                                                       >  1:117486/1‑147 (MQ=255)
                                            accaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAGTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAg                                                                                                       <  2:1346101/149‑1 (MQ=255)
                                             ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                      >  2:1251920/1‑149 (MQ=255)
                                              caGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCGCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                     <  2:272734/149‑1 (MQ=255)
                                                gCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCtat                                                                                                   <  2:1129819/149‑1 (MQ=255)
                                                 cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGGCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCtata                                                                                                  <  1:893697/149‑1 (MQ=255)
                                                   cTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATc                                                                                                <  2:659788/149‑1 (MQ=255)
                                                      ccGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCACTTCCCGAATTTATTGCTCATCCTGGTACCCCGGCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcaca                                                                                             <  2:926867/148‑1 (MQ=255)
                                                      aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcaca                                                                                             >  2:1497299/1‑149 (MQ=255)
                                                            gCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCTCACGctct                                                                                       >  1:669272/1‑149 (MQ=255)
                                                            gCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGctct                                                                                       >  1:1522204/1‑149 (MQ=255)
                                                                  agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTGCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                 <  1:1135072/149‑1 (MQ=255)
                                                                     gtgtGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTATCCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGaa                                                                              <  1:304752/149‑1 (MQ=255)
                                                                      tgtgGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAg                                                                             >  1:255006/1‑149 (MQ=255)
                                                                         ggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTcc                                                                          <  1:489950/149‑1 (MQ=255)
                                                                          gtgATCGCCGCACATGAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCTCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTccc                                                                         <  2:28638/149‑1 (MQ=255)
                                                                           tgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTcccc                                                                        <  2:367916/149‑1 (MQ=255)
                                                                            gATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGACAGGCTGGACTAAGCTATATCACACGCTCTTCAGCGGATGTCCCCt                                                                       >  1:1372813/1‑149 (MQ=255)
                                                                             aTCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGACCGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                      >  2:235056/1‑149 (MQ=255)
                                                                               cgccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTcca                                                                    >  1:111743/1‑149 (MQ=255)
                                                                                 ccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccacc                                                                  >  1:635012/1‑149 (MQ=255)
                                                                                  cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccacca                                                                 >  2:958694/1‑149 (MQ=255)
                                                                                  cgcACATCAGGCATTGTTACAGAAATTCACTAATTTATTGCTCATCCTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACTCTCTTCAGGGGAAGTCCCCTccacca                                                                 <  1:704131/149‑1 (MQ=255)
                                                                                    cacaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCCCACGCTCTTCCGGGGAAGTCCCCTCCACAAGc                                                               >  2:138943/1‑149 (MQ=255)
                                                                                    cacaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGc                                                               <  2:443587/149‑1 (MQ=255)
                                                                                     acaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCa                                                              >  1:1440923/1‑149 (MQ=255)
                                                                                       aTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTGCCCGGATGCGATTAACCTTTTCCGCCAGGCTGGACTAAGCTATATGACACGCTCTTCAGGGGAAGTCCACTCAAACAGAACg                                                            >  2:276016/1‑149 (MQ=255)
                                                                                        tCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCAAGGCTGGACTAAGCTATATCACACGCTCTTCAGCGGAAGTCCCATCCACCAGCACGc                                                           >  2:336443/1‑149 (MQ=255)
                                                                                         taGGCATTGTTACAGCAATTCCCGATTTTATTGCTCATGCTGGTACCCCGTCATCGGGAACGCTTCCTGGATGCGATTATCCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCa                                                          <  2:1052975/148‑1 (MQ=255)
                                                                                            gCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGt                                                       >  2:803000/1‑149 (MQ=255)
                                                                                             cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGtt                                                      <  1:784031/149‑1 (MQ=255)
                                                                                               ttGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTtgt                                                    >  1:1201991/1‑149 (MQ=255)
                                                                                                tGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTtgtg                                                   <  1:1102168/149‑1 (MQ=255)
                                                                                                   tACAGCAATTCACGAATTTATTTCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGtt                                                <  2:1390367/149‑1 (MQ=255)
                                                                                                    aCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCAAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTg                                               >  1:1283273/1‑149 (MQ=255)
                                                                                                      aGCAAATCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGGCCCCTCCACCAGCACGCAGGTTGTGGTTGGc                                             <  2:175170/149‑1 (MQ=255)
                                                                                                        cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGa                                           >  2:769368/1‑149 (MQ=255)
                                                                                                          aTTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATa                                         <  1:7279/149‑1 (MQ=255)
                                                                                                          aTTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATa                                         <  1:1312287/149‑1 (MQ=255)
                                                                                                          aTTCCCGAATTTATTGCTCATCCTGGTACCCCGGCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATa                                         <  2:142616/149‑1 (MQ=255)
                                                                                                              ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCCCACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGAt                                     >  2:816706/1‑149 (MQ=255)
                                                                                                              ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGAt                                     <  1:1413034/149‑1 (MQ=255)
                                                                                                              ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCACTCCTCCAGCACGCAGGTTGTGGTTGGCGATACGAt                                     >  1:123034/1‑149 (MQ=255)
                                                                                                                gAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATgg                                   >  1:1190747/1‑149 (MQ=255)
                                                                                                                 aaTTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATggg                                  >  1:57937/1‑149 (MQ=255)
                                                                                                                      aTTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGt                             >  2:1169043/1‑149 (MQ=255)
                                                                                                                         gCTGATCCTGGTACCCCGTCATCCGGAACGATGCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGTGTGTGGTTGGCGATACGATGGGCGAGTTGa                          <  2:951642/149‑1 (MQ=255)
                                                                                                                          cTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGAt                         <  1:57658/149‑1 (MQ=255)
                                                                                                                           tCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATg                        <  2:1522204/149‑1 (MQ=255)
                                                                                                                           tCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATg                        >  2:824762/1‑149 (MQ=255)
                                                                                                                           tCATCCTGGTACCCCGGCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATg                        <  2:1245965/149‑1 (MQ=255)
                                                                                                                              tCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTa                     <  1:1354843/149‑1 (MQ=255)
                                                                                                                              tCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTa                     <  1:1497299/149‑1 (MQ=255)
                                                                                                                              tCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGGCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTa                     <  1:1372799/149‑1 (MQ=255)
                                                                                                                               ccTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTAc                    >  1:744104/1‑149 (MQ=255)
                                                                                                                                  ggTACCCCCTCATCCGGAACGCTTCCCGGTTGCGATTAACGTTGTCCTCCAGGCTTGACTAAGCTATATCACACTCTCTTCATGGGACGTCCACTCCACCAGCACGCAGGTTGTGGTTGGCGTTACGATGGTCGAGTTGATGTTACTGt                 <  2:692819/149‑1 (MQ=255)
                                                                                                                                       cccGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGGTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGGATGGc            >  2:1070903/1‑149 (MQ=255)
                                                                                                                                            cATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCCCCCGCACGCAGGTTGTGGTGGGCGATACGCTGGGCGAGTTGATGTTACTGTATGGCATTGc       >  1:375803/1‑149 (MQ=255)
                                                                                                                                            cATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGc       >  2:1150715/1‑149 (MQ=255)
                                                                                                                                               ccGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGa    >  1:1543096/1‑149 (MQ=255)
                                                                                                                                                cGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGTGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAt   <  1:1104884/149‑1 (MQ=255)
                                                                                                                                                 ggAACGCTTCCCGGATGCGATTCACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCATCACGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAtc  <  2:1254818/149‑1 (MQ=255)
                                                                                                                                                 ggAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAtc  >  2:1073787/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATC  >  NC_000913/3809207‑3809500

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