Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,8780TA21.6% 251.8 / 72.2 125intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (51/47);  new base A (12/15);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 128 total aligned reads displayed.

TGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAT  >  NC_000913/3477779‑3477977
                                                                                                   |                                                                                                    
tgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                      >  2:312181/1‑100 (MQ=255)
    gACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTccc                                                                                                  <  2:222766/100‑1 (MQ=255)
      ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCg                                                                                                >  1:2321231/1‑100 (MQ=255)
      ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGAGTAAAAAAAGCGGGGAATCCCCg                                                                                                <  2:348169/100‑1 (MQ=255)
       cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                               >  2:535051/1‑100 (MQ=255)
         agaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tt                                                                                            <  1:1812651/100‑1 (MQ=255)
           aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tttt                                                                                          <  1:1040342/100‑1 (MQ=255)
           aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tttt                                                                                          <  1:2963669/100‑1 (MQ=255)
            aGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttt                                                                                          <  2:1721335/99‑1 (MQ=255)
             gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tttttt                                                                                        >  1:1104928/1‑100 (MQ=255)
             gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tttttt                                                                                        >  2:553762/1‑100 (MQ=255)
             gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttt                                                                                        <  1:837949/100‑1 (MQ=255)
             gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttt                                                                                        <  2:1878452/100‑1 (MQ=255)
             atctaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑tttttt                                                                                        <  2:1731987/96‑1 (MQ=255)
               tAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGt                                                                                      >  2:1322717/1‑100 (MQ=255)
                aaaaGCCAGCGAGCCGTCGAACAGCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑ttttttt                                                                                       <  2:2673552/98‑2 (MQ=255)
                aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTc                                                                                     <  2:1745479/100‑1 (MQ=255)
                aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTc                                                                                     >  2:2879195/1‑100 (MQ=255)
                aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttttac                                                                                     <  1:2679826/100‑4 (MQ=255)
                aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttttac                                                                                     <  1:181989/100‑4 (MQ=255)
                aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑ttttttta                                                                                      <  2:522122/99‑3 (MQ=255)
                 aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑ttttttagat                                                                                    >  1:1732006/1‑96 (MQ=255)
                 aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCa                                                                                    >  2:728338/1‑100 (MQ=255)
                 aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttttacg                                                                                    <  1:1078973/100‑5 (MQ=255)
                 aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGC‑tttttttacg                                                                                    <  2:1676824/100‑5 (MQ=255)
                  aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGCAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCAc                                                                                   <  1:884074/100‑1 (MQ=255)
                  aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCAc                                                                                   >  2:526544/1‑100 (MQ=255)
                   agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACt                                                                                  >  1:2655876/1‑100 (MQ=255)
                   agccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACt                                                                                  >  2:1209418/1‑100 (MQ=255)
                    gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACTCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACtt                                                                                 >  1:1247209/1‑100 (MQ=255)
                        gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTc                                                                             <  1:1418171/100‑1 (MQ=255)
                         ccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCg                                                                            <  1:372239/100‑1 (MQ=255)
                          cagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCgg                                                                           >  2:2120378/1‑100 (MQ=255)
                          cagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCgg                                                                           <  2:2155296/100‑1 (MQ=255)
                          cagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCgg                                                                           <  1:2540794/100‑1 (MQ=255)
                           agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGt                                                                          <  1:337965/100‑1 (MQ=255)
                           agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGt                                                                          >  1:2470725/1‑100 (MQ=255)
                           agccGCCGAACATCGCGTCGCAGACTGAAGAAGCGCCACCGCCACATTATTGAGGCGTAAAAAAGGCGGGGATTCCCCGC‑TTTTTTGTCGCTTTTTCgg                                                                           <  2:878935/99‑1 (MQ=255)
                            gccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTa                                                                         <  2:2321208/100‑1 (MQ=255)
                             ccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTAt                                                                        >  2:1452567/1‑100 (MQ=255)
                              cggcGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTAtt                                                                       <  1:684172/97‑1 (MQ=255)
                               gTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTa                                                                      >  1:1258738/1‑100 (MQ=255)
                                 cGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAg                                                                     <  2:1247198/99‑1 (MQ=255)
                                   aaCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTgg                                                                  <  2:2060253/100‑1 (MQ=255)
                                   aaCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTgg                                                                  <  2:154358/100‑1 (MQ=255)
                                     cATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgca                                                                >  2:306876/1‑100 (MQ=255)
                                       tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaac                                                              <  1:720076/100‑1 (MQ=255)
                                       tCGCGTCGCAGGCTAAAGAAACGCCTCCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTGGTCACTTTTTCAGTATTAGTGgcaac                                                              <  2:2600439/100‑1 (MQ=255)
                                        cgcgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaacc                                                             <  1:863444/100‑1 (MQ=255)
                                        cgcgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaacc                                                             >  1:29144/1‑100 (MQ=255)
                                          cgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTCGTCACTTTTTCGGTATTAGTGgcaaccgc                                                           <  1:1749336/100‑1 (MQ=255)
                                            tCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgcaa                                                         <  2:958267/100‑1 (MQ=255)
                                             cGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgcaac                                                        >  1:2704145/1‑100 (MQ=255)
                                              gCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgcaacc                                                       <  1:2549350/100‑1 (MQ=255)
                                               cAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGcaaccgcaaccg                                                      >  1:2071322/1‑100 (MQ=255)
                                                  gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgccg                                                   >  2:397932/1‑100 (MQ=255)
                                                     gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTg                                                >  2:902259/1‑100 (MQ=255)
                                                     gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTg                                                >  1:2799962/1‑100 (MQ=255)
                                                         aaatgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCtt                                            <  1:1242486/96‑1 (MQ=255)
                                                         aaaCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCtt                                            >  1:2392155/1‑100 (MQ=255)
                                                          aaCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCttt                                           >  1:884030/1‑100 (MQ=255)
                                                           aCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTa                                          <  2:1258727/100‑1 (MQ=255)
                                                            cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                         >  2:2514185/1‑100 (MQ=255)
                                                             gccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAcc                                        <  2:1063470/100‑1 (MQ=255)
                                                              ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCg                                       <  2:3014155/100‑1 (MQ=255)
                                                              ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGCATCCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCg                                       <  1:1341876/100‑1 (MQ=255)
                                                               caccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                      >  2:2924656/1‑100 (MQ=255)
                                                               caccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                      >  2:715396/1‑100 (MQ=255)
                                                                 ccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                    <  1:315182/100‑1 (MQ=255)
                                                                  cgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                   <  1:1669641/100‑1 (MQ=255)
                                                                   gccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcaca                                  <  1:57763/100‑1 (MQ=255)
                                                                    ccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                 <  1:1557373/100‑1 (MQ=255)
                                                                       caTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagcag                              <  2:1326387/100‑1 (MQ=255)
                                                                         ttattGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGcc                            <  1:346521/100‑1 (MQ=255)
                                                                          tattGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCt                           >  2:1398893/1‑100 (MQ=255)
                                                                              gAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAAACGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGc                       >  2:162295/1‑100 (MQ=255)
                                                                                  tgacaaaaaaGCGGGGAATCCCCGCTTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAc                    >  1:2675850/5‑100 (MQ=38)
                                                                                  tgacaaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                   >  2:2304535/5‑100 (MQ=255)
                                                                                   gacaaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCg                  >  2:834082/4‑100 (MQ=255)
                                                                                   gTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCg                  >  1:914083/1‑100 (MQ=255)
                                                                                    tAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                 >  2:664420/1‑100 (MQ=255)
                                                                                     naaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTTGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtg                >  2:658508/2‑100 (MQ=255)
                                                                                     naaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtg                >  2:2026164/2‑100 (MQ=255)
                                                                                     caaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtg                >  1:712065/2‑100 (MQ=255)
                                                                                     caaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtg                >  1:1222988/2‑100 (MQ=255)
                                                                                     caaaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCATTTACTGCCACAGCAGCCTTCGCCACCgtg                >  2:1831133/2‑100 (MQ=39)
                                                                                     aaaaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtg                >  2:1813801/1‑100 (MQ=255)
                                                                                      aaaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtgt               <  1:118131/100‑1 (MQ=255)
                                                                                       aaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt              <  1:470654/100‑1 (MQ=255)
                                                                                       aaaaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt              >  2:653631/1‑100 (MQ=255)
                                                                                       aaaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt              >  1:2286961/1‑100 (MQ=255)
                                                                                        aaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCTTTTCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTc             >  2:253778/1‑100 (MQ=255)
                                                                                        aaaaGCGGGGAATCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTc             >  1:491959/1‑100 (MQ=255)
                                                                                         aaaGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCa            <  2:2799935/100‑1 (MQ=255)
                                                                                            gCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGa         >  2:166883/1‑100 (MQ=255)
                                                                                            gCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGa         >  1:1817869/1‑100 (MQ=255)
                                                                                              ggggATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCTTGTTCATGAcc       >  1:1328265/1‑100 (MQ=255)
                                                                                               gggATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCg      >  2:1502219/1‑100 (MQ=255)
                                                                                               gggATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCg      >  2:1502319/1‑100 (MQ=255)
                                                                                                   ttCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAACGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAt  <  1:2981943/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
TGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGC‑TTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAT  >  NC_000913/3477779‑3477977

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