Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5240CT48.5% 420.7 / inf 375I156I (ATC→ATT) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base C (92/101);  new base T (89/93);  total (181/194)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 460 total aligned reads displayed.

ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGC  >  NC_000913/580438‑580624
                                                                                      |                                                                                                    
              tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                                          <  1:9171237/97‑1 (MQ=0)
               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                                         <  1:2445343/98‑1 (MQ=0)
               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                                         >  2:6090144/4‑101 (MQ=0)
                                  gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                                      <  1:4220264/96‑1 (MQ=0)
                                   cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                                     <  1:7755946/97‑1 (MQ=0)
                                   cttcccGGAGCTGGAAAACCGACATATGGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                                     >  1:11155711/5‑101 (MQ=0)
                                    ttcccGGAGCTGGAAAACCGACATATTGATTTCATCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCt                                                    <  2:3122214/98‑1 (MQ=0)
                                                gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                                        >  1:5071534/1‑101 (MQ=11)
                                                gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                                        <  1:4652625/101‑1 (MQ=11)
                                                gAAAACCGACATATTGATTTCCTCAAGAAGGACCTCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                                        <  1:5189711/101‑1 (MQ=11)
                                                 aaaaCCGACATATTGTTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                                       >  1:13383260/1‑101 (MQ=11)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                                      <  1:10068588/101‑1 (MQ=11)
                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                                      >  1:2120361/1‑101 (MQ=11)
                                                   aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatatnnn                                     >  2:14200387/1‑98 (MQ=11)
                                                   aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                                     <  1:9945328/101‑1 (MQ=33)
                                                   aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGTTGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                                     >  1:3368707/1‑101 (MQ=11)
                                                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATcnnn                                    >  2:13478024/1‑98 (MQ=14)
                                                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                                    <  2:1797250/101‑1 (MQ=35)
                                                     ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCTCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                                   <  1:13240345/101‑1 (MQ=25)
                                                     ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                                   >  2:6380918/1‑101 (MQ=35)
                                                     ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc                                   >  1:12422827/1‑101 (MQ=11)
                                                      nnncATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                                  <  2:12681847/98‑1 (MQ=38)
                                                      nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaaaa                                  <  2:14720517/98‑3 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                                  >  1:5826556/1‑101 (MQ=11)
                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCGGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                                  >  1:3691924/1‑101 (MQ=11)
                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaannn                                 >  2:9043237/1‑98 (MQ=25)
                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                                 >  1:10930731/1‑101 (MQ=38)
                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                                 <  1:10299462/101‑1 (MQ=38)
                                                       gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                                 >  1:1732193/1‑101 (MQ=14)
                                                        gcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                                <  2:5743167/100‑1 (MQ=14)
                                                        aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGAACAGt                                >  1:13112047/1‑101 (MQ=255)
                                                        aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAcnnn                                >  2:9195353/1‑98 (MQ=2)
                                                        aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                                >  1:344436/1‑101 (MQ=9)
                                                        aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                                >  1:1255696/1‑101 (MQ=14)
                                                         nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                               <  2:10697805/98‑1 (MQ=255)
                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCGGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                               >  1:7001120/1‑101 (MQ=255)
                                                         cATATTGATTTCGTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACannn                               >  2:7616095/1‑98 (MQ=11)
                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                              >  1:5638108/4‑101 (MQ=14)
                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                              <  2:6846006/98‑1 (MQ=14)
                                                          aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                              <  1:13455354/101‑1 (MQ=255)
                                                          aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                              >  2:3894571/1‑101 (MQ=255)
                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                             <  1:13539510/99‑1 (MQ=14)
                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                             <  1:1819050/99‑1 (MQ=14)
                                                           nnntGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                             <  2:14461727/98‑1 (MQ=14)
                                                            nnngATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCatg                            <  2:10967530/98‑3 (MQ=14)
                                                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                            <  1:6277127/101‑1 (MQ=255)
                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgt                            >  2:3342580/2‑100 (MQ=14)
                                                             ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt                           <  1:10560909/101‑1 (MQ=14)
                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa                         <  2:1388016/101‑1 (MQ=17)
                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGACCAGTCAGGTTaa                        <  2:4648692/101‑1 (MQ=255)
                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn                        >  2:11991787/1‑98 (MQ=255)
                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa                        >  1:12663679/1‑101 (MQ=255)
                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc                       >  2:953008/1‑101 (MQ=255)
                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc                       <  1:8516342/101‑1 (MQ=255)
                                                                  ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa                      <  2:5544011/101‑1 (MQ=255)
                                                                  nnncTGAAACGGGATATCATCAAAGCCACGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa                      <  2:10390509/98‑1 (MQ=255)
                                                                   nnntGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAccg                     <  2:12433543/98‑3 (MQ=255)
                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAcnnn                    >  2:8452882/1‑98 (MQ=255)
                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg                    <  1:14556728/101‑1 (MQ=255)
                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGCGAATATATCGAACAGTCAGGGTAACAgg                    >  1:10155597/1‑101 (MQ=255)
                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc                   <  1:8842837/101‑1 (MQ=255)
                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc                   <  1:14188602/101‑1 (MQ=255)
                                                                       gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg                 <  1:4848505/101‑1 (MQ=255)
                                                                        nnncGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc                <  2:11561535/98‑1 (MQ=255)
                                                                        aagaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGc                >  1:11461133/6‑101 (MQ=17)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc                <  1:9275113/101‑1 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc                <  1:8211750/101‑1 (MQ=255)
                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc                >  1:732498/1‑101 (MQ=255)
                                                                         agaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGca               >  1:6022764/5‑100 (MQ=17)
                                                                         agaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGca               <  1:4357695/97‑2 (MQ=17)
                                                                          nnnggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCgg              <  2:8233155/98‑1 (MQ=255)
                                                                          gaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag              >  1:306071/4‑99 (MQ=17)
                                                                          gaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag              >  1:11803802/4‑99 (MQ=17)
                                                                          gaaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcag              >  1:14741330/4‑99 (MQ=17)
                                                                           gaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc             <  1:5490938/99‑1 (MQ=17)
                                                                           aaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc             >  1:1215277/3‑101 (MQ=17)
                                                                           aaggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGc             >  2:4272701/3‑101 (MQ=17)
                                                                            aggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagca            >  1:12979216/2‑101 (MQ=17)
                                                                            aggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagca            <  1:7193501/100‑1 (MQ=17)
                                                                             ggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCggnnn           >  2:13638499/1‑98 (MQ=255)
                                                                             ggATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAt           <  1:13616201/101‑1 (MQ=255)
                                                                             ggACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcac           >  2:6490755/1‑100 (MQ=17)
                                                                              gacatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcact          <  1:2603328/98‑3 (MQ=17)
                                                                              gATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAtt          <  2:5352606/101‑1 (MQ=255)
                                                                              gATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAtt          <  1:14431444/101‑1 (MQ=255)
                                                                               acatcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactc         >  1:2978176/3‑98 (MQ=17)
                                                                                tatCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAtttt        >  2:500185/1‑101 (MQ=255)
                                                                                tatCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCAtttt        <  2:2776332/101‑1 (MQ=255)
                                                                                catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactct        >  1:5863190/2‑97 (MQ=17)
                                                                                catcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactct        >  1:13912831/2‑97 (MQ=17)
                                                                                 atcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGTGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctg       >  1:4482491/1‑96 (MQ=17)
                                                                                 atcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctg       >  1:13979375/1‑96 (MQ=17)
                                                                                 atcattAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctg       >  1:8127765/1‑96 (MQ=17)
                                                                                 atcatcAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTg       <  1:9507500/101‑1 (MQ=255)
                                                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGTGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGt      <  1:659651/101‑1 (MQ=255)
                                                                                  nnntcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGt      <  2:7441012/98‑1 (MQ=255)
                                                                                   catcaAAGCCATGAACAAAGCAGTCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGTCAGGTTAACAGGCTGTGGCATTTTGTc     >  2:3614419/1‑101 (MQ=255)
                                                                                   catcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTc     >  2:2359934/1‑101 (MQ=255)
                                                                                     tcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCcg   >  1:9834424/1‑101 (MQ=255)
                                                                                      caAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCcgc  >  2:2104154/1‑101 (MQ=255)
                                                                                      |                                                                                                    
ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGC  >  NC_000913/580438‑580624

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