Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 926,180 A→G 15.0% intergenic (‑209/+45) serW ← / ← infA tRNA‑Ser/translation initiation factor IF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913926,1800AG15.0% 452.8 / 67.0 199intergenic (‑209/+45)serW/infAtRNA‑Ser/translation initiation factor IF‑1
Reads supporting (aligned to +/- strand):  ref base A (74/94);  new base G (14/16);  total (89/110)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 201 total aligned reads displayed.

AACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCA  >  NC_000913/926083‑926279
                                                                                                 |                                                                                                   
aaCAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAACGGCCGGTTAAACCGAcc                                                                                                   >  2:786004/1‑100 (MQ=255)
aaCAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGAcc                                                                                                   <  2:887823/100‑1 (MQ=255)
       tgtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACt                                                                                            <  1:733225/100‑1 (MQ=255)
        gtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttatt                                                                                           <  1:3883306/100‑5 (MQ=38)
        gtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTc                                                                                           >  1:2444266/1‑100 (MQ=255)
        gtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTc                                                                                           >  2:3307464/1‑100 (MQ=255)
         tgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattt                                                                                          <  1:2884335/100‑6 (MQ=38)
         tgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCg                                                                                          <  1:1418947/100‑1 (MQ=255)
          gCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttatttt                                                                                         <  2:4085556/100‑7 (MQ=37)
          gCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttatttt                                                                                         <  1:2224066/100‑7 (MQ=37)
          gCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGt                                                                                         >  2:2524177/1‑100 (MQ=255)
            aGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttgt                                                                                       <  2:2684191/100‑9 (MQ=37)
             gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttct                                                                                      >  1:4120549/1‑100 (MQ=255)
             gTATATCTACATCGAGACAAGTTACGGACTTAGACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttgt                                                                                       <  1:1845363/99‑9 (MQ=25)
                tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                                   <  1:1711436/100‑1 (MQ=255)
                tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                                   <  2:2945557/100‑1 (MQ=255)
                tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                                   <  1:2090711/100‑1 (MQ=255)
                tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                                   <  1:1028854/100‑1 (MQ=255)
                tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                                   <  1:1613309/100‑1 (MQ=255)
                tatCTACATCGAGACAAGTTACAGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctt                                                                                     <  2:548433/98‑1 (MQ=255)
                 atCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttct                                                                                  >  1:2054014/1‑100 (MQ=255)
                  tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAATGCCGGTTAAACCGACCTTTTACTCGTTCTTtctc                                                                                 <  2:1278384/100‑1 (MQ=255)
                  tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTtctc                                                                                 <  2:3056819/100‑1 (MQ=255)
                  tCTACATCGAGACAAGTTACGGACTCATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAACCCGCCCTTTTACTCGTTCTTtatc                                                                                 >  2:142781/1‑100 (MQ=255)
                   cTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCGCAAAATAAAACCACCGGATAAACCGACCTTTTACTAGTTCTTTctct                                                                                >  2:2624493/1‑100 (MQ=255)
                     aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTc                                                                              <  1:3379422/100‑1 (MQ=255)
                      cATCTAGACAAGTTATGGACTTATACTTCGAAAGTACATCATACATACCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCg                                                                             <  2:1900369/100‑1 (MQ=255)
                       aTCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGc                                                                            <  2:3907420/100‑1 (MQ=255)
                          gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                         >  2:1197519/1‑100 (MQ=255)
                          gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                         <  2:124945/100‑1 (MQ=255)
                          gagaCAAGTTACGGACTTATACTTACAAAGTACTTTATACATATCACAAAATAAAAAGGACGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                         >  2:44909/1‑100 (MQ=255)
                            gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATc                                                                       <  1:2193100/100‑1 (MQ=255)
                             aCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCa                                                                      <  2:1078259/100‑1 (MQ=255)
                              cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAg                                                                     >  2:2536145/1‑100 (MQ=255)
                               aaGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATACAAACGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAgg                                                                    <  2:4086545/100‑1 (MQ=255)
                               aaGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAgg                                                                    >  1:1858228/1‑100 (MQ=255)
                                aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                                   >  1:1905490/1‑100 (MQ=255)
                                aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                                   <  2:2135727/100‑1 (MQ=255)
                                 gTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCg                                                                  >  1:3269251/1‑100 (MQ=255)
                                 gTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCg                                                                  <  1:3533345/100‑1 (MQ=255)
                                  ttACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                 <  1:996447/100‑1 (MQ=255)
                                  ttACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                 <  1:2158425/100‑1 (MQ=255)
                                  ttACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                 <  2:3671200/100‑1 (MQ=255)
                                   tACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGt                                                                <  2:1734510/100‑1 (MQ=255)
                                      ggACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa                                                             <  2:3841860/100‑1 (MQ=255)
                                        aCTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAAc                                                           <  2:525391/100‑1 (MQ=255)
                                          ttATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                         <  1:3605180/100‑1 (MQ=255)
                                          ttATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                         >  1:3560072/1‑100 (MQ=255)
                                            ntaCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc                                                       >  2:3417602/2‑100 (MQ=255)
                                             taCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                      <  2:420107/100‑1 (MQ=255)
                                             taCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                      >  1:1606454/1‑100 (MQ=255)
                                             taCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                      <  1:462187/100‑1 (MQ=255)
                                               cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                                    >  1:2661425/1‑100 (MQ=255)
                                               cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                                    >  1:1844654/1‑100 (MQ=255)
                                               cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                                    >  1:2361481/1‑100 (MQ=255)
                                                ttCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCg                                                   <  2:2317629/100‑1 (MQ=255)
                                                  ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGAc                                                 >  1:627017/1‑100 (MQ=255)
                                                  ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGAc                                                 >  2:2865353/1‑100 (MQ=255)
                                                  ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGAc                                                 <  1:3159113/100‑1 (MQ=255)
                                                        tatttCATACATATTACATAATAAAAAGTCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTTGTCCATCAGGCGGTAAAACAATCAGCGACTACg                                             <  2:68224/95‑1 (MQ=255)
                                                        tACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                           <  2:3255354/100‑1 (MQ=255)
                                                         aCTTCATACATATCTCCAAATAAAAAGGCCGGTTAAATCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                           <  2:1350745/99‑1 (MQ=255)
                                                           ttCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGa                                        <  1:4055266/100‑1 (MQ=255)
                                                            tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                                       <  1:653061/100‑1 (MQ=255)
                                                            tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                                       >  1:1926615/1‑100 (MQ=255)
                                                             catacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACa                                      <  1:237880/100‑1 (MQ=255)
                                                               tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAAt                                    >  1:1423227/1‑100 (MQ=255)
                                                                acataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg                                   >  2:3519028/1‑100 (MQ=255)
                                                                 cataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGc                                  >  1:224539/1‑100 (MQ=255)
                                                                    aaaaaaaaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                               <  2:324036/95‑1 (MQ=255)
                                                                     tCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGcc                              >  2:1048699/1‑100 (MQ=255)
                                                                      cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                             <  1:1741317/100‑1 (MQ=255)
                                                                      cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                             >  2:1880367/1‑100 (MQ=255)
                                                                       acaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                            >  1:1142315/1‑100 (MQ=255)
                                                                       acaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                            <  1:2336502/100‑1 (MQ=255)
                                                                        aacgagtaaaaGGTCGGTTTAAACGGCTTTATACTCGTTCTTTCTCTTCGCCAATAAGGCGGTTAAACAATACGCAACTACGGAAAACACTGCGGccat                            >  2:594157/8‑97 (MQ=255)
                                                                         aaaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTg                          <  1:1566772/100‑1 (MQ=255)
                                                                          cgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                         >  2:504148/6‑100 (MQ=38)
                                                                          cgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                         >  1:3742632/6‑100 (MQ=38)
                                                                           gagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCt                        >  1:3011641/5‑100 (MQ=38)
                                                                           aaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGCAAGACAATGCGGCCTTTGCt                        >  2:2000985/1‑100 (MQ=255)
                                                                            aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc                       >  1:3349497/1‑100 (MQ=255)
                                                                             tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa                      <  1:3770502/100‑1 (MQ=255)
                                                                             tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa                      <  1:530611/100‑1 (MQ=255)
                                                                             gtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa                      >  1:3119272/3‑100 (MQ=39)
                                                                              taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGCCTACGGAAGACAATGCGGCCTTTGCTCAg                     >  2:806736/2‑100 (MQ=39)
                                                                               aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg                    >  2:4021703/1‑100 (MQ=255)
                                                                                aaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTATGGAAGACAATGCGGCCTTTGCTCAGGt                   >  1:1244510/1‑100 (MQ=255)
                                                                                aaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGt                   >  2:3072657/1‑100 (MQ=255)
                                                                                aaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGt                   <  1:2987138/100‑1 (MQ=255)
                                                                                 aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTc                  >  1:2148207/1‑100 (MQ=255)
                                                                                  aGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTc                  <  2:3122207/99‑1 (MQ=255)
                                                                                  aGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTc                  <  1:304218/99‑1 (MQ=255)
                                                                                  aGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCg                 >  1:3965090/1‑100 (MQ=255)
                                                                                   ggCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGt                >  2:2364806/1‑100 (MQ=255)
                                                                                       ggTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg            >  1:2701676/1‑100 (MQ=255)
                                                                                           taaCCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTc        <  1:504152/99‑1 (MQ=255)
                                                                                           aaaCCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTc        <  2:1578666/100‑1 (MQ=255)
                                                                                              ccGACCGGTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGCTCGTACGGGGTCAg      <  2:2439964/99‑1 (MQ=255)
                                                                                                gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTc   <  2:1463802/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
AACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCA  >  NC_000913/926083‑926279

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