Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 779,661 T→G 5.2% W22G (TGG→GGG)  ybgF → periplasmic TolA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913779,6610TG 5.2% 348.1 / 6.0 118W22G (TGG→GGG) ybgFperiplasmic TolA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (59/52);  new base G (3/4);  total (62/56)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 118 total aligned reads displayed.

TAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATT  >  NC_000913/779513‑779807
                                                                                                                                                    |                                                                                                                                                  
tAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTg                                                                                                                                                   <  1:1355597/150‑1 (MQ=255)
   aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                >  1:1051042/1‑150 (MQ=255)
   aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                >  2:454721/1‑150 (MQ=255)
   aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                >  1:1430298/1‑150 (MQ=255)
   aGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGc                                                                                                                                                >  1:543077/1‑150 (MQ=255)
     aaaaaCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCTACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCg                                                                                                                                              >  1:315189/1‑150 (MQ=255)
      aaaaCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGc                                                                                                                                             >  2:1322030/1‑150 (MQ=255)
         aCCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtat                                                                                                                                          <  1:793446/150‑3 (MQ=255)
          ccTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCGTGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCAGGGGGCCGCtata                                                                                                                                         <  2:664166/150‑4 (MQ=255)
           cTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTACAGCGGCCCAGGGGGCCGCtatac                                                                                                                                        <  1:384723/150‑5 (MQ=255)
                 tACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAgg                                                                                                                                  <  1:378322/150‑1 (MQ=255)
                  aCTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGc                                                                                                                                 >  2:1138074/1‑150 (MQ=255)
                     ggggCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCAcc                                                                                                                              <  1:19751/150‑1 (MQ=255)
                     gggTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCAcc                                                                                                                              >  1:617623/1‑150 (MQ=255)
                     gggTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCAcc                                                                                                                              >  1:71082/1‑150 (MQ=255)
                     gggTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCTGGGCCGCTTTTGCTCAGGCACCa                                                                                                                             >  2:5191/1‑150 (MQ=255)
                     gggTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGGCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCAcc                                                                                                                              <  1:13633/150‑1 (MQ=255)
                         cATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATATCGGCCCCCTGGGCCGCTTTTGCTCAGTCattacca                                                                                                                          >  2:913863/1‑144 (MQ=255)
                         cATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCattacca                                                                                                                          >  2:919963/1‑144 (MQ=255)
                         cATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATc                                                                                                                          <  2:1424699/150‑1 (MQ=255)
                           tGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCag                                                                                                                        <  2:13057/150‑1 (MQ=255)
                            gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagt                                                                                                                       <  1:1181455/150‑1 (MQ=255)
                            gACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagt                                                                                                                       >  1:1168732/1‑150 (MQ=255)
                             aCGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagta                                                                                                                      <  1:1644150/150‑1 (MQ=255)
                               gAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCagtagt                                                                                                                    >  2:1540550/1‑150 (MQ=255)
                                aaGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCATTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCCACAATCAGTAgtg                                                                                                                   >  1:198921/1‑150 (MQ=255)
                                     gCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGc                                                                                                              >  2:1392111/1‑150 (MQ=255)
                                     gCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCGCAGGCCCCAATCAGTAGTGTCGGc                                                                                                              >  1:877415/1‑150 (MQ=255)
                                      cATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCt                                                                                                             >  2:1340979/1‑150 (MQ=255)
                                        tACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCa                                                                                                           <  1:207108/150‑1 (MQ=255)
                                         aCTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAg                                                                                                          >  1:717543/1‑150 (MQ=255)
                                           tCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCGCAGGc                                                                                                        >  2:468501/1‑150 (MQ=255)
                                            ccAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGcg                                                                                                       >  1:1281150/1‑149 (MQ=255)
                                            ccAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCGCAGGcg                                                                                                       >  2:477898/1‑149 (MQ=255)
                                              aaaaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCg                                                                                                     >  2:513474/1‑150 (MQ=255)
                                                aaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGt                                                                                                   <  2:273075/150‑1 (MQ=255)
                                                aaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGt                                                                                                   <  1:112513/150‑1 (MQ=255)
                                                aaaCCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGt                                                                                                   <  2:1403821/150‑1 (MQ=255)
                                                   ccGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGa                                                                                                <  2:1261229/150‑1 (MQ=255)
                                                     gtcgtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAg                                                                                              >  2:293818/1‑150 (MQ=255)
                                                        gtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGAcc                                                                                           <  1:1623481/150‑1 (MQ=255)
                                                        gtGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGAcc                                                                                           >  1:820129/1‑150 (MQ=255)
                                                          gCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACcgc                                                                                         <  2:1158930/150‑1 (MQ=255)
                                                             gTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTc                                                                                      >  2:1246266/1‑150 (MQ=255)
                                                              tACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCa                                                                                     >  2:1003457/1‑150 (MQ=255)
                                                               aCTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCAc                                                                                    <  2:464106/150‑1 (MQ=255)
                                                                cTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGCCCGCGTCACt                                                                                   >  2:1606422/1‑150 (MQ=255)
                                                                 tGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCAGAGGCTCGGGCGAAGACCGCGTCAcgc                                                                                  >  2:1662044/1‑148 (MQ=255)
                                                                  ggTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCCCCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCCCTCa                                                                                 >  2:959399/1‑150 (MQ=255)
                                                                   gTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTTTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCaa                                                                                >  1:1080934/1‑150 (MQ=255)
                                                                    tttACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAAc                                                                               <  1:1278396/150‑1 (MQ=255)
                                                                         tAAGAGAATTGCTTGAGCAGTAACTTCAGACATAAACTATTGAGGCTGTCGTTACTGGGTGGTATTGCGGCCCCCGTGGCCGCTTTTGCTCAGGCACCATTCAGTAGTGTCGGCTCAGTCTCGGTCGAAGACCGCGTCACTCAACTTGAg                                                                          <  1:1059334/150‑1 (MQ=255)
                                                                         tAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAg                                                                          >  2:109540/1‑150 (MQ=255)
                                                                          aaGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         >  2:526898/1‑150 (MQ=255)
                                                                          aaGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         <  2:1069039/150‑1 (MQ=255)
                                                                          aaGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGc                                                                         >  1:209846/1‑150 (MQ=255)
                                                                          aaGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCACAGGCACAAATCAGTAGTGTCGGCTCAGGCTCGGACGAATACCGCGTCACTCAACTTGAgg                                                                         >  2:178224/1‑149 (MQ=255)
                                                                                  tGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCt                                                                 <  2:1427632/150‑1 (MQ=255)
                                                                                   gCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTa                                                                >  1:448040/1‑150 (MQ=255)
                                                                                    cATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGCCTCAGGCTCGGTCGAAGCCCGCGTCACTCAACTTGAGCGTATTTCTaa                                                               >  2:1292949/1‑150 (MQ=255)
                                                                                       gAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCGCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGc                                                            <  1:1589753/150‑1 (MQ=255)
                                                                                          cAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTca                                                         <  1:1392111/150‑1 (MQ=255)
                                                                                           aGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTcac                                                        <  1:1138074/150‑1 (MQ=255)
                                                                                            gTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTTGTGTCGGCTCAGGCTCGGTAGAAGACCGCGTAAATCAACTTGAGCGGATTTGTAATGCTaacc                                                       >  1:1101462/1‑149 (MQ=255)
                                                                                            gTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGAACGCGTCACGGAACTTGAGCGTTTTTCTAATGCTcacc                                                       >  1:859944/1‑149 (MQ=255)
                                                                                              aaCTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagc                                                     <  1:1246266/150‑1 (MQ=255)
                                                                                               aCTTCAGACATCAACTATTGAGTCTGTCGTTACTTGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCGCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcc                                                    <  1:278688/150‑1 (MQ=255)
                                                                                                cTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCAcagcca                                                   <  1:1198956/150‑1 (MQ=255)
                                                                                                    aGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCtt                                               <  2:1020969/150‑1 (MQ=255)
                                                                                                     gACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCttt                                              <  2:67052/150‑1 (MQ=255)
                                                                                                       cATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGTCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTa                                            >  1:198151/1‑150 (MQ=255)
                                                                                                       cATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTa                                            >  1:775705/1‑150 (MQ=255)
                                                                                                       cATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCACGGTCGAAGACCGCGTGACTCAACTTGAGCGTATTTGTAATGCTCACAGCCACCttttt                                            >  2:768643/1‑149 (MQ=255)
                                                                                                           aaCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAAccc                                        <  2:1430298/150‑1 (MQ=255)
                                                                                                            aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCa                                       <  2:543077/150‑1 (MQ=255)
                                                                                                            aCTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCa                                       <  2:1051042/150‑1 (MQ=255)
                                                                                                               aTTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACt                                    <  1:85724/150‑1 (MQ=255)
                                                                                                               aTTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTc                                   <  2:874499/150‑1 (MQ=255)
                                                                                                                 tGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTcc                                  <  2:1673275/150‑1 (MQ=255)
                                                                                                                   aGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAg                                <  2:771901/150‑1 (MQ=255)
                                                                                                                      cTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaa                             <  2:1469222/150‑1 (MQ=255)
                                                                                                                      cTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACGTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaa                             >  1:910895/1‑150 (MQ=255)
                                                                                                                        gTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGcaaca                           <  2:496720/150‑1 (MQ=255)
                                                                                                                         tCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTATAATGCTCACAGCCAGCTTTTAACCCAACTCACGcaccaa                          >  2:1093065/1‑150 (MQ=255)
                                                                                                                           gTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACt                        <  1:1340979/150‑1 (MQ=255)
                                                                                                                            ttACTGGTTGGTATAGCGGCCCCCTTGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                       <  1:1691503/150‑1 (MQ=255)
                                                                                                                            ttACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACtt                       <  2:717543/150‑1 (MQ=255)
                                                                                                                               cTGGGTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCt                    <  2:1576977/150‑1 (MQ=255)
                                                                                                                                 ggttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGa                  <  1:1322030/150‑1 (MQ=255)
                                                                                                                                 ggttggtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGa                  >  1:1195332/1‑150 (MQ=255)
                                                                                                                                     gcaaaaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAAATTTCTGATaat              >  2:484161/6‑150 (MQ=255)
                                                                                                                                      gtATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATaatc             >  2:966051/1‑150 (MQ=255)
                                                                                                                                        ataGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGGTCACAGCCCGCTTTTAACCCAACTCCAGCAGCAACGTTCGGATaatcca           >  2:320467/1‑148 (MQ=255)
                                                                                                                                        aaaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGACACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAATCCAACTCCAGCAACAACTTTCTGATaatcaa           >  1:552087/3‑150 (MQ=255)
                                                                                                                                        aaaGCGGCCCAGGGGGCCGCTTTTGCTCAGGCACCAAGCAGTAGGGTCGGCGCAGGCTCGGTCGAAGCCCGCGTCACTCAACTTGAGCGTAGTTCTAATGCTCACAGCCAGCTTGTAACCCCCACGCAGCACCAACTTTCTGATattcaa           >  1:301320/3‑150 (MQ=255)
                                                                                                                                            cGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCg       <  2:209846/150‑1 (MQ=255)
                                                                                                                                             ggCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCGAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGa      >  1:777701/1‑150 (MQ=255)
                                                                                                                                              gCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGat     <  2:745492/150‑1 (MQ=255)
                                                                                                                                              gCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGat     <  2:1160605/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATAtt  <  2:454052/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAAAGAAAAACCTGCAGTACTGGGTCATGACGAAGCGGCATACTCCAAAAACCGTCGTGCGGTACTGGTTTACTAAGAGAATTGCATGAGCAGTAACTTCAGACATCAACTATTGAGTCTGTCGTTACTGGTTGGTATAGCGGCCCCCTGGGCCGCTTTTGCTCAGGCACCAATCAGTAGTGTCGGCTCAGGCTCGGTCGAAGACCGCGTCACTCAACTTGAGCGTATTTCTAATGCTCACAGCCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATT  >  NC_000913/779513‑779807

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