Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,021,309 C→T 5.7% T74T (ACC→ACT fliN → flagellar motor switch protein FliN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,021,3090CT 5.7% 219.8 / 2.2 70T74T (ACC→ACTfliNflagellar motor switch protein FliN
Reads supporting (aligned to +/- strand):  ref base C (28/38);  new base T (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCG  >  NC_000913/2021160‑2021457
                                                                                                                                                     |                                                                                                                                                    
gAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGAcc                                                                                                                                                      >  2:650659/1‑150 (MQ=255)
gAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGAcc                                                                                                                                                      >  1:1407367/1‑150 (MQ=255)
 aaCAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCa                                                                                                                                                     <  1:795693/150‑1 (MQ=255)
 aaCAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGAGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCa                                                                                                                                                     >  1:1390745/1‑150 (MQ=255)
     aaaaTCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCaa                                                                                                                                                 >  1:452446/1‑150 (MQ=255)
         tCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAagag                                                                                                                                             <  2:363526/150‑1 (MQ=255)
             ccAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGt                                                                                                                                         <  2:452446/150‑1 (MQ=255)
               agcagcAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTg                                                                                                                                       <  1:967323/150‑1 (MQ=255)
                gcagcaAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACTATCAAAGAGCTGTTg                                                                                                                                       <  2:1167509/149‑1 (MQ=255)
                gcagcaAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACTATCAAAGAGCTGTTg                                                                                                                                       <  2:1162600/149‑1 (MQ=255)
                gcagcaAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACTATCAAAGAGCTGTTg                                                                                                                                       >  1:1167509/1‑149 (MQ=255)
                gcagcaAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACTATCAAAGAGCTGTTg                                                                                                                                       >  1:1162600/1‑149 (MQ=255)
                 cagcaAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCg                                                                                                                                     <  2:192567/150‑1 (MQ=255)
                     aaaaGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTg                                                                                                                                 >  1:1339616/1‑150 (MQ=255)
                         gcgcTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGc                                                                                                                             <  2:227971/150‑1 (MQ=255)
                             tGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAgg                                                                                                                         <  1:873857/150‑1 (MQ=255)
                             tGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAgg                                                                                                                         >  1:550382/1‑150 (MQ=255)
                              gcagagaCGGTGTTACAGAAATTTGGCGTTGGTGATGTCAGCGTAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAggg                                                                                                                        >  2:191711/4‑150 (MQ=255)
                              gCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAggg                                                                                                                        <  2:84005/150‑1 (MQ=255)
                                ggagaCTGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATGGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCTGATGACCATCAAAGAGCTGTTGCGTCAGTCGCaacagct                                                                                                                      <  2:849307/149‑6 (MQ=255)
                                  agaCGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCcg                                                                                                                    <  1:1251026/150‑1 (MQ=255)
                                        tgtTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATCGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTgg                                                                                                              >  2:770151/1‑150 (MQ=255)
                                            ccAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCt                                                                                                          <  1:918423/150‑1 (MQ=255)
                                            ccAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCt                                                                                                          >  2:1420068/1‑150 (MQ=255)
                                            ccAGCAATTTGGCGGTGGAGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGATCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGAGTCTGACGCAAGGGTCCGTCGTGGCGCt                                                                                                          >  1:798133/1‑150 (MQ=255)
                                             cAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTg                                                                                                         >  1:1038698/1‑150 (MQ=255)
                                                  aTTTGGCGGTGGTTATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTGCGGTCAATCTTACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACgg                                                                                                    <  1:703220/150‑1 (MQ=255)
                                                   tttGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGt                                                                                                   <  2:642384/150‑1 (MQ=255)
                                                    ttGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTc                                                                                                  <  1:425163/150‑1 (MQ=255)
                                                       gCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTgg                                                                                               <  2:109975/150‑1 (MQ=255)
                                                         ggtggtGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTggcg                                                                                             <  1:682965/150‑1 (MQ=255)
                                                            ggtgATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTggcgggc                                                                                          <  1:461247/150‑1 (MQ=255)
                                                            ggtgATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTggcgggc                                                                                          >  1:40627/1‑150 (MQ=255)
                                                                aTGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTAGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAAc                                                                                      >  1:1286543/1‑150 (MQ=255)
                                                                 tGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAAcc                                                                                     >  1:1354970/1‑150 (MQ=255)
                                                                  gTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCa                                                                                    >  2:144270/1‑150 (MQ=255)
                                                                     aGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGGCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGGCGTGGCGCTGGACGGTCTGACGGGCGAACCACTg                                                                                 <  1:1118838/150‑1 (MQ=255)
                                                                      gCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTgg                                                                                >  2:845227/1‑150 (MQ=255)
                                                                      gCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTgg                                                                                <  2:850928/150‑1 (MQ=255)
                                                                       cGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGa                                                                               <  1:773276/150‑1 (MQ=255)
                                                                         gAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGata                                                                             >  1:1184244/1‑150 (MQ=255)
                                                                            cGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTc                                                                          >  2:465256/1‑150 (MQ=255)
                                                                            cGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTc                                                                          <  2:1286543/150‑1 (MQ=255)
                                                                             gTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCt                                                                         >  1:878664/1‑150 (MQ=255)
                                                                               tGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGa                                                                       >  2:1398932/1‑150 (MQ=255)
                                                                                     atatCGACCTTATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACg                                                                 <  2:1285224/150‑1 (MQ=255)
                                                                                      tatCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACgg                                                                <  1:191711/150‑1 (MQ=255)
                                                                                      tatCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACgg                                                                >  1:713372/1‑150 (MQ=255)
                                                                                       atCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGt                                                               <  1:845227/150‑1 (MQ=255)
                                                                                        tCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGtt                                                              >  1:1061768/1‑150 (MQ=255)
                                                                                                 ttATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCg                                                     <  1:1420068/150‑1 (MQ=255)
                                                                                                    tGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGccc                                                  >  2:228009/1‑150 (MQ=255)
                                                                                                    tGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGccc                                                  >  2:353411/1‑150 (MQ=255)
                                                                                                         atTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGc                                             >  2:1309445/1‑150 (MQ=255)
                                                                                                         atTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGc                                             >  1:554840/1‑150 (MQ=255)
                                                                                                                   aGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAgtggt                                     >  2:150170/1‑148 (MQ=255)
                                                                                                                   aGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAgtggt                                     <  1:150170/148‑1 (MQ=255)
                                                                                                                         ccGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTg                             >  1:437225/1‑150 (MQ=255)
                                                                                                                            tCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCg                          <  2:269491/150‑1 (MQ=255)
                                                                                                                                gCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTAGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATaa                      >  1:369829/1‑150 (MQ=255)
                                                                                                                                 cTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATaaa                     <  2:1354970/150‑1 (MQ=255)
                                                                                                                                 cTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATaaa                     <  2:1038698/150‑1 (MQ=255)
                                                                                                                                      ccGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATgg                <  1:465256/150‑1 (MQ=255)
                                                                                                                                      ccGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATgg                <  2:1390745/150‑1 (MQ=255)
                                                                                                                                           cgcgGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGc           <  1:1304393/150‑1 (MQ=255)
                                                                                                                                            gcgGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCg          <  1:770151/150‑1 (MQ=255)
                                                                                                                                             cgGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCgg         <  1:1398932/150‑1 (MQ=255)
                                                                                                                                              ggATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGa        <  2:40627/150‑1 (MQ=255)
                                                                                                                                              ggATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGa        <  2:56587/150‑1 (MQ=255)
                                                                                                                                              ggATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGa        <  2:1184244/150‑1 (MQ=255)
                                                                                                                                                 tGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCa     >  2:799699/1‑150 (MQ=255)
                                                                                                                                                  gACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCAc    <  1:353411/150‑1 (MQ=255)
                                                                                                                                                    ccATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCg  >  2:971544/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACAAAAATCAACCAGCAGCAAAAGCGCTGCCGAGACGGTGTTCCAGCAATTTGGCGGTGGTGATGTCAGCGGAACGTTGCAGGATATCGACCTGATTATGGATATTCCGGTCAAGCTGACCGTCGAGCTGGGCCGTACGCGGATGACCATCAAAGAGCTGTTGCGTCTGACGCAAGGGTCCGTCGTGGCGCTGGACGGTCTGGCGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCG  >  NC_000913/2021160‑2021457

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