Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5120AG53.0% 138.4 / 223.4 151K152K (AAA→AAG) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (40/40);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 160 total aligned reads displayed.

CCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGG  >  NC_000913/580452‑580612
                                                            |                                                                                                    
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                              >  1:3861989/5‑101 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                              >  1:1336976/5‑101 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                              <  1:2814133/97‑1 (MQ=0)
tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTgg                                                              <  1:6384548/97‑1 (MQ=0)
 cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                             >  1:5618320/4‑101 (MQ=0)
                 ggggTTTTCCGGAACTAGAACACCGACATGTTGATTCCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCgggg                                             <  2:6658891/98‑1 (MQ=1)
                    gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                          <  1:1974867/96‑1 (MQ=0)
                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                         >  2:1701065/5‑101 (MQ=0)
                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCt                                        <  2:5190469/98‑1 (MQ=1)
                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGa                                  >  1:6859580/4‑101 (MQ=0)
                                  gAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                            <  2:459295/101‑1 (MQ=11)
                                  gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                            >  1:818442/1‑101 (MQ=11)
                                  gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                            >  2:2518279/1‑101 (MQ=11)
                                  gAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                            <  1:6094903/101‑1 (MQ=11)
                                   aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                           <  2:1988320/101‑1 (MQ=31)
                                   aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                           <  1:5963648/101‑1 (MQ=31)
                                   aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  1:1709959/101‑1 (MQ=11)
                                   aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  1:1636009/101‑1 (MQ=11)
                                   aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  2:2766750/101‑1 (MQ=11)
                                   aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           <  2:5428262/101‑1 (MQ=11)
                                   aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                           >  1:940585/1‑101 (MQ=11)
                                    aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                          <  1:723133/101‑1 (MQ=11)
                                    aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                          >  2:5350728/1‑101 (MQ=11)
                                    aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                          >  2:4697227/1‑101 (MQ=11)
                                     aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATTCCGGGGTTGCTGAGTGAATATATCGa                         <  2:407251/101‑1 (MQ=18)
                                     aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                         >  2:2706289/1‑101 (MQ=11)
                                     aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                         >  2:788575/1‑101 (MQ=11)
                                      aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                        >  1:5548232/1‑101 (MQ=35)
                                      aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                        <  2:2937587/101‑1 (MQ=35)
                                       ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                       <  1:2389707/101‑1 (MQ=35)
                                        cNACNNGTTGATTTNCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAct                      >  2:1124158/1‑100 (MQ=255)
                                        cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                      >  1:1385385/1‑101 (MQ=38)
                                        cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                      >  1:26144/1‑101 (MQ=11)
                                        cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                      <  1:5094170/101‑1 (MQ=11)
                                        cGACATATTGATTTACTCAAGAAGGACATCATTAAAGCCGAGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGAGTTGCTGAGTGAATATATCGAAct                      >  2:3314966/1‑100 (MQ=2)
                                         gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                     >  1:1637257/1‑101 (MQ=38)
                                         gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                     >  1:177507/1‑101 (MQ=14)
                                         gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCGGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg                     <  1:2495019/101‑1 (MQ=14)
                                          aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGGTGCTGAGAGAATATATCGAACAGt                    >  2:3107281/1‑101 (MQ=255)
                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                    >  2:5130946/1‑101 (MQ=14)
                                           cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                   <  2:2347671/101‑1 (MQ=40)
                                           cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                   <  2:2355169/101‑1 (MQ=40)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                   <  1:3518435/101‑1 (MQ=14)
                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                   <  2:598525/101‑1 (MQ=14)
                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                  >  1:4223862/4‑101 (MQ=14)
                                            aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                  >  1:6975672/1‑101 (MQ=255)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:592119/99‑1 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 >  2:3528753/3‑101 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:376161/99‑1 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  2:581850/99‑1 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 >  2:1315506/3‑101 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:2143100/99‑1 (MQ=14)
                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                 <  1:1100005/99‑1 (MQ=14)
                                             tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTCGCTGAGTGAATATATCGAACAGCCAg                 <  2:6925260/101‑1 (MQ=255)
                                              gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  2:6419177/101‑1 (MQ=255)
                                              gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  2:6471093/101‑1 (MQ=255)
                                              attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  2:3944054/100‑1 (MQ=14)
                                              attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                <  1:3911692/100‑1 (MQ=14)
                                               ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  1:1142733/101‑1 (MQ=255)
                                               ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  2:2631700/101‑1 (MQ=255)
                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGATCAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTTAATATATCGATCCGTCAGGt               <  2:5190463/101‑1 (MQ=9)
                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt               <  1:2582601/101‑1 (MQ=14)
                                                tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              >  2:2243521/1‑101 (MQ=255)
                                                tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              >  2:5006605/1‑101 (MQ=255)
                                                tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              <  1:5850247/101‑1 (MQ=255)
                                                tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              >  2:2308063/1‑101 (MQ=255)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              <  1:4905975/101‑1 (MQ=14)
                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGAGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt              >  2:1621258/1‑101 (MQ=14)
                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             <  1:3125118/101‑1 (MQ=255)
                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             <  1:2539804/101‑1 (MQ=255)
                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             >  1:1555563/1‑101 (MQ=255)
                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             >  2:2833556/1‑101 (MQ=255)
                                                 gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa             >  2:937329/1‑101 (MQ=17)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGACCAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  2:1041258/1‑101 (MQ=255)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  2:887219/1‑101 (MQ=255)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  1:5930770/1‑101 (MQ=255)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  1:4157475/1‑101 (MQ=255)
                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  1:2803826/1‑101 (MQ=255)
                                                  aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            <  2:2568269/101‑1 (MQ=17)
                                                  aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa            >  1:6221983/1‑101 (MQ=17)
                                                   tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc           <  1:3811359/101‑1 (MQ=255)
                                                   tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc           >  2:6240805/1‑101 (MQ=255)
                                                    ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa          >  1:5405178/1‑101 (MQ=255)
                                                    ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa          <  2:4898090/101‑1 (MQ=255)
                                                    tgccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa          <  2:2634664/99‑1 (MQ=255)
                                                     tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg         >  2:1992928/1‑101 (MQ=255)
                                                      ccTGAAACGGGATATCGTCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACAGATCCCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg        <  1:3341468/101‑1 (MQ=255)
                                                      ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg        <  1:4953126/101‑1 (MQ=255)
                                                      ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg        <  1:5071082/101‑1 (MQ=255)
                                                       cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc       <  2:5452777/101‑1 (MQ=255)
                                                       cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc       >  1:4337777/1‑101 (MQ=255)
                                                       cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc       >  2:2885619/1‑101 (MQ=255)
                                                        tGAAACGGGATATCATCAAAGCCATGAGCAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGCGTGAATATATCTAACAGTCAGGTTAACAgggt      <  2:6960285/101‑3 (MQ=17)
                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCt      >  1:3575998/1‑101 (MQ=255)
                                                         gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg     >  1:2090127/1‑101 (MQ=255)
                                                          aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc    >  2:6052403/1‑101 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg   <  1:3273531/101‑1 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg   <  1:3674642/101‑1 (MQ=255)
                                                           aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg   >  1:4574714/1‑101 (MQ=255)
                                                            aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTGAACAGGCTGCgg  <  2:6238118/101‑1 (MQ=255)
                                                            |                                                                                                    
CCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGG  >  NC_000913/580452‑580612

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