Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,133,864 A→T 7.6% V21V (GTA→GTT flgF → flagellar basal‑body rod protein FlgF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,8640AT 7.6% 318.3 / 11.0 104V21V (GTA→GTTflgFflagellar basal‑body rod protein FlgF
Reads supporting (aligned to +/- strand):  ref base A (56/40);  new base T (6/2);  total (62/42)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAA  >  NC_000913/1133716‑1134011
                                                                                                                                                    |                                                                                                                                                    
tCAGTCTAACGCCCAGACCATCAAAACCCAGGCCCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGCATCAACAGGCGGTaa                                                                                                                                                     >  1:667795/1‑150 (MQ=255)
         cGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTa                                                                                                                                            >  2:1512880/1‑150 (MQ=255)
               gACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTgg                                                                                                                                      <  1:694828/150‑1 (MQ=255)
                aCCATCAAAACACAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACTGGATAGCTCAATGGATCACGCAATTTATACAGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGAGCTAACCGCCAGTAATCTGGc                                                                                                                                     >  2:1375680/1‑150 (MQ=255)
                 ccATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGTGGTAACCGCCAGTAATCTGGcc                                                                                                                                    >  1:122855/1‑150 (MQ=255)
                 ccATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGcc                                                                                                                                    >  1:145205/1‑150 (MQ=255)
                 ccATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGcc                                                                                                                                    >  1:444844/1‑150 (MQ=255)
                   aTCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCaa                                                                                                                                  <  1:249186/150‑1 (MQ=255)
                    tCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCATTGGATGAATCTATTTATACCGCGATGGGAGCATCAAGATAGTCACTGAATGAACTTGAGGAATCAGGCAGTATTCTggataat                                                                                                                                 <  2:321059/150‑6 (MQ=255)
                    tCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAAt                                                                                                                                 >  1:605769/1‑150 (MQ=255)
                      aaaaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGc                                                                                                                               >  1:1234697/1‑150 (MQ=255)
                      aaaaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGc                                                                                                                               >  2:764798/1‑150 (MQ=255)
                         aCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                                            <  1:1054873/150‑1 (MQ=255)
                         aCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                                            <  1:606037/150‑1 (MQ=255)
                           ccAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCaa                                                                                                                          >  1:708985/1‑150 (MQ=255)
                              ggACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGc                                                                                                                       >  2:1517847/1‑150 (MQ=255)
                              ggACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGc                                                                                                                       >  1:1015807/1‑150 (MQ=255)
                               gACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGcc                                                                                                                      >  2:63265/1‑150 (MQ=255)
                                aCCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGccc                                                                                                                     <  2:15150/150‑1 (MQ=255)
                                  cAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCgg                                                                                                                   <  1:53970/150‑1 (MQ=255)
                                   aGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGt                                                                                                                  >  2:1266546/1‑150 (MQ=255)
                                   aGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGt                                                                                                                  >  2:1317094/1‑150 (MQ=255)
                                    gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGtt                                                                                                                 >  2:1214951/1‑150 (MQ=255)
                                    gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGtt                                                                                                                 >  1:407382/1‑150 (MQ=255)
                                     aTCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGAGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGttt                                                                                                                >  1:935209/1‑150 (MQ=255)
                                        cTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcg                                                                                                             <  1:422235/150‑1 (MQ=255)
                                         tCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgc                                                                                                            >  1:47054/1‑150 (MQ=255)
                                            acacGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgcg                                                                                                         >  2:1134096/1‑150 (MQ=255)
                                               cGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGcc                                                                                                                              >  1:1339447/1‑126 (MQ=255)
                                               cGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGcc                                                                                                                              <  2:1339447/126‑1 (MQ=255)
                                               cGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAg                                                                                                      >  2:1348937/1‑150 (MQ=255)
                                                   ggTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGa                                                                                                  <  2:667795/150‑1 (MQ=255)
                                                     ttAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAAt                                                                                                >  1:251092/1‑150 (MQ=255)
                                                        aCTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCt                                                                                             >  1:528077/1‑150 (MQ=255)
                                                          ttACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                                           <  2:528077/150‑1 (MQ=255)
                                                          ttACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                                           <  1:1134096/150‑1 (MQ=255)
                                                            aCGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAc                                                                                         >  1:81020/1‑150 (MQ=255)
                                                             cGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcg                                                                                        <  1:1317094/150‑1 (MQ=255)
                                                             cGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcg                                                                                        >  2:1084954/1‑150 (MQ=255)
                                                              gCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCTGTTTTCGCGCGCAGTTGAATGCTTTAcgc                                                                                       <  2:586297/150‑1 (MQ=255)
                                                               cTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                                      >  2:26107/1‑150 (MQ=255)
                                                                   tCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGt                                                                                  <  1:1544676/150‑1 (MQ=255)
                                                                      cTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGcc                                                                               >  2:1264156/1‑150 (MQ=255)
                                                                       tGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCa                                                                              <  1:516436/150‑1 (MQ=255)
                                                                       tGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCa                                                                              <  2:1020254/150‑1 (MQ=255)
                                                                          cGGGATAGCTCTATGGAGCACGCAATTTATTCCGCGATGGGAGCAGCCAGCCAGACTCTGAATCTAGTGTCGGTAACCGCCAGTAATCTGGCCTATGCCTCTACGCCCGGTTTTCGCGCGCAGTTGATTGCTTTACGCGCGGTGCCAGTg                                                                           <  2:476698/150‑1 (MQ=255)
                                                                              aTAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAg                                                                       <  1:764798/150‑1 (MQ=255)
                                                                                 gCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGc                                                                    >  1:748432/1‑150 (MQ=255)
                                                                                 gCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGc                                                                    >  1:929423/1‑150 (MQ=255)
                                                                                    cAATTGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCttt                                                                 <  2:1397499/150‑1 (MQ=255)
                                                                                        ggATCACGAAAATTATACCGCAATCGGATAAGCCAGACAAACACTGAATCAACACCCGCTAACAACAAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTctct                                                             >  2:1541402/1‑150 (MQ=255)
                                                                                          aTCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGc                                                           >  1:1401463/1‑150 (MQ=255)
                                                                                           tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGcc                                                          <  1:1512880/150‑1 (MQ=255)
                                                                                            cACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGccc                                                         >  2:111165/1‑150 (MQ=255)
                                                                                            cACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGACTTTCTCTGccc                                                         >  1:293317/1‑150 (MQ=255)
                                                                                              cGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAc                                                       >  2:203858/1‑150 (MQ=255)
                                                                                              cGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAc                                                       >  1:1290767/1‑150 (MQ=255)
                                                                                                cAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                                     <  2:708985/150‑1 (MQ=255)
                                                                                                  aTTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACgcgc                                                   >  1:1053461/1‑150 (MQ=255)
                                                                                                   tttATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCa                                                  <  1:1517847/150‑1 (MQ=255)
                                                                                                    ttATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCAc                                                 <  1:492248/150‑1 (MQ=255)
                                                                                                       taCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                                              >  1:368774/1‑150 (MQ=255)
                                                                                                       taCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                                              <  1:734365/150‑1 (MQ=255)
                                                                                                        aCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTg                                             <  2:1015807/150‑1 (MQ=255)
                                                                                                        aCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTg                                             <  1:978750/150‑1 (MQ=255)
                                                                                                         ccGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                                            >  1:1402905/1‑150 (MQ=255)
                                                                                                          cgcgATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGTCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTTGAAGGGCTTTCTCTGCCCACGCGCACGTTGGt                                           <  2:758643/150‑1 (MQ=255)
                                                                                                                gggAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGc                                     <  1:26107/150‑1 (MQ=255)
                                                                                                                 ggAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTg                                                            >  2:1196461/1‑126 (MQ=255)
                                                                                                                 ggAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTg                                                            <  1:1196461/126‑1 (MQ=255)
                                                                                                                 ggAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCg                                    <  2:1129474/150‑1 (MQ=255)
                                                                                                                       gccagccaGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAcc                              >  1:123913/1‑149 (MQ=255)
                                                                                                                       gccagccaGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACg                              >  1:776858/1‑150 (MQ=255)
                                                                                                                        ccagccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGc                             >  2:860555/1‑150 (MQ=255)
                                                                                                                          agccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCAGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCg                           >  2:1456410/1‑150 (MQ=255)
                                                                                                                             cagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgc                        >  1:675298/1‑150 (MQ=255)
                                                                                                                              agACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcg                       <  1:938371/150‑1 (MQ=255)
                                                                                                                               gACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGATTGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcgc                     >  1:1345284/1‑150 (MQ=255)
                                                                                                                                 cacTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAg                   >  2:67144/1‑150 (MQ=255)
                                                                                                                                 cacTGAATAAACAGGCGGTAACCGCCAGTAATCTGGCCAATGACTCAACGCCCGGATTTCGCGCGAAGTTGAATGCTTTACGAGCGGTCCCAGTGGAAGGGCTTTCTCTGCCAAAGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAg                   >  2:1469427/1‑150 (MQ=255)
                                                                                                                                   cTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGat                 >  1:166566/1‑150 (MQ=255)
                                                                                                                                    tGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGata                >  2:732001/1‑150 (MQ=255)
                                                                                                                                     gAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatat               >  1:1027356/1‑150 (MQ=255)
                                                                                                                                     gAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                                            >  2:325539/1‑122 (MQ=255)
                                                                                                                                     gAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                                            <  1:325539/122‑1 (MQ=255)
                                                                                                                                       gattactggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGa             >  2:1459716/8‑150 (MQ=255)
                                                                                                                                        attactggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGAc            >  1:506382/7‑150 (MQ=255)
                                                                                                                                         cAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACg           >  2:638478/1‑150 (MQ=255)
                                                                                                                                         cAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACg           >  2:1219772/1‑150 (MQ=255)
                                                                                                                                          tactggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAATGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatat               <  1:568398/141‑1 (MQ=255)
                                                                                                                                          tactggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAATGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatat               >  2:568398/5‑145 (MQ=255)
                                                                                                                                          tactggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAATGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatat               <  1:1178498/141‑1 (MQ=255)
                                                                                                                                           actggCGGTTACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGcc         >  1:866035/4‑150 (MQ=255)
                                                                                                                                           actggCGGTTACCGCCAGAAATCTGGCCAATGCCTCAACGCACGCTTTTCGCGCGCAGTTGAATCCTTTACGCGAGGTCCCAGTGGAAGCGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCACATATGACGcc         >  2:509177/4‑150 (MQ=255)
                                                                                                                                           aCAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGACTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGcc         <  2:293317/150‑1 (MQ=255)
                                                                                                                                             aGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCg       <  2:1323948/150‑1 (MQ=255)
                                                                                                                                             aGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTAACGGCGTCAACGCCGGGCGCAGATATGACGCCCg       <  1:1375680/150‑1 (MQ=255)
                                                                                                                                               gCGGTAATCTACAGTAATCTTTTCAATGACTCAACGCTCGGTTTTCGCTCGCAGTTGATTGGTTTACGCGCGGTGGCAGTGTAAGGGATTTCTCTTCCGACGCGCACGTTGGTCACGGCTTCAACGCCGGGCGTCATATCTGCGCCCGGc     <  2:145205/150‑1 (MQ=255)
                                                                                                                                               gCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGc     <  2:675298/150‑1 (MQ=255)
                                                                                                                                                  gTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCaaa  >  2:480879/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAA  >  NC_000913/1133716‑1134011

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