Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,151,636 G→A R117C (CGT→TGT)  exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,151,6360GA99.2% 490.3 / NA 122R117C (CGT→TGT) exbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (58/63);  minor base C (0/1);  total (58/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 122 total aligned reads displayed.

TAGTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAA  >  NC_000913/3151489‑3151783
                                                                                                                                                   |                                                                                                                                                   
tAGTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTc                                                                                                                                                  >  2:203616/1‑151 (MQ=255)
  gTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTcc                                                                                                                                                 >  2:663462/1‑150 (MQ=255)
  gTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTcc                                                                                                                                                 >  4:15727/1‑150 (MQ=255)
  gTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCa                                                                                                                                                <  1:902935/151‑1 (MQ=255)
   tGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTcc                                                                                                                                                 >  2:78680/1‑149 (MQ=255)
       cTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACgg                                                                                                                                           >  2:908493/1‑151 (MQ=255)
         gCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGCCCAACGACCGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACgg                                                                                                                                           <  2:11284/149‑1 (MQ=255)
              tGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTgg                                                                                                                                    >  2:980521/1‑151 (MQ=255)
                 gcgATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTAc                                                                                                                                 <  1:589128/151‑1 (MQ=255)
                    aTACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGtt                                                                                                                              <  1:1012128/151‑1 (MQ=255)
                    aTACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGtt                                                                                                                              >  2:1023868/1‑151 (MQ=255)
                    aTACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGtt                                                                                                                              <  2:87589/151‑1 (MQ=255)
                      aCCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCt                                                                                                                            <  2:948535/151‑1 (MQ=255)
                       ccAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCtt                                                                                                                           >  1:248410/1‑151 (MQ=255)
                        cAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCttt                                                                                                                          <  2:387432/151‑1 (MQ=255)
                         aaTAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCttt                                                                                                                          >  2:1193998/1‑150 (MQ=255)
                          aTAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTaa                                                                                                                        <  3:318604/151‑1 (MQ=255)
                            aaaGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAAt                                                                                                                       >  3:375338/1‑150 (MQ=255)
                             aaGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATAc                                                                                                                     >  2:1270999/1‑151 (MQ=255)
                                  gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTc                                                                                                                 <  1:276738/150‑1 (MQ=255)
                                  gTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCg                                                                                                                >  2:379304/1‑151 (MQ=255)
                                   ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGt                                                                                                               >  1:421251/1‑151 (MQ=255)
                                   ttCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGt                                                                                                               >  1:560685/1‑151 (MQ=255)
                                    tCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTGCCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGt                                                                                                               <  2:828881/150‑1 (MQ=255)
                                        gatgCCCCGGGCCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTTCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTAtc                                                                                                           <  2:1123467/150‑1 (MQ=255)
                                           gCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTAtcgtc                                                                                                        <  2:412983/150‑1 (MQ=255)
                                            ccccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGc                                                                                                      >  2:88022/1‑151 (MQ=255)
                                            ccccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGc                                                                                                      >  3:102829/1‑151 (MQ=255)
                                            ccccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGc                                                                                                      <  1:1023529/151‑1 (MQ=255)
                                             cccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCt                                                                                                     >  3:46459/1‑151 (MQ=255)
                                             cccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCt                                                                                                     <  2:1152843/151‑1 (MQ=255)
                                              ccAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTg                                                                                                    <  4:194687/151‑1 (MQ=255)
                                                aGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGcc                                                                                                  >  2:457360/1‑151 (MQ=255)
                                                    cGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCt                                                                                              <  1:825587/151‑1 (MQ=255)
                                                     gTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTg                                                                                             <  2:493082/151‑1 (MQ=255)
                                                     gTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCAATGCGGCGCCCCGACCCTCCCGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTg                                                                                             <  2:1049785/151‑1 (MQ=255)
                                                      tACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTg                                                                                             >  2:63673/1‑150 (MQ=255)
                                                          aaaCAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGt                                                                                        >  2:5504/1‑151 (MQ=255)
                                                           aaCAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTCCCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGtt                                                                                       <  3:227125/151‑1 (MQ=255)
                                                           aaCAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGtt                                                                                       >  1:1009611/1‑151 (MQ=255)
                                                           aaCAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCGTCTGACAGtt                                                                                       <  4:335347/151‑1 (MQ=255)
                                                             cAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTc                                                                                      <  2:577893/150‑1 (MQ=255)
                                                              aGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCa                                                                                    <  1:1172873/151‑1 (MQ=255)
                                                                   aacgaacgGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTc                                                                                >  1:642801/1‑150 (MQ=255)
                                                                     cgaacgGAGAAATCGCGCCGATGGTAGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCAtt                                                                             >  3:340935/1‑151 (MQ=255)
                                                                       aacgGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCt                                                                           <  1:312467/151‑1 (MQ=255)
                                                                        acgGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTg                                                                          >  4:174673/1‑151 (MQ=255)
                                                                         cgGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgc                                                                         <  3:68209/151‑1 (MQ=255)
                                                                          ggAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgcg                                                                        >  2:1058851/1‑151 (MQ=255)
                                                                          ggAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgcg                                                                        >  4:246952/1‑151 (MQ=255)
                                                                           gagaAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgcgc                                                                       <  1:1187842/151‑1 (MQ=255)
                                                                            agaAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTgcgat                                                                      <  1:522365/151‑3 (MQ=255)
                                                                            agaAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCt                                                                      >  1:901285/1‑151 (MQ=255)
                                                                               aaTCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCa                                                                   <  2:954003/151‑1 (MQ=255)
                                                                                 tCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGAAGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCAtt                                                                 >  2:503968/1‑151 (MQ=255)
                                                                                  cgcgCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGGCCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTg                                                                <  4:123112/151‑1 (MQ=255)
                                                                                   gcgcCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTg                                                                >  2:833988/1‑150 (MQ=255)
                                                                                   gcgcCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTg                                                                >  4:412666/1‑150 (MQ=255)
                                                                                    cgcCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAg                                                              <  2:790231/151‑1 (MQ=255)
                                                                                      ccGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCa                                                            <  3:389518/151‑1 (MQ=255)
                                                                                       cGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCa                                                            >  2:770129/1‑150 (MQ=255)
                                                                                       cGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCa                                                            >  3:161557/1‑150 (MQ=255)
                                                                                        gATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCaaa                                                          <  3:429381/151‑1 (MQ=255)
                                                                                         aTGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAAt                                                         >  2:438116/1‑151 (MQ=255)
                                                                                          tGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATg                                                        >  4:439291/1‑151 (MQ=255)
                                                                                          tGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTGCCAGCTCATTCTGCGCTTCATTGAGCAAATg                                                        >  1:1144262/1‑151 (MQ=255)
                                                                                           ggTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATg                                                        <  2:129310/150‑1 (MQ=255)
                                                                                            gTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCa                                                      <  1:649413/151‑1 (MQ=255)
                                                                                               gCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAg                                                     >  1:1044753/1‑149 (MQ=255)
                                                                                                ccAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCt                                                  <  3:47826/151‑1 (MQ=255)
                                                                                                   gATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTaa                                               <  1:1234529/151‑1 (MQ=255)
                                                                                                    aTAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTa                                                >  1:544881/1‑149 (MQ=255)
                                                                                                    aTAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAg                                              <  1:1068039/151‑1 (MQ=255)
                                                                                                      aGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCt                                            <  4:428856/151‑1 (MQ=255)
                                                                                                       gCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCtt                                           <  2:1116941/151‑1 (MQ=255)
                                                                                                        ccGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCttt                                          >  2:671316/1‑151 (MQ=255)
                                                                                                         cGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCtttt                                         <  4:164552/151‑1 (MQ=255)
                                                                                                            tACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCt                                      >  2:470245/1‑151 (MQ=255)
                                                                                                             aCCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCt                                      <  2:756002/150‑1 (MQ=255)
                                                                                                               cgcgGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTAcc                                   <  3:207659/151‑1 (MQ=255)
                                                                                                                gcgGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTAcc                                   >  1:1206725/1‑150 (MQ=255)
                                                                                                                 cgGCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCa                                  >  4:357475/1‑150 (MQ=255)
                                                                                                                  ggCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCaaa                                <  1:229373/151‑1 (MQ=255)
                                                                                                                   gCCCATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCaaa                                >  2:1082702/1‑150 (MQ=255)
                                                                                                                     ccATTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATc                             >  2:37449/1‑151 (MQ=255)
                                                                                                                       aTTTGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAg                           <  1:1271231/151‑1 (MQ=255)
                                                                                                                        tttGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAgc                          >  3:185658/1‑151 (MQ=255)
                                                                                                                         ttGACGACCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAgc                          <  4:228573/150‑1 (MQ=255)
                                                                                                                              gacCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGat                    >  4:158191/1‑151 (MQ=255)
                                                                                                                              gacCCACTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGat                    <  2:1063169/151‑1 (MQ=255)
                                                                                                                                   aCTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGt               >  1:489368/1‑151 (MQ=255)
                                                                                                                                   aCTGCGGCGACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGt               >  3:189444/1‑151 (MQ=255)
                                                                                                                                     tGCGGCGACCCGACACTCCAGACGGAAGCTGGCACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGtt              >  1:1198594/1‑150 (MQ=255)
                                                                                                                                      gcggcgACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTgg            <  1:640131/151‑1 (MQ=255)
                                                                                                                                       cggcgACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGc           <  1:996981/151‑1 (MQ=255)
                                                                                                                                         gcgACCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCt         <  1:938393/151‑1 (MQ=255)
                                                                                                                                            aCCCGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAGGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGt      <  1:663595/151‑1 (MQ=255)
                                                                                                                                             cccGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGtt     >  1:1224143/1‑151 (MQ=255)
                                                                                                                                               cGACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATAGTTGGCCTGGttt    >  1:192320/1‑150 (MQ=255)
                                                                                                                                                gACACTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTaa  <  4:450932/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAGTGGTCTGCGTTTGCGCGATACCAATAAAGCTGTTCATGATGCCCCAGACCGTACCAAACAGACCAACGAACGGAGAAATCGCGCCGATGGTTGCCAGATAGCCGTTACCGCGGCCCATTTGACGACCCACTGCGGCGACCCGACGCTCCAGACGGAAGCTGGTACGTTCTTTAATACCTTCGTTATCGTCGCTGCCTTCTGACAGTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAA  >  NC_000913/3151489‑3151783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: