Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,400,849 A→T 5.9% V80D (GTT→GAT)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,8490AT 5.9% 729.0 / 25.2 202V80D (GTT→GAT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (98/91);  new base T (7/5);  total (106/96)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 205 total aligned reads displayed.

GGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGAT  >  NC_000913/3400701‑3400997
                                                                                                                                                    |                                                                                                                                                    
ggATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATaa                                                                                                                                                      <  1:393140/149‑1 (MQ=255)
    tCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                                                                                                                  <  3:347499/149‑1 (MQ=255)
      aCGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGt                                                                                                                                                >  4:1121606/1‑149 (MQ=255)
        gAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                                                              <  1:845019/149‑1 (MQ=255)
        gAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                                                              >  3:918819/1‑149 (MQ=255)
        gAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                                                              <  1:317319/149‑1 (MQ=255)
               gcggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCttt                                                                                                                                            >  4:1604397/1‑144 (MQ=255)
                cggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                                                              >  4:1412742/1‑141 (MQ=255)
                cggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                                                              <  3:1412742/141‑1 (MQ=255)
                     ttCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGaa                                                                                                                                 <  4:1746826/149‑1 (MQ=255)
                      tCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                                                                                                                                <  2:1485825/149‑1 (MQ=255)
                          ccTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                                                                                                                            <  2:1776265/149‑1 (MQ=255)
                             gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcagca                                                                                                                         >  3:827960/1‑149 (MQ=255)
                                    gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGc                                                                                                                  >  4:1004575/1‑149 (MQ=255)
                                      cAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTgg                                                                                                                                      <  2:1881296/127‑1 (MQ=255)
                                      cAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGccc                                                                                                                >  2:424709/1‑149 (MQ=255)
                                                acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTa                                                                                                          <  1:1103407/145‑1 (MQ=255)
                                                 caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATAtt                                                                                                                    >  1:45994/1‑134 (MQ=255)
                                                 caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATAtt                                                                                                                    <  2:45994/134‑1 (MQ=255)
                                                 caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCg                                                                                                               <  4:1454807/139‑1 (MQ=255)
                                                 caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGcc                                                                                                     <  4:674354/149‑1 (MQ=255)
                                                   aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTgg                                                                                                                                      >  2:335170/1‑114 (MQ=255)
                                                    aCCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGc                                                                                                                  <  4:775418/133‑1 (MQ=255)
                                                    aCCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAg                                                                                                  >  3:563586/1‑149 (MQ=255)
                                                      cAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCa                                                                                                >  4:1354447/1‑149 (MQ=255)
                                                      cAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCa                                                                                                <  1:622250/149‑1 (MQ=255)
                                                         aaCGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCt                                                                                             >  1:1695043/1‑149 (MQ=255)
                                                              gcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAata                                                                                                                      <  2:614897/119‑1 (MQ=255)
                                                                 ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTcgcg                                                                                     >  3:441140/1‑149 (MQ=255)
                                                                    ttGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCa                                                                                  <  2:1695043/149‑1 (MQ=255)
                                                                    ttGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTAGGCGAAt                                                                                                                                <  2:711032/103‑1 (MQ=255)
                                                                     tGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGTAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCAt                                                                                 <  1:996034/149‑1 (MQ=255)
                                                                       gACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGa                                                                               >  1:546785/1‑149 (MQ=255)
                                                                        aCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGAc                                                                              <  2:206993/149‑1 (MQ=255)
                                                                        aCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGAGTCATGAc                                                                              <  4:918819/149‑1 (MQ=255)
                                                                            aTAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGc                                                                                                                                                   >  3:6437/1‑76 (MQ=255)
                                                                            aTAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGc                                                                                                                                                   <  4:6437/76‑1 (MQ=255)
                                                                             tAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGGAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                                                                         >  2:1322598/1‑149 (MQ=255)
                                                                             tAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcagc                                                                                                                          <  1:1553908/100‑1 (MQ=255)
                                                                             tAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                                                                         <  4:563586/149‑1 (MQ=255)
                                                                              aaaGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                                                                   <  4:226503/138‑1 (MQ=255)
                                                                              aaaGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                                                                   >  3:226503/1‑138 (MQ=255)
                                                                               aaGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGc                                                                       >  3:313497/1‑149 (MQ=255)
                                                                                aGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTc                                                                                                                                           >  3:736940/1‑80 (MQ=255)
                                                                                   tGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                                                                   >  3:682385/1‑133 (MQ=255)
                                                                                    gTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTa                                                                  >  2:942606/1‑149 (MQ=255)
                                                                                     ttGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGc                                                                                           <  2:1052924/123‑1 (MQ=255)
                                                                                       gCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGc                                                               >  1:103806/1‑149 (MQ=255)
                                                                                       gCTGTGCACTTGTTTGATAAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                                                                    >  4:1786995/1‑112 (MQ=255)
                                                                                       gCTGTGCACTTGTTTGATAAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                                                                    <  3:1786995/112‑1 (MQ=255)
                                                                                         tgtgCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCtt                                                             >  3:1486934/1‑149 (MQ=255)
                                                                                          gtgCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGGAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCttt                                                            <  1:1322598/149‑1 (MQ=255)
                                                                                           tgCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                                                                    <  1:1186142/108‑1 (MQ=255)
                                                                                           tgCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTAc                                                                 >  1:213367/1‑143 (MQ=255)
                                                                                             cACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCg                                                         <  1:559700/149‑1 (MQ=255)
                                                                                              aCTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                                                                   >  4:1254093/1‑122 (MQ=255)
                                                                                              aCTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                                                                   >  4:343075/1‑122 (MQ=255)
                                                                                              aCTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTTgg                                                        <  2:1163346/149‑1 (MQ=255)
                                                                                               cTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt                                                       >  2:691841/1‑149 (MQ=255)
                                                                                               cTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt                                                       >  1:1840300/1‑149 (MQ=255)
                                                                                                ttGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCAt                                                                                 <  1:702302/122‑1 (MQ=255)
                                                                                                ttGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTg                                                      >  1:1291680/1‑149 (MQ=255)
                                                                                                   tttGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAAc                                                   <  4:441140/149‑1 (MQ=255)
                                                                                                     tgatgaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCg                                                 >  3:1346515/1‑149 (MQ=255)
                                                                                                       atgaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                                                                                                                             <  3:1560282/71‑1 (MQ=255)
                                                                                                       atgaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                                                                                                                             >  4:1560282/1‑71 (MQ=255)
                                                                                                        tgaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCa                                              >  2:1447620/1‑149 (MQ=255)
                                                                                                         gaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGa                                                     <  1:1432253/141‑1 (MQ=255)
                                                                                                         gaAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGa                                                     >  2:1432253/1‑141 (MQ=255)
                                                                                                           aGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt                                                       >  3:468676/1‑137 (MQ=255)
                                                                                                           aGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt                                                       <  4:468676/137‑1 (MQ=255)
                                                                                                                      aTTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTg                                                                                            >  3:1657085/1‑89 (MQ=255)
                                                                                                                      aTTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTg                                                                                            <  4:1657085/89‑1 (MQ=255)
                                                                                                                       ttttttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATc                                         <  3:1438029/139‑1 (MQ=255)
                                                                                                                           ttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGcca                           <  4:951368/149‑1 (MQ=255)
                                                                                                                              aGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGaa                                 >  1:72848/1‑140 (MQ=255)
                                                                                                                                tCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACg                      >  1:1584334/1‑149 (MQ=255)
                                                                                                                                acACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCTGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTc                                                                                        >  2:1592426/2‑83 (MQ=255)
                                                                                                                                   cGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGaa                   <  3:1213030/149‑1 (MQ=255)
                                                                                                                                    gaagaaGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGa                                     <  2:814664/130‑1 (MQ=255)
                                                                                                                                    gaagaaGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAg                  >  4:1198228/1‑149 (MQ=255)
                                                                                                                                      gaaagACGGCGTTATCGCCGACTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGaa                   <  2:1572012/144‑1 (MQ=255)
                                                                                                                                      agaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCAc                                             <  1:491370/120‑1 (MQ=255)
                                                                                                                                      agaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCAc                                             >  2:491370/1‑120 (MQ=255)
                                                                                                                                      agaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGc                >  1:201270/1‑149 (MQ=255)
                                                                                                                                       gaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCt               >  2:999754/1‑149 (MQ=255)
                                                                                                                                         agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCa             <  3:1354447/149‑1 (MQ=255)
                                                                                                                                         agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCa             >  2:1291584/1‑149 (MQ=255)
                                                                                                                                         agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCa             <  3:1103994/149‑1 (MQ=255)
                                                                                                                                             ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAg                                                                     >  3:1504760/1‑89 (MQ=255)
                                                                                                                                             ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAg                                                                     <  4:1504760/89‑1 (MQ=255)
                                                                                                                                             ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCa             <  4:695409/145‑1 (MQ=255)
                                                                                                                                             ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCa         >  3:1765544/1‑149 (MQ=255)
                                                                                                                                              gCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCaa        >  4:297679/1‑149 (MQ=255)
                                                                                                                                               cGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAAt       <  2:1840300/149‑1 (MQ=255)
                                                                                                                                                 aTAGCGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt                                                       >  4:316560/1‑99 (MQ=255)
                                                                                                                                                 aTAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                                                                    >  2:1269981/1‑54 (MQ=255)
                                                                                                                                                 aTAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                                                                    <  1:1269981/54‑1 (MQ=255)
                                                                                                                                                 aTAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATAc     >  1:1794468/1‑149 (MQ=255)
                                                                                                                                                    aCGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGAt  <  2:546785/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGAT  >  NC_000913/3400701‑3400997

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