Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913603,3950CG21.2% 330.6 / 86.8 138intergenic (+60/+21)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (52/56);  new base G (15/14);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 146 total aligned reads displayed.

ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494
                                                                                                   |                                                                                                   
aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                     >  2:1387979/1‑100 (MQ=255)
aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                     >  2:2913485/1‑100 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  <  1:1123313/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACGTCCGACAAAACGATGTCAACCATCCGAAACCGCTCGCAGCCATTCGATg                                                                                                 <  2:2833779/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  1:2634933/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATACATTc                                                                                                     <  2:1358669/96‑1 (MQ=255)
     tGTTTATGTCATTTAAAACGCTGCGTCGGAAATAAGGCATGCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                                   <  2:1497621/97‑1 (MQ=255)
      gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                               >  1:926067/1‑100 (MQ=255)
        ttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag                                                                                             <  1:2567648/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  1:2557070/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  1:2846060/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  1:1185233/100‑1 (MQ=39)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  2:474186/100‑1 (MQ=38)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  1:1669376/100‑1 (MQ=255)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                          >  1:1510576/1‑100 (MQ=255)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                          >  1:118972/1‑99 (MQ=255)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                         >  2:465280/1‑100 (MQ=255)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                         >  1:949233/1‑100 (MQ=255)
             gCATTTAAAACGCTGCGTCGGAAATAAGTCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                          <  2:1102981/98‑1 (MQ=38)
               aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                      >  1:605213/1‑100 (MQ=255)
                tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttcg                                                                                     <  1:883895/100‑3 (MQ=37)
                 ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                    >  2:1284229/1‑100 (MQ=255)
                 ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                    >  1:2413761/1‑100 (MQ=255)
                  tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                   >  1:1129726/1‑100 (MQ=255)
                     aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttcggatgg                                                                                <  1:2953560/100‑8 (MQ=25)
                        gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTan                                                                             <  2:1033320/100‑2 (MQ=255)
                         cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                            >  1:227851/1‑100 (MQ=255)
                         cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                            <  1:1147027/100‑1 (MQ=255)
                             gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                        >  2:269438/1‑100 (MQ=255)
                             gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                        >  1:73612/1‑100 (MQ=255)
                              tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTa                                                                       <  1:1968813/100‑1 (MQ=255)
                              tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTa                                                                       >  1:647964/1‑100 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                     >  1:578325/1‑100 (MQ=255)
                                 gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                    >  2:330827/1‑100 (MQ=255)
                                  aactAAGGCATACACGCTACATCCGACATAACGACGTCAACCTTCCCAAACCGCTCTCATTCATTAGATGAGTGCGATTTTTTTAATTACTCCTTAAATg                                                                   <  2:607356/97‑1 (MQ=255)
                                  aaaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                   <  2:2829459/100‑1 (MQ=255)
                                      aaGGCATTCACGCTACATCCGACAAAACGATTTCAACCATCACAAACCACTCTCATCCATTCGTTGAGAGCGGATTTTTTAATTACTGCTTAAATTCAc                                                                >  2:1179929/1‑99 (MQ=255)
                                      aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                               <  2:949224/100‑1 (MQ=255)
                                      aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                               >  2:2361840/1‑100 (MQ=255)
                                          cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                           <  2:60800/100‑1 (MQ=255)
                                            ttCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCag                                                         <  1:2933647/100‑1 (MQ=255)
                                            ttCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCAGTTGTTTTAATTACTGCTTAAATGCACCCGCCa                                                          >  2:2351812/1‑99 (MQ=255)
                                              cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                       <  1:628862/100‑1 (MQ=255)
                                                 gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                    <  2:368017/100‑1 (MQ=255)
                                                    aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                                 <  2:2086871/100‑1 (MQ=255)
                                                      aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                               >  1:2760979/1‑100 (MQ=255)
                                                        ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                             >  1:1575608/1‑100 (MQ=255)
                                                        ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                             >  2:702237/1‑100 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            <  1:269441/100‑1 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            <  1:2070357/100‑1 (MQ=255)
                                                              aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                       >  2:1025104/1‑100 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  2:2963161/100‑1 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  1:2913511/100‑1 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  2:1781169/100‑1 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCACTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  1:1757741/100‑1 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCACTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  1:1757777/100‑1 (MQ=255)
                                                                  gATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                   <  1:3027296/100‑1 (MQ=255)
                                                                    tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTgg                                 <  2:2886005/100‑1 (MQ=255)
                                                                      tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                               <  1:2436120/100‑1 (MQ=255)
                                                                      tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                               <  2:926059/100‑1 (MQ=255)
                                                                      tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                               <  1:2436392/100‑1 (MQ=255)
                                                                        aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                             <  1:1243259/100‑1 (MQ=255)
                                                                         cgatctcGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            <  1:1462446/94‑1 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            <  1:521556/100‑1 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  1:2793849/1‑100 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  2:2294227/1‑100 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  1:1962244/1‑100 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  2:1214568/1‑100 (MQ=255)
                                                                          ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                           <  2:118971/100‑1 (MQ=255)
                                                                          attaaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                           >  1:1447735/8‑100 (MQ=25)
                                                                            taaaaaaaCCGCTTTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAGATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                          <  2:1529963/94‑1 (MQ=255)
                                                                             tCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                        <  1:1146748/100‑1 (MQ=255)
                                                                             aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                        >  2:1109090/5‑100 (MQ=37)
                                                                              ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACttcatc                       <  2:2680034/100‑6 (MQ=255)
                                                                              aaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                       >  1:513895/4‑100 (MQ=37)
                                                                               cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                      <  2:2025866/100‑1 (MQ=255)
                                                                                gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACttcatcag                     >  1:2021769/1‑93 (MQ=255)
                                                                                gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                     <  2:241840/100‑1 (MQ=255)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:1017404/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACCGATGAAGTCGc                 >  2:2608363/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 <  1:2016744/100‑1 (MQ=255)
                                                                                     cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  1:2068733/1‑100 (MQ=38)
                                                                                     cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  2:2019750/1‑100 (MQ=39)
                                                                                     cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  1:878272/1‑100 (MQ=255)
                                                                                     cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  2:1158289/1‑100 (MQ=255)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:2508497/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:2865755/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:1406166/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:2281939/1‑100 (MQ=39)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:1032445/1‑100 (MQ=255)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCACCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               <  1:2021080/100‑1 (MQ=255)
                                                                                       ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa              >  2:2606980/1‑100 (MQ=39)
                                                                                          tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa           <  1:1695783/100‑1 (MQ=255)
                                                                                          tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa           >  2:680566/1‑100 (MQ=255)
                                                                                           cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa          >  1:1064929/1‑100 (MQ=255)
                                                                                           cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGCTGGCGACTGATGAAGTCGCAGACCaa          <  1:1847367/100‑1 (MQ=255)
                                                                                             tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc        >  1:1262140/1‑100 (MQ=255)
                                                                                              ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAct       >  1:3030792/1‑100 (MQ=255)
                                                                                                aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct     >  1:1896604/1‑100 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGTCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  2:1866759/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494

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