Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,507,076 G→C 5.5% L87L (CTC→CTG yhfT ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,507,0760GC 5.5% 669.8 / 25.3 182L87L (CTC→CTGyhfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (98/74);  new base C (5/5);  total (103/79)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 182 total aligned reads displayed.

GCGAGCTTAATTCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAG  >  NC_000913/3506931‑3507222
                                                                                                                                                 |                                                                                                                                                  
gCGAGCTTAATTCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGa                                                                                                                                                <  2:333532/150‑1 (MQ=255)
     cTTAATTCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCg                                                                                                                                           <  2:55804/150‑1 (MQ=255)
        aaTTCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTa                                                                                                                                        <  2:1135321/150‑1 (MQ=255)
           tCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGga                                                                                                                                     >  2:504434/1‑150 (MQ=255)
             ccccAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCCAGGATATCGGTAGgaag                                                                                                                                   >  1:1111455/1‑150 (MQ=255)
             ccccAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCCAGGATATCGGTAGgaag                                                                                                                                   >  1:341393/1‑150 (MQ=255)
             ccccAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAAAGCCATCAGGCTGTTTATCGCCAGGACGCCCAGGATATCGGTAGgaag                                                                                                                                   >  1:112157/1‑150 (MQ=255)
                  ggCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGaagaaaga                                                                                                                              >  1:1647563/1‑150 (MQ=255)
                  ggCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGaagaaaga                                                                                                                              >  1:1142647/1‑150 (MQ=255)
                   gCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAagag                                                                                                                             <  2:884200/150‑1 (MQ=255)
                   gCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAagag                                                                                                                             <  1:1448051/150‑1 (MQ=255)
                     tACCTAATACATCCACCGGCCGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTa                                                                                                                           >  1:1561827/1‑150 (MQ=255)
                      aCCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTaa                                                                                                                          <  2:289362/150‑1 (MQ=255)
                      aCCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTaa                                                                                                                          <  2:1316622/150‑1 (MQ=255)
                        cTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAAcc                                                                                                                        >  2:402687/1‑150 (MQ=255)
                        cTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCAACAGGCTGTTTATCGCCAGGACGCCGAGGAAATCGGTAGGAAGAAAGAGTAAcc                                                                                                                        >  2:402649/1‑150 (MQ=255)
                         tAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAATCGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCa                                                                                                                       >  2:275933/1‑150 (MQ=255)
                         tAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCa                                                                                                                       >  1:311895/1‑150 (MQ=255)
                         tAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCa                                                                                                                       >  2:522518/1‑150 (MQ=255)
                         tAATACATCCACCGGCAGCGCGCTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCa                                                                                                                       >  1:1067702/1‑150 (MQ=255)
                          aaTACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCAt                                                                                                                      >  1:2008866/1‑150 (MQ=255)
                           aTACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATg                                                                                                                     >  2:691183/1‑150 (MQ=255)
                           aTACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATg                                                                                                                     >  2:1965188/1‑150 (MQ=255)
                                 ccACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGa                                                                                                               >  1:556385/1‑150 (MQ=255)
                                 ccACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGa                                                                                                               <  1:2140684/150‑1 (MQ=255)
                                 ccACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGa                                                                                                               <  1:2166219/150‑1 (MQ=255)
                                 ccACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGa                                                                                                               <  1:471087/150‑1 (MQ=255)
                                  cACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGag                                                                                                              >  1:1899504/1‑150 (MQ=255)
                                  cACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGag                                                                                                              >  2:472337/1‑150 (MQ=255)
                                  cACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGag                                                                                                              >  1:396403/1‑150 (MQ=255)
                                  cACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGag                                                                                                              >  1:1549554/1‑150 (MQ=255)
                                    ccGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGagca                                                                                                            <  2:570599/150‑1 (MQ=255)
                                           cgcgGTCAGCAGCTGGTTTACTGGCAACAGTCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGt                                                                                                     <  1:328535/150‑1 (MQ=255)
                                           cgcgGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGt                                                                                                     <  1:694400/150‑1 (MQ=255)
                                           cgcgGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGt                                                                                                     <  1:1631348/150‑1 (MQ=255)
                                            gcgGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGtt                                                                                                    <  2:820885/150‑1 (MQ=255)
                                            gcgGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGAACAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGtt                                                                                                    <  1:647722/150‑1 (MQ=255)
                                                 cagcagCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAg                                                                                               >  1:1371598/1‑150 (MQ=255)
                                                  agcagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAgg                                                                                              >  1:1464070/1‑150 (MQ=255)
                                                  agcagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAgg                                                                                              <  1:2373996/150‑1 (MQ=255)
                                                  agcagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAgg                                                                                              >  1:1018964/1‑150 (MQ=255)
                                                   gcagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAggg                                                                                             >  1:1618913/1‑150 (MQ=255)
                                                    cagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGt                                                                                            >  1:1182651/1‑150 (MQ=255)
                                                    cagcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGt                                                                                            >  1:1182683/1‑150 (MQ=255)
                                                     agcTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTa                                                                                           >  2:586177/1‑150 (MQ=255)
                                                            ttACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAgaga                                                                                    <  2:1810694/150‑1 (MQ=255)
                                                                tGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATccc                                                                                >  2:691138/1‑150 (MQ=255)
                                                                tGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATccc                                                                                >  2:1415032/1‑150 (MQ=255)
                                                                tGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATccc                                                                                >  2:573971/1‑150 (MQ=255)
                                                                tGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATccc                                                                                >  2:1939807/1‑150 (MQ=255)
                                                                tGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATccc                                                                                >  2:1503881/1‑150 (MQ=255)
                                                                       aGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAg                                                                         >  1:1837319/1‑150 (MQ=255)
                                                                         gCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGc                                                                       >  1:52467/1‑150 (MQ=255)
                                                                            aGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCcac                                                                    >  1:1807665/1‑150 (MQ=255)
                                                                                   aTCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCCCAAAACCa                                                             >  2:1263422/1‑150 (MQ=255)
                                                                                   aTCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCa                                                             >  2:902694/1‑150 (MQ=255)
                                                                                   aTCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCa                                                             <  1:1629923/150‑1 (MQ=255)
                                                                                     cAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAAt                                                           <  1:731384/150‑1 (MQ=255)
                                                                                        cacGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCt                                                        >  1:333176/1‑150 (MQ=255)
                                                                                          cGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCa                                                      <  1:1424434/150‑1 (MQ=255)
                                                                                          cGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCa                                                      <  1:1929349/150‑1 (MQ=255)
                                                                                                 gATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCaaa                                               >  1:749218/1‑150 (MQ=255)
                                                                                                   tAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAg                                             <  2:741951/150‑1 (MQ=255)
                                                                                                   tAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAg                                             <  2:1561827/150‑1 (MQ=255)
                                                                                                    aGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGc                                            <  2:2118242/150‑1 (MQ=255)
                                                                                                    aGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGc                                            >  2:2043222/1‑150 (MQ=255)
                                                                                                     gcgcCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCg                                           >  2:1250919/1‑150 (MQ=255)
                                                                                                     gcgcCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCg                                           >  2:2289179/1‑150 (MQ=255)
                                                                                                      cgcCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGa                                          <  1:2378982/150‑1 (MQ=255)
                                                                                                             aCCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGGTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTccc                                   >  2:2100716/1‑150 (MQ=255)
                                                                                                               cAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCg                                 >  1:1815936/1‑150 (MQ=255)
                                                                                                               cAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCg                                 >  2:1360213/1‑150 (MQ=255)
                                                                                                               cAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCg                                 >  2:1354025/1‑150 (MQ=255)
                                                                                                                aaaCGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGc                                >  1:1955550/1‑150 (MQ=255)
                                                                                                                  aCGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCt                              <  1:2246684/150‑1 (MQ=255)
                                                                                                                   cGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTc                             <  1:586177/150‑1 (MQ=255)
                                                                                                                    gCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTcg                            <  2:556385/150‑1 (MQ=255)
                                                                                                                           ggCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGt                     >  2:1667302/1‑150 (MQ=255)
                                                                                                                           ggCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGt                     >  2:2017027/1‑150 (MQ=255)
                                                                                                                            gCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGtt                    >  1:900138/1‑150 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  2:1602852/150‑1 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  2:222592/150‑1 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  1:1711866/150‑1 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   >  1:2048076/1‑150 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCCAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  2:1111455/150‑1 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCCAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  2:341393/150‑1 (MQ=255)
                                                                                                                             cTGTTTATCGCCAGGACGCCCAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTc                   <  2:1111467/150‑1 (MQ=255)
                                                                                                                              tGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCa                  <  1:270942/150‑1 (MQ=255)
                                                                                                                              tGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCa                  <  1:867284/150‑1 (MQ=255)
                                                                                                                                  tATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCATCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCAtata              <  1:275933/150‑1 (MQ=255)
                                                                                                                                      gCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCt          >  1:2067295/1‑150 (MQ=255)
                                                                                                                                      gCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCt          >  2:2273935/1‑150 (MQ=255)
                                                                                                                                      gCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCt          >  1:904119/1‑150 (MQ=255)
                                                                                                                                      gCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCt          >  1:904144/1‑150 (MQ=255)
                                                                                                                                      gCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCt          >  1:1508193/1‑150 (MQ=255)
                                                                                                                                       ccAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCtt         >  2:619290/1‑150 (MQ=255)
                                                                                                                                        cAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTc        >  1:57418/1‑150 (MQ=255)
                                                                                                                                          ggACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCaa      <  1:1390198/150‑1 (MQ=255)
                                                                                                                                             cGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCa   <  1:2029944/150‑1 (MQ=255)
                                                                                                                                              gCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAg  >  2:2354915/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCGAGCTTAATTCCCCCAGGCTACCTAATACATCCACCGGCAGCGCGGTCAGCAGCTGGTTTACTGGCAACAGGCAAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAG  >  NC_000913/3506931‑3507222

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