Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 994,630 C→T 7.8% R68Q (CGG→CAG)  ssuC ← aliphatic sulfonate ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913994,6300CT 7.8% 230.8 / 6.8 77R68Q (CGG→CAG) ssuCaliphatic sulfonate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (3/3);  total (39/38)
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.47e-01

CACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCC  >  NC_000913/994485‑994776
                                                                                                                                                 |                                                                                                                                                  
cACCAGCGGGATCAGCCCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCATCAGCCGTTCTCCCCAGCGATATAGCCAGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCa                                                                                                                                                >  2:847085/1‑150 (MQ=255)
 aCCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATACAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATAGAACAGCCAATCAGCGCCCGCCAg                                                                                                                                               >  2:1069865/1‑150 (MQ=255)
     gCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGc                                                                                                                                           >  2:1062175/1‑150 (MQ=255)
      cGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCt                                                                                                                                          >  2:1513682/1‑150 (MQ=255)
         gATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAt                                                                                                                                       <  2:1772150/150‑1 (MQ=255)
            cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGc                                                                                                                                    >  2:1172282/1‑150 (MQ=255)
             aGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGcc                                                                                                                                   >  1:1085770/1‑150 (MQ=255)
             aGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGcc                                                                                                                                   >  1:613146/1‑150 (MQ=255)
             aGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGcc                                                                                                                                   <  1:513657/150‑1 (MQ=255)
               cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAACGCCAATCAGCTCCTGGCGGGCGCTGAttggcgt                                                                                                                                 <  1:1574939/150‑7 (MQ=255)
                gcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGa                                                                                                                                >  1:1668199/1‑150 (MQ=255)
                         ggCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCcac                                                                                                                       >  2:1456087/1‑150 (MQ=255)
                           cACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAg                                                                                                                     <  2:326144/150‑1 (MQ=255)
                                tGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAgg                                                                                                                                              >  2:1926841/1‑120 (MQ=255)
                                tGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAgg                                                                                                                                              <  1:1926841/120‑1 (MQ=255)
                                  cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccg                                                                                                              <  1:1315190/150‑1 (MQ=255)
                                             aaTTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGagag                                                                                                   <  1:886463/150‑1 (MQ=255)
                                              aTTGACGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGc                                                                                                  >  2:1756972/1‑150 (MQ=255)
                                                  aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                              <  2:1419579/150‑1 (MQ=255)
                                                  aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGAGCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                              <  2:1256167/150‑1 (MQ=255)
                                                     tATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGa                                                                                           <  1:1713606/150‑1 (MQ=255)
                                                        ccAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACg                                                                                        >  1:760609/1‑150 (MQ=255)
                                                           gcagcCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCg                                                                                     >  1:1608376/1‑150 (MQ=255)
                                                             agcCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTc                                                                                   <  2:963754/150‑1 (MQ=255)
                                                               ccGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTcac                                                                                 >  1:91731/1‑150 (MQ=255)
                                                                cGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTcacc                                                                                <  1:253770/150‑1 (MQ=255)
                                                                cGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCACTGGCGGAGAGCGTCCAGAACGCCGTcacc                                                                                <  2:843191/150‑1 (MQ=255)
                                                                   tctcCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACt                                                                             <  2:85044/150‑1 (MQ=255)
                                                                      ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                          >  1:1251885/1‑150 (MQ=255)
                                                                      ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                          >  2:1744416/1‑150 (MQ=255)
                                                                      ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                          >  1:1924234/1‑150 (MQ=255)
                                                                         cAGCGAGATAGCCCGCTAATCAGGCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCATTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGGCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTcc                                                                       <  2:1973749/150‑1 (MQ=255)
                                                                          aGCGAGATAGTCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCg                                                                      <  1:1164568/150‑1 (MQ=255)
                                                                          aGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCg                                                                      >  1:1342247/1‑150 (MQ=255)
                                                                             gagaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTg                                                                   <  1:1308711/150‑1 (MQ=255)
                                                                              agaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGa                                                                  <  1:799060/150‑1 (MQ=255)
                                                                               taTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                                 <  1:847085/149‑1 (MQ=255)
                                                                                   gCCCGCTAATCAGTCCCAAAATCAACCACAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCCCCAGGAGCTGAACGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAgccc                                                             >  1:1291133/1‑147 (MQ=255)
                                                                                     ccGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCaaa                                                           <  2:1139479/150‑1 (MQ=255)
                                                                                       gCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAAt                                                         <  1:953164/150‑1 (MQ=255)
                                                                                          aaTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                                 >  2:966830/1‑139 (MQ=255)
                                                                                          aaTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                                 <  1:966830/139‑1 (MQ=255)
                                                                                          aaTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATAcg                                                      <  2:1739161/150‑1 (MQ=255)
                                                                                               gTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCg                                                 >  2:1152193/1‑150 (MQ=255)
                                                                                               gTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCg                                                 >  2:1663147/1‑150 (MQ=255)
                                                                                                tCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGa                                                <  1:1456087/150‑1 (MQ=255)
                                                                                                  ccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAca                                              >  1:1072164/1‑150 (MQ=255)
                                                                                                  ccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAca                                              >  2:1389937/1‑150 (MQ=255)
                                                                                                     aaaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagcc                                           <  2:613146/150‑1 (MQ=255)
                                                                                                       aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGccagccag                                         >  1:541779/1‑150 (MQ=255)
                                                                                                         cAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACagccagcc                                       >  2:1915122/1‑150 (MQ=255)
                                                                                                          aaCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAgccagcca                                      >  1:776918/1‑150 (MQ=255)
                                                                                                           aCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCaa                                     >  1:1222503/1‑150 (MQ=255)
                                                                                                            ccccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAAc                                    >  1:1066591/1‑150 (MQ=255)
                                                                                                            ccccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAAc                                    >  2:245948/1‑150 (MQ=255)
                                                                                                             cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGCCAGCCAGCCAAcc                                   >  2:312156/1‑150 (MQ=255)
                                                                                                              ccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCg                                  <  2:1668199/150‑1 (MQ=255)
                                                                                                               cAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCga                                 >  1:1148182/1‑150 (MQ=255)
                                                                                                                   gATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCgaggag                             <  2:1148182/150‑1 (MQ=255)
                                                                                                                      ccGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGcc                          <  1:1744416/150‑1 (MQ=255)
                                                                                                                      ccGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACACTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGcc                          <  1:1513682/150‑1 (MQ=255)
                                                                                                                       cGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCa                         >  2:1005779/1‑150 (MQ=255)
                                                                                                                               gAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCcac                 >  1:1536248/1‑150 (MQ=255)
                                                                                                                                  ctggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACAcc              >  1:832528/7‑150 (MQ=255)
                                                                                                                                   aGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGGCGTCACCACTCCTTGCGTTGAAGGCAAAATACGGGTCGACATCCAGCCAACGGAGGAGGCCAGTTGCCACACCg             <  1:1143225/150‑1 (MQ=255)
                                                                                                                                    ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGtt                      >  2:938827/5‑140 (MQ=255)
                                                                                                                                    ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGtt                      <  1:938827/136‑1 (MQ=255)
                                                                                                                                    ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGc            >  2:1363072/5‑150 (MQ=255)
                                                                                                                                     ccAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGcc           <  2:417565/150‑1 (MQ=255)
                                                                                                                                     ccAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGcc           >  1:162454/1‑150 (MQ=255)
                                                                                                                                       aaTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCAc         <  1:643632/150‑1 (MQ=255)
                                                                                                                                       aaTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCAc         <  1:1389937/150‑1 (MQ=255)
                                                                                                                                       aaTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCAc         >  1:606528/1‑150 (MQ=255)
                                                                                                                                        aTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACg        <  2:606528/150‑1 (MQ=255)
                                                                                                                                         tCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGa       <  1:139101/150‑1 (MQ=255)
                                                                                                                                            gcgcCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGc    <  2:1924234/150‑1 (MQ=255)
                                                                                                                                            gcgcCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGc    <  2:1251885/150‑1 (MQ=255)
                                                                                                                                              gcccgccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCGGCTGGCGGAGAGCGTCCAGAACGCCGTCTCCACTCCTTCCGGTGAAGGCAAAATAGGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGccc  <  2:1608376/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCC  >  NC_000913/994485‑994776

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: