Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 105,046 A→C 5.4% Q355H (CAA→CAC ftsA → ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913105,0460AC 5.4% 376.7 / 8.8 123Q355H (CAA→CACftsAATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base A (61/55);  new base C (5/2);  total (66/57)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 123 total aligned reads displayed.

TGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCA  >  NC_000913/104901‑105193
                                                                                                                                                 |                                                                                                                                                   
tGGTCAACGCAGAGATATTGCAGTTGCATGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTATCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGc                                                                                                                                                 <  1:115790/150‑1 (MQ=255)
tGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGc                                                                                                                                                 <  1:1472485/150‑1 (MQ=255)
  gTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGAAGATCGAAGGTCTTGCAGCCTGTGCACAGCGCGTGTTTCATACGCAAGTGCGt                                                                                                                                               >  2:913302/1‑150 (MQ=255)
    cAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTAt                                                                                                                                             >  1:1167597/1‑150 (MQ=255)
    cAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTAt                                                                                                                                             <  2:205057/150‑1 (MQ=255)
     aaCGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCTCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATc                                                                                                                                            >  2:152782/1‑150 (MQ=255)
     aaCGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATc                                                                                                                                            >  2:1561270/1‑150 (MQ=255)
     aaCGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATc                                                                                                                                            <  2:1691696/150‑1 (MQ=255)
       cGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCgg                                                                                                                                          >  2:1621438/1‑150 (MQ=255)
         aaGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACTAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGgcg                                                                                                                                        >  1:117587/1‑150 (MQ=255)
         aaGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCACTTGCGTAtgaaac                                                                                                                                        <  2:656500/150‑6 (MQ=255)
           gagaTATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGgcgcg                                                                                                                                      >  2:1671931/1‑150 (MQ=255)
            agaTATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGgcgcgc                                                                                                                                     >  2:1432271/1‑150 (MQ=255)
               tatTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGcgccgc                                                                                                                                  >  1:1427947/1‑150 (MQ=255)
                atTGCAGTTGCAGGAATAGCTTAGCCAAAAAGGGGTTAAACATAACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGAGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGAAAGTGCGTATCGGCGCGCCGCt                                                                                                                                 <  2:1053268/150‑1 (MQ=255)
                 ttgcagttgcagGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTg                                                                                                                                >  1:1514069/1‑150 (MQ=255)
                        tgcagGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTa                                                                                                                         >  2:955703/1‑150 (MQ=255)
                          cagGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCACAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTAcc                                                                                                                       <  2:29227/150‑1 (MQ=255)
                           agGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCg                                                                                                                      >  1:1660089/1‑150 (MQ=255)
                            ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCAGGCGTCAGGCATTGTATAAACCGGTGGCGCAGCGAAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCCTATCGGAGCGCCGCTGAACATTACCgg                                                                                                                     >  1:1021233/1‑150 (MQ=255)
                             gAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGt                                                                                                                    <  1:837515/150‑1 (MQ=255)
                                   cTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTAGCGGTTTAACg                                                                                                              >  1:259754/1‑150 (MQ=255)
                                   cTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACg                                                                                                              <  1:955703/150‑1 (MQ=255)
                                   cTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACg                                                                                                              <  1:826358/150‑1 (MQ=255)
                                   cTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCGTGCAGCCTGTGCTCAGCGCGTGTTTCAAACGCACGTGCGTATCGGCGCGCCGCTGAACATTCCCGGTTTACCg                                                                                                              >  1:263585/1‑150 (MQ=255)
                                     tCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGTCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGa                                                                                                            <  1:1545270/150‑1 (MQ=255)
                                     tCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGa                                                                                                            <  2:576249/150‑1 (MQ=255)
                                      cGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGAt                                                                                                           <  1:824445/150‑1 (MQ=255)
                                         caacaaGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTAt                                                                                                        <  1:497968/150‑1 (MQ=255)
                                              aGGGGTTAAACATCACCTGGCGGCAGGATTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGAATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCa                                                                                                   <  1:1538526/150‑1 (MQ=255)
                                                      aaCATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGAGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGt                                                                                           >  2:1086146/1‑150 (MQ=255)
                                                        cttCAGCTGGCGGCAGGAATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCTGCCTGTCCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGtat                                                                                         <  2:1626002/148‑1 (MQ=255)
                                                         aTCCCCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGTAGGTCTTGCAGCCTGTGCTGAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGtatt                                                                                        <  1:913302/150‑1 (MQ=255)
                                                               tGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGa                                                                                  <  2:392362/150‑1 (MQ=255)
                                                                 gcggcAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGAcg                                                                                >  2:906596/1‑150 (MQ=255)
                                                                   ggcaggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGAcggc                                                                              <  1:1635466/150‑1 (MQ=255)
                                                                     caggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGAcggcgg                                                                            <  1:1621438/150‑1 (MQ=255)
                                                                       ggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTg                                                                          <  1:1671931/150‑1 (MQ=255)
                                                                       ggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTg                                                                          >  1:478164/1‑150 (MQ=255)
                                                                       ggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGGCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTg                                                                          <  1:314188/150‑1 (MQ=255)
                                                                         caTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTATTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTggg                                                                        >  2:557284/1‑150 (MQ=255)
                                                                            tGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTTACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGAtt                                                                     >  2:591286/1‑150 (MQ=255)
                                                                            tGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGAtt                                                                     >  2:188453/1‑150 (MQ=255)
                                                                            tGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGAtt                                                                     >  2:681558/1‑150 (MQ=255)
                                                                            tGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGAtt                                                                     >  1:367286/1‑150 (MQ=255)
                                                                                 tAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTc                                                                <  1:1587014/150‑1 (MQ=255)
                                                                                 tAACCGGTGGCGCAGCGCAGATCGAAGGGCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTc                                                                <  1:105976/150‑1 (MQ=255)
                                                                                  aaCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCa                                                               >  1:1389968/1‑150 (MQ=255)
                                                                                  aaCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCa                                                               <  1:1068970/150‑1 (MQ=255)
                                                                                   aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCAc                                                              <  1:1000656/150‑1 (MQ=255)
                                                                                    ccGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACt                                                             <  1:3329/150‑1 (MQ=255)
                                                                                    ccGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACt                                                             >  2:372524/1‑150 (MQ=255)
                                                                                     cGGTGGCGCAGCGCAGATCGAAGGTCTTGTAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTa                                                            >  1:65703/1‑150 (MQ=255)
                                                                                     cGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTa                                                            <  2:1427947/150‑1 (MQ=255)
                                                                                       gTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCGCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATg                                                          <  2:1422209/150‑1 (MQ=255)
                                                                                        tGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATgg                                                         >  1:1304406/1‑150 (MQ=255)
                                                                                           cgcagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGaa                                                      >  1:713014/1‑150 (MQ=255)
                                                                                             cagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAag                                                    >  1:725832/1‑150 (MQ=255)
                                                                                             cagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAag                                                    >  1:1346518/1‑150 (MQ=255)
                                                                                               gcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAagag                                                  >  1:492082/1‑150 (MQ=255)
                                                                                                cgcagAGCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGt                                                 <  1:160351/150‑1 (MQ=255)
                                                                                                   agATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTcac                                              <  2:593547/150‑1 (MQ=255)
                                                                                                    gATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTcaca                                             >  2:1216106/1‑150 (MQ=255)
                                                                                                       cGAAGGTCTTGCAGCCTGTGTTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCt                                          >  1:797018/1‑150 (MQ=255)
                                                                                                       cGAAGGTCTTGCAGCCTGTGCTCAGCGGGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACTGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCt                                          <  1:706120/150‑1 (MQ=255)
                                                                                                        gAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCATCACTATGGGAAAGAGTCACATCtt                                         <  1:1432271/150‑1 (MQ=255)
                                                                                                        gAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGGGGGATTGCTTCCCTATGGGAAAGAGTCACATCtt                                         >  2:542217/1‑150 (MQ=255)
                                                                                                              cTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGt                                   <  2:855834/150‑1 (MQ=255)
                                                                                                                tGCAGCCTGTGCTCCGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGa                                 >  1:1239598/1‑150 (MQ=255)
                                                                                                                tGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGa                                 >  2:407528/1‑150 (MQ=255)
                                                                                                                 gCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGaa                                >  2:136030/1‑150 (MQ=255)
                                                                                                                  cAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAg                               <  1:85505/150‑1 (MQ=255)
                                                                                                                    gCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGCACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCCCATCTTAACGGTGAAGCt                             >  1:883374/1‑150 (MQ=255)
                                                                                                                      cTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGa                           <  1:803125/150‑1 (MQ=255)
                                                                                                                        gtgCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAg                         >  1:1704850/1‑150 (MQ=255)
                                                                                                                            tCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGa                     <  2:117587/150‑1 (MQ=255)
                                                                                                                             cAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaa                    >  2:801543/1‑150 (MQ=255)
                                                                                                                             cAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaa                    >  1:985801/1‑150 (MQ=255)
                                                                                                                              aGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaaa                   <  2:348549/150‑1 (MQ=255)
                                                                                                                                cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaaaaa                 >  1:337607/1‑150 (MQ=255)
                                                                                                                                  cgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACg               <  1:681558/150‑1 (MQ=255)
                                                                                                                                  cgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACg               <  2:375776/150‑1 (MQ=255)
                                                                                                                                   cgcgccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgt              >  1:1003769/8‑150 (MQ=255)
                                                                                                                                    tgtTTGATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTTACGGATTATGCTGATGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTATAAAAACgtg             <  2:1021233/150‑1 (MQ=255)
                                                                                                                                    tgtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGTGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgtg             >  2:133778/1‑150 (MQ=255)
                                                                                                                                     gtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgtgt            <  1:136075/150‑1 (MQ=255)
                                                                                                                                      tttCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTAGCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGGAAGAGGCACATCTTAACGGTGATGCTGAAGTAGAAAACGGTGtt           >  1:311544/1‑150 (MQ=255)
                                                                                                                                      tttCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGtt           >  2:672220/1‑150 (MQ=255)
                                                                                                                                      tttCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGtt           >  1:595793/1‑150 (MQ=255)
                                                                                                                                      gccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGtt           >  1:479551/5‑150 (MQ=255)
                                                                                                                                       ccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTa          >  1:1259992/4‑150 (MQ=255)
                                                                                                                                       ccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCAAATCTTAACGGTGAAGCTGAAGTAGAAAAACGTTTTa          >  2:251681/4‑150 (MQ=255)
                                                                                                                                           tACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGc      >  1:1340025/1‑150 (MQ=255)
                                                                                                                                           tACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGc      <  2:1389968/150‑1 (MQ=255)
                                                                                                                                           tACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGc      >  1:728114/1‑150 (MQ=255)
                                                                                                                                            aCGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCCGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCa     <  2:215101/150‑1 (MQ=255)
                                                                                                                                             cGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGa                           >  1:117356/1‑127 (MQ=255)
                                                                                                                                             cGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGa                           <  2:117356/127‑1 (MQ=255)
                                                                                                                                              gCAAGGGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATc   <  1:618879/150‑1 (MQ=255)
                                                                                                                                               cAAGTGCGTAGCGGCGCGCCGCGGAACATTACCGGTTTAACGGATTATGCGCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCa  <  2:286003/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCA  >  NC_000913/104901‑105193

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