Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,732,563 (C)8→9 coding (492/732 nt) yfiH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,732,5551.C94.9% 404.7 / 13.6 137coding (500/732 nt)yfiHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base C (68/62);  total (72/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 184 total aligned reads displayed.

CAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCG  >  NC_000913/2732463‑2732642
                                                                                               |                                                                                        
cAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                           >  1:1991547/1‑93 (MQ=255)
 aaGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTc                                                                                         >  2:2619105/1‑93 (MQ=255)
 aaGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTc                                                                                         <  2:3314883/93‑1 (MQ=255)
 aaGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTc                                                                                         >  1:2309483/1‑93 (MQ=255)
  aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                        >  2:3882063/1‑93 (MQ=255)
   gCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTccc                                                                                       <  2:204733/93‑1 (MQ=255)
     tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTccccc                                                                                    >  1:1157212/1‑93 (MQ=255)
      gATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccc                                                                                   >  1:1524743/1‑93 (MQ=255)
      gATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccc                                                                                   >  2:591345/1‑93 (MQ=255)
      gATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccc                                                                                   >  1:3088928/1‑93 (MQ=255)
        tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccccc                                                                                 >  2:3713505/1‑93 (MQ=255)
        tAAATATCCCCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccccc                                                                                 >  1:361209/1‑93 (MQ=255)
            tatCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccccccact                                                                             >  1:4406700/1‑89 (MQ=255)
                 gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                         <  2:3488697/93‑1 (MQ=255)
                 gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                         <  1:4549725/93‑1 (MQ=255)
                  ccAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                        <  1:4336408/93‑1 (MQ=255)
                    aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgc                                                                      >  1:4472564/1‑93 (MQ=255)
                      atccTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGACCCTCCCCCCCCCACTTCGAAcgcgc                                                                    <  2:3121453/90‑1 (MQ=255)
                      aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcgc                                                                    <  1:758391/93‑1 (MQ=255)
                       tACTTATCCCCATGCTGAATTAAAGCTTCACCTGCTTTAGCGGCTACTGCCATAAACGCCCCGCGAACCCCCCCCCCCCACTTCGAAcgcgcg                                                                   <  1:3152889/93‑1 (MQ=255)
                       tACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcgcg                                                                   >  1:1294884/1‑93 (MQ=255)
                        aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGt                                                                  >  1:4631660/1‑93 (MQ=255)
                        aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGt                                                                  <  1:1395675/93‑1 (MQ=255)
                          ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTgg                                                                >  1:2071327/1‑93 (MQ=255)
                          ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTgg                                                                <  2:4622568/93‑1 (MQ=255)
                            aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGAc                                                              <  2:4489230/93‑1 (MQ=255)
                            aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGAc                                                              >  1:3891199/1‑93 (MQ=255)
                             tCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTACCCCCCCCACTTCCACCGCGCGTCGAcc                                                             >  2:1300910/1‑93 (MQ=255)
                               aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCaa                                                           >  1:4361642/1‑93 (MQ=255)
                               aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCaa                                                           >  1:394382/1‑93 (MQ=255)
                               aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAACCGCGTGGACCaa                                                           >  1:293639/1‑93 (MQ=255)
                                ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGCACGCGCGTGGACCAAt                                                          >  2:4116410/1‑93 (MQ=255)
                                ccATGCTGAATGAAAGCGCCACTTGCTTTAGCGTCTACTGCCATAAACGCCCCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAAt                                                          <  2:1423478/93‑1 (MQ=255)
                                 cATGCTGAATGAAAGCTGCACTTGCTTTGGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                         <  1:4529813/93‑1 (MQ=255)
                                 cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                         <  2:911983/93‑1 (MQ=255)
                                  aTGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGGGCGTGGACAAAttc                                                        >  1:3268230/1‑92 (MQ=255)
                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                       <  1:3709715/93‑1 (MQ=255)
                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                       <  1:1316257/93‑1 (MQ=255)
                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGcc                                                      >  2:4531530/1‑93 (MQ=255)
                                    gCTGAATGAAAGCTGCACTTGCTTTAGCCTCTACTGCCATAAACGCCTCGCGAGCCTCCCCGCCCCACTTCGAACGCGCGTGGACCAATTGcc                                                      <  2:3755203/93‑1 (MQ=255)
                                      tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCgg                                                    >  2:661041/1‑93 (MQ=255)
                                       gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                   >  1:594924/1‑93 (MQ=255)
                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                  >  2:2163411/1‑93 (MQ=255)
                                          tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCt                                                <  2:4595847/93‑1 (MQ=255)
                                          tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCt                                                >  2:771022/1‑93 (MQ=255)
                                          tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACCTCGAACGCGCGTGGACCAATTGCCCGCCCt                                                >  1:2854270/1‑93 (MQ=255)
                                          tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTTCAAACCGCGTGGGCCAAATGCCCGCCCTa                                               >  1:4369939/1‑93 (MQ=255)
                                           gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTa                                               >  2:1244453/1‑93 (MQ=255)
                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAACCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                             <  1:1739618/93‑1 (MQ=255)
                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                             >  2:3377950/1‑93 (MQ=255)
                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                             >  2:3808568/1‑93 (MQ=255)
                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTAGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                             >  2:3080450/1‑93 (MQ=255)
                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                            <  1:1170224/93‑1 (MQ=255)
                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                            >  1:2626184/1‑93 (MQ=255)
                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                            <  2:1903964/93‑1 (MQ=255)
                                               gCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                           <  1:2620892/93‑1 (MQ=255)
                                                cTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                          <  2:3203773/93‑1 (MQ=255)
                                                cTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                          <  1:638937/93‑1 (MQ=255)
                                                cTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                          <  2:4481631/93‑1 (MQ=255)
                                                cTGCACGTGCTTTAGCTTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                          <  1:3458415/93‑1 (MQ=255)
                                                  tcACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGc                                        <  2:2820305/92‑1 (MQ=255)
                                                  gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACCTCCAACCCGGGGGGACCAAATGGCGGCCCCAAACAggg                                        >  1:1057685/1‑92 (MQ=255)
                                                    aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                      >  1:3367649/1‑93 (MQ=255)
                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgag                                     >  1:251638/1‑93 (MQ=255)
                                                      ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                    >  1:1592389/1‑93 (MQ=255)
                                                      ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                    >  2:1044861/1‑93 (MQ=255)
                                                      ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                    >  1:3318962/1‑93 (MQ=255)
                                                      ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCAAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                    >  1:2647014/1‑93 (MQ=255)
                                                       tGCTTTAGCGTCTCCTGCCATAAACGCCTCGCGAACCCCCCCCCCCCACTACGAACGCGCGTGGGCCAATTGCCGGCCCTAACCAGGCGAGaa                                   <  1:4230694/93‑1 (MQ=255)
                                                       tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                   <  2:2434179/93‑1 (MQ=255)
                                                       tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                   <  2:1886240/93‑1 (MQ=255)
                                                        gCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAata                                 >  1:3694318/1‑93 (MQ=255)
                                                          tttAGCGTCTACTGCCATAAACGCCTCGCGACCCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAatat                                <  2:1520564/93‑1 (MQ=255)
                                                          tttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAaaat                                >  1:2415982/1‑93 (MQ=255)
                                                              gCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTc                            >  2:943063/1‑93 (MQ=255)
                                                              gCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTc                            <  2:2111778/93‑1 (MQ=255)
                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                           <  2:4328085/93‑1 (MQ=255)
                                                                 tCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                         >  1:2247322/1‑93 (MQ=255)
                                                                 tCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCCAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcctg                         >  2:1733295/1‑91 (MQ=255)
                                                                     cTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTa                     <  1:805212/93‑1 (MQ=255)
                                                                      tGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAt                    >  2:2068674/1‑93 (MQ=255)
                                                                      tGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAt                    >  2:606922/1‑93 (MQ=255)
                                                                       gCCATAAACGCCTCGCGATCCCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATc                   <  1:4488871/93‑1 (MQ=255)
                                                                        ccATAAACGCCTCGCGACCCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                  <  1:2104877/93‑1 (MQ=255)
                                                                        ccATAAACGCCTCGCGACCCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                  <  2:5484/93‑1 (MQ=255)
                                                                         cATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                 >  2:4063764/1‑93 (MQ=255)
                                                                         cATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                 <  1:969813/93‑1 (MQ=255)
                                                                          aTAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAga                >  2:1355940/1‑92 (MQ=255)
                                                                           tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTCACCAGGCGAGAATATTTTTCGTATTATCAGCa               >  1:4281383/1‑93 (MQ=255)
                                                                           tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGGGGACCAATTGCCGGCCCTAACCAGCCGAGAATATTTTCGGGATTATCAGCa               >  2:1658008/1‑93 (MQ=255)
                                                                            aaaCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa              <  2:2979560/93‑1 (MQ=255)
                                                                             aaCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa             >  2:3239337/1‑93 (MQ=255)
                                                                                 ccTCGCGACCCCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg         <  1:933363/93‑1 (MQ=255)
                                                                                 ccTCGCGAGCCTCCCCCCCCCGCTGCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg         <  2:2371602/93‑1 (MQ=255)
                                                                                  cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg        <  1:3721285/93‑1 (MQ=255)
                                                                                   tCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa              >  1:2989449/1‑86 (MQ=255)
                                                                                   tCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa       >  2:1681830/1‑93 (MQ=255)
                                                                                       gAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAcc   <  1:1243531/93‑1 (MQ=255)
                                                                                        aaCCTCCCCCCCCCACTTCGACCGCGCCTGTACCAATTTACGGCCCTAACCAGGCGATAATATTTTTCGCATTATCCGCAAAACAGGCAAcct  >  2:730819/1‑92 (MQ=255)
                                                                                        aaCCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATTAGCAAAACAGGAAACCg  >  2:1445621/1‑93 (MQ=255)
                                                                                               |                                                                                        
CAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCG  >  NC_000913/2732463‑2732642

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