Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,792,490 T→G E81A (GAA→GCA)  nlpC ← putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,4900TG97.4% 666.8 / 5.7 195E81A (GAA→GCA) nlpCputative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (108/82);  total (110/85)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 196 total aligned reads displayed.

AAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACG  >  NC_000913/1792391‑1792586
                                                                                                   |                                                                                                
aaaCCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGttt                                                                                                  >  2:986190/1‑100 (MQ=255)
   ccATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCAc                                                                                               >  1:603359/1‑100 (MQ=255)
     aTTTTGTCCGGACCGCGTTTTGAAAAAGACCAGGTCAGCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACgg                                                                                             <  2:3123081/100‑1 (MQ=255)
     aTTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACgg                                                                                             >  1:5192670/1‑100 (MQ=255)
     aTTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACgg                                                                                             >  2:2140732/1‑100 (MQ=255)
      ttttGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACggg                                                                                            <  2:4702262/100‑1 (MQ=255)
         tGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCa                                                                                         <  1:2202765/100‑1 (MQ=255)
         tGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCa                                                                                         >  2:4427599/1‑100 (MQ=255)
          gTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAg                                                                                        >  1:1768278/1‑100 (MQ=255)
           tCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGc                                                                                       <  1:1492645/100‑1 (MQ=255)
            ccGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCt                                                                                      <  1:875096/100‑1 (MQ=255)
             cGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTg                                                                                     >  1:1665556/1‑100 (MQ=255)
             cGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTg                                                                                     >  1:5005327/1‑100 (MQ=255)
             cGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTg                                                                                     >  1:812843/1‑100 (MQ=255)
                 ccccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCaaa                                                                                 >  1:4778215/1‑100 (MQ=255)
                 ccccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCaaa                                                                                 >  2:890112/1‑100 (MQ=255)
                  cccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAAt                                                                                >  2:701845/1‑100 (MQ=255)
                  cccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAAt                                                                                <  2:4991810/100‑1 (MQ=255)
                  cccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAAt                                                                                <  2:596164/100‑1 (MQ=255)
                   ccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATc                                                                               <  1:4501260/100‑1 (MQ=255)
                   ccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATc                                                                               >  2:1289180/1‑100 (MQ=255)
                   ccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATc                                                                               <  2:262693/100‑1 (MQ=255)
                   ccGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATc                                                                               <  1:2068878/100‑1 (MQ=255)
                    cGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCg                                                                              >  1:2441350/1‑100 (MQ=255)
                     gTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTCTCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCGCGGGGCAGCTGCAAATCGa                                                                             >  2:197506/1‑100 (MQ=255)
                        ttGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCGGCTGCAAATCGAAAc                                                                          <  2:4173471/100‑1 (MQ=255)
                        ttGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAAc                                                                          >  2:4289285/1‑100 (MQ=255)
                         tGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACg                                                                         <  2:323462/100‑1 (MQ=255)
                         tGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACg                                                                         >  1:1020619/1‑100 (MQ=255)
                          gAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTACAAATCGAAACGa                                                                        >  2:1584659/1‑100 (MQ=255)
                           aaaaaGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGAt                                                                       >  2:2308025/1‑100 (MQ=255)
                               aGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATcgcg                                                                   <  1:999460/100‑1 (MQ=255)
                               aGACCAGGTCACCTGGCAGAGACTCGTCTTAATCAATTTGCGTGCCAATAGAGGCTGGTTGTTTGGTTGCACGGGACAGCTGCAAATCGAAACGATcgcg                                                                   <  2:2338724/100‑1 (MQ=255)
                                gACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCgcgc                                                                  >  1:4301345/1‑100 (MQ=255)
                                gACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCgcgc                                                                  <  2:2933505/100‑1 (MQ=255)
                                   cAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATc                                                               <  2:774448/100‑1 (MQ=255)
                                    aGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCg                                                              >  2:3162599/1‑100 (MQ=255)
                                    aGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCg                                                              <  2:1747545/100‑1 (MQ=255)
                                     ggTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGt                                                             >  1:5483557/1‑100 (MQ=255)
                                     ggTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGt                                                             >  2:4577920/1‑100 (MQ=255)
                                      gTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTc                                                            <  1:951025/100‑1 (MQ=255)
                                        cACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTcac                                                          <  2:3600710/100‑1 (MQ=255)
                                         aCCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCGTCGTcaca                                                         <  2:5274949/100‑1 (MQ=255)
                                          ccTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACaa                                                        <  2:1849705/100‑1 (MQ=255)
                                           cTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAAc                                                       <  1:4219306/100‑1 (MQ=255)
                                            tGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAAcc                                                      <  2:2845496/100‑1 (MQ=255)
                                            tGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAAcc                                                      >  1:4810665/1‑100 (MQ=255)
                                               cagcaACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACcaca                                                   <  1:4664348/100‑1 (MQ=255)
                                                 gcaACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACaaa                                                 >  2:3812631/1‑100 (MQ=255)
                                                    aCTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATcc                                              >  2:3671128/1‑100 (MQ=255)
                                                      tCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCg                                            >  1:2781518/1‑100 (MQ=255)
                                                        gTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGaa                                          >  1:4246492/1‑100 (MQ=255)
                                                        gTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGaa                                          >  2:1667214/1‑100 (MQ=255)
                                                         tCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAAc                                         <  1:2459677/100‑1 (MQ=255)
                                                           tttATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAg                                       <  1:3693109/100‑1 (MQ=255)
                                                           tttATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAg                                       >  2:3389982/1‑100 (MQ=255)
                                                           tttATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAg                                       >  2:3887579/1‑100 (MQ=255)
                                                            ttATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGt                                      >  2:4048315/1‑100 (MQ=255)
                                                             tATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTc                                     >  1:2196740/1‑100 (MQ=255)
                                                             tATCAATTTGCGTGCAAATAGAGGCTTGTTGCTTGGTTGCACGGGGCAGCCGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTc                                     <  1:538262/100‑1 (MQ=255)
                                                              aTCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTcc                                    <  2:5402735/100‑1 (MQ=255)
                                                                 aaTTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCcaca                                 >  2:2904276/1‑100 (MQ=255)
                                                                 aaTTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCcaca                                 >  1:5504728/1‑100 (MQ=255)
                                                                 aaTTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCcaca                                 >  1:4640671/1‑100 (MQ=255)
                                                                 aaTTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCcaca                                 <  1:1042106/100‑1 (MQ=255)
                                                                  aTTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCacac                                >  1:1955765/1‑100 (MQ=255)
                                                                       cGTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACAccgcgc                           >  1:4887538/1‑100 (MQ=255)
                                                                        gTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCNTCNTCACAANCACAAATCCCGNACAGTCCACAccgcgcc                          >  1:1198681/1‑100 (MQ=255)
                                                                        gTGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACAccgcgcc                          >  1:3251026/1‑100 (MQ=255)
                                                                         tGCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACAccgcgccg                         >  1:3392228/1‑100 (MQ=255)
                                                                          gCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACCAATCCCGAACAGTCCACAccgcgccgc                        >  2:4477195/1‑100 (MQ=255)
                                                                          gCCAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACAccgcgccgc                        >  1:4335519/1‑100 (MQ=255)
                                                                           ccAATAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACAccgcgccgcg                       <  2:5194853/100‑1 (MQ=255)
                                                                             aaTAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTc                     >  1:1244688/1‑100 (MQ=255)
                                                                               tAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCAt                   <  1:5487302/100‑1 (MQ=255)
                                                                               tAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCAt                   <  1:1306460/100‑1 (MQ=255)
                                                                               tAGAGGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCAt                   >  1:4120160/1‑100 (MQ=255)
                                                                                agagGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATg                  <  2:1511999/100‑1 (MQ=255)
                                                                                 gagGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGc                 <  1:4522878/100‑1 (MQ=255)
                                                                                 gagGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGc                 >  1:1196622/1‑100 (MQ=255)
                                                                                  agGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcc                >  2:1139656/1‑100 (MQ=255)
                                                                                  agGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcc                <  2:3707574/100‑1 (MQ=255)
                                                                                  agGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcc                <  2:3536498/100‑1 (MQ=255)
                                                                                  agGCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcc                <  1:1008821/100‑1 (MQ=255)
                                                                                   ggCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcca               >  2:13846/1‑100 (MQ=255)
                                                                                   ggCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGcca               >  1:3865604/1‑100 (MQ=255)
                                                                                    gCTTGTTGTTTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGccac              >  2:2093674/1‑100 (MQ=255)
                                                                                       tgttgtTTGGTTGCAAGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCAt           >  2:3771894/1‑100 (MQ=255)
                                                                                         ttgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATaa         >  1:4610514/1‑100 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        <  1:3561322/100‑1 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        >  2:5499498/1‑100 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        >  1:3964980/1‑100 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        <  1:1589487/100‑1 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        >  2:1382003/1‑100 (MQ=255)
                                                                                          tgttTGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAAc        >  2:2206651/1‑100 (MQ=255)
                                                                                           gttTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACg       >  1:2595988/1‑100 (MQ=255)
                                                                                            tttGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCACAACGa      >  1:3175830/1‑100 (MQ=255)
                                                                                             ttGGTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCTCGCCGCGTCATGCCACCATAACGAt     <  2:4382213/100‑1 (MQ=255)
                                                                                               ggTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATAc   >  2:4864334/1‑100 (MQ=255)
                                                                                                gTTGCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACg  >  2:111559/1‑100 (MQ=255)
                                                                                                   |                                                                                                
AAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAATCCCGAACAGTCCACACCGCGCCGCGTCATGCCACCATAACGATACG  >  NC_000913/1792391‑1792586

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