Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,159,481 A→G M402T (ATG→ACG)  sufI ← suppressor of ftsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,159,4810AG99.4% 1007.8 / ‑3.6 328M402T (ATG→ACG) sufIsuppressor of ftsI
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (143/183);  total (144/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 342 total aligned reads displayed.

AGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATACGG  >  NC_000913/3159391‑3159570
                                                                                          |                                                                                          
agcTCCACTTGTCCGTCACCCCAAACGGTAGCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGtt                                                                                          <  1:7904914/93‑1 (MQ=255)
agcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGtt                                                                                          >  1:130181/1‑93 (MQ=255)
 gcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTa                                                                                         >  2:4892863/1‑93 (MQ=255)
 gcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTa                                                                                         <  1:7755072/93‑1 (MQ=255)
 gcTCCACTTGTCCGTCAACCCAAACGGTATCTTCCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTa                                                                                         <  1:2915144/93‑1 (MQ=255)
  cTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTAc                                                                                        <  1:4492814/93‑1 (MQ=255)
  cTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTAc                                                                                        >  1:2995231/1‑93 (MQ=255)
   tCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACg                                                                                       >  1:4000149/1‑93 (MQ=255)
   tCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACg                                                                                       >  2:6307436/1‑93 (MQ=255)
   tCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACg                                                                                       >  2:2458241/1‑93 (MQ=255)
    ccACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGc                                                                                      >  1:6760017/1‑93 (MQ=255)
     cACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGcc                                                                                     >  1:6538023/1‑93 (MQ=255)
      aCTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCt                                                                                    >  2:7623451/1‑93 (MQ=255)
      aCTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCt                                                                                    >  2:1242196/1‑93 (MQ=255)
      aCTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCt                                                                                    >  2:4584087/1‑93 (MQ=255)
       cTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCtt                                                                                   >  1:885847/1‑93 (MQ=255)
       cTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCtt                                                                                   >  2:5671983/1‑93 (MQ=255)
        ttGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTc                                                                                  >  1:7360469/1‑93 (MQ=255)
        ttGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTc                                                                                  >  1:6924192/1‑93 (MQ=255)
        ttGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTc                                                                                  >  2:3760005/1‑93 (MQ=255)
        ttGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTc                                                                                  >  1:1097365/1‑93 (MQ=255)
         tGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACAGATCTGGAACGTTACGCCTTCa                                                                                 >  2:6451145/1‑93 (MQ=255)
           tCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAat                                                                               >  2:1794283/1‑93 (MQ=255)
           tCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAat                                                                               >  2:735943/1‑93 (MQ=255)
            ccGTCACCCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAata                                                                              <  1:4226861/93‑1 (MQ=255)
               tCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATgg                                                                           >  2:5558998/1‑93 (MQ=255)
                cAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGa                                                                          <  2:7500065/93‑1 (MQ=255)
                cAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGa                                                                          <  1:1862232/93‑1 (MQ=255)
                 accccAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGaa                                                                         <  2:1339484/91‑1 (MQ=255)
                 aaCCCAAACGGTATTTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGaa                                                                         <  1:2611696/93‑1 (MQ=255)
                 aaCCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGaa                                                                         <  2:4343077/93‑1 (MQ=255)
                  aCCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAAc                                                                        <  1:5162066/93‑1 (MQ=255)
                    ccAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGc                                                                      <  1:534706/93‑1 (MQ=255)
                    ccAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGc                                                                      <  1:1538686/93‑1 (MQ=255)
                     cAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCt                                                                     >  1:6002703/1‑93 (MQ=255)
                     cAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCt                                                                     >  1:4982904/1‑93 (MQ=255)
                     cAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCt                                                                     <  2:7459701/93‑1 (MQ=255)
                     cAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCt                                                                     >  2:6573975/1‑93 (MQ=255)
                      aaaCGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCtt                                                                    >  1:1815416/1‑93 (MQ=255)
                       agcGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTg                                                                   <  1:7662209/91‑1 (MQ=255)
                       aaCGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTg                                                                   <  2:7174660/93‑1 (MQ=255)
                        aCGGTATCTTCCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGc                                                                  <  1:484374/93‑1 (MQ=255)
                         cGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCg                                                                 <  2:41420/93‑1 (MQ=255)
                         cGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCg                                                                 >  1:5734616/1‑93 (MQ=255)
                          ggTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCgg                                                                <  2:10160/93‑1 (MQ=255)
                          ggTATCTTTCCAGCCTCTGTCTTCCGGAAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCgg                                                                <  2:3415268/93‑1 (MQ=255)
                           gTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGc                                                               >  2:7114325/1‑93 (MQ=255)
                            tATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCt                                                              >  1:5314810/1‑93 (MQ=255)
                              tCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTcc                                                            >  2:6432586/1‑92 (MQ=255)
                               ttttCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGt                                                           <  1:1132969/92‑1 (MQ=255)
                               cTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTccc                                                           >  2:4005706/1‑91 (MQ=255)
                                tttCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGctcctc                                                          >  2:7715984/1‑93 (MQ=255)
                                tttCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTc                                                          <  1:5639020/93‑1 (MQ=255)
                                 ttCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGttc                                                         >  1:5863552/1‑91 (MQ=255)
                                 ttCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTcc                                                         <  1:7298581/93‑1 (MQ=255)
                                     aGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCcgcgc                                                     <  2:6091175/93‑1 (MQ=255)
                                     aGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCcgcgc                                                     <  2:2811117/93‑1 (MQ=255)
                                      gCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCcccgcg                                                    >  1:3017369/1‑93 (MQ=255)
                                         tctGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCcgcccgcgc                                                 >  2:2606050/1‑90 (MQ=255)
                                           tGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCCCCCCGCGGACCg                                               >  1:4306794/1‑93 (MQ=255)
                                               ttCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCa                                           <  2:237924/93‑1 (MQ=255)
                                                tCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAg                                          <  2:1937652/93‑1 (MQ=255)
                                                 ccGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGCCCGCCCGTACCGTcccac                                         >  2:3040444/1‑90 (MQ=255)
                                                  cGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCg                                        <  1:2200756/93‑1 (MQ=255)
                                                      aaCGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCCCCCCGACCCTCCAGCggttc                                    >  1:5473770/1‑89 (MQ=255)
                                                       aCGGCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTccc                                  <  2:7254616/93‑1 (MQ=255)
                                                        cGGCATCGCTCCATTCACGTTAGGGAACTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCa                                 <  2:2349734/93‑1 (MQ=255)
                                                         ggCATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCAc                                <  1:214582/93‑1 (MQ=255)
                                                          gCATCGCGCCATTCACGTTATGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCGCCCCGCCGACCGTCCACCGCTTccccc                                >  2:1885255/1‑91 (MQ=255)
                                                           cATCGCGCCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGt                              <  2:362480/93‑1 (MQ=255)
                                                               gcgcCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCt                          <  1:362965/93‑1 (MQ=255)
                                                               gcgcCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCt                          <  1:3279949/93‑1 (MQ=255)
                                                               gcgcCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCt                          <  1:7512761/93‑1 (MQ=255)
                                                               ccgcCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCt                          <  2:6758782/92‑1 (MQ=255)
                                                                cgcCATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCtt                         <  2:783057/93‑1 (MQ=255)
                                                                  ccATTCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCGCGCGCAGCGTCCCGCGGTTCCACCTTcccttc                       >  2:5293896/1‑89 (MQ=255)
                                                                     ttCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTtgctgc                    <  2:4996972/93‑1 (MQ=255)
                                                                     ttCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCGTACCGTCCAGCGGTTCCCCCCTCCTttctgc                    >  1:2815441/1‑93 (MQ=255)
                                                                     ttCACGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCGCGCGGACCGTCCACCGGTCCCCCGTTCCTtgctgc                    >  1:1722248/1‑93 (MQ=255)
                                                                     ttCACGTTACGGAGCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGGCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTtgctgc                    <  1:6382462/93‑1 (MQ=255)
                                                                       cACGTTACGGATCTGGACCGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTgcgc                  <  2:4279892/93‑1 (MQ=255)
                                                                        aCGTTACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCGCGCGTACCGTCCAGCGTTCCCACCCTCCTTGCTGcgcg                 >  2:4733570/1‑93 (MQ=255)
                                                                           ttACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTg              <  2:4815781/93‑1 (MQ=255)
                                                                            tACGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGa             <  2:3920334/93‑1 (MQ=255)
                                                                              cGGATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACa           <  1:6241708/93‑1 (MQ=255)
                                                                               ggATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACAt          <  1:1028390/93‑1 (MQ=255)
                                                                                gATCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATc         <  2:2859597/93‑1 (MQ=255)
                                                                                  tCTGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCaa       <  1:3169775/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  2:1568497/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  2:5858347/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  2:3734980/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  1:4878739/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  1:337325/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  2:5391/93‑1 (MQ=255)
                                                                                    tGGAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATa     <  2:1119869/93‑1 (MQ=255)
                                                                                     ggAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATAc    <  2:7056666/93‑1 (MQ=255)
                                                                                     ggAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATAc    <  1:1239790/93‑1 (MQ=255)
                                                                                      gAACGTTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCCCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGGGGCATCAATACg   >  2:2682656/1‑93 (MQ=255)
                                                                                       aaCGTTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATACgg  <  1:1604111/93‑1 (MQ=255)
                                                                                       aaCGTTACGCCTTCAATATGGAACGCTTGCGGCTCGCCCGCGCGGACCGTCCCGCGGTCCCCCCTTCCTTTCTGCGCGGGGGCATCAATAAgg  >  2:215767/1‑93 (MQ=255)
                                                                                          |                                                                                          
AGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCTGCGCGGTGACATCAATACGG  >  NC_000913/3159391‑3159570

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