Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 580,509 G→C L151L (CTG→CTC nohB → bacteriophage DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,5090GC75.9% 43.3 / 144.4 191L151L (CTG→CTCnohBbacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base C (68/77);  total (91/100)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 196 total aligned reads displayed.

GCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTA  >  NC_000913/580420‑580601
                                                                                         |                                                                                            
gccggTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGaaa                                                                                           <  2:4756334/90‑1 (MQ=1)
                                tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                           <  1:5503973/89‑1 (MQ=0)
                                tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                           <  1:524023/89‑1 (MQ=0)
                                tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                           >  1:4534568/5‑93 (MQ=0)
                                tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                           <  2:1025070/89‑1 (MQ=0)
                                tcctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagc                                                           >  2:580076/5‑93 (MQ=0)
                                 cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagct                                                          >  1:4733406/4‑92 (MQ=0)
                                 cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagct                                                          <  2:133955/90‑2 (MQ=0)
                                 cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagct                                                          <  2:5999696/90‑2 (MQ=0)
                                 cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcagct                                                          <  1:3429191/90‑2 (MQ=0)
                                          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                 >  2:2287180/4‑93 (MQ=0)
                                          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                 <  2:1521116/90‑1 (MQ=0)
                                          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                 <  2:4924481/90‑1 (MQ=0)
                                          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                 <  1:2047916/90‑1 (MQ=0)
                                          gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGAt                                                 >  2:606042/4‑93 (MQ=0)
                                                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                      >  1:2313508/5‑93 (MQ=0)
                                                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                      <  1:2725513/89‑1 (MQ=0)
                                                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                      <  1:5162767/89‑1 (MQ=0)
                                                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                      <  1:2916303/89‑1 (MQ=0)
                                                     cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc                                      <  1:3462938/89‑1 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  1:2229145/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  2:2861312/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     <  2:257586/90‑1 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  1:4682155/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  2:3574517/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  2:4148379/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     <  2:586894/90‑1 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     <  1:3083305/90‑1 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  2:5648296/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     <  2:4489409/90‑1 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCg                                     >  2:4629390/4‑93 (MQ=0)
                                                      ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAACTGCGCTGGATGAAATAATACCg                                     >  1:3483239/4‑93 (MQ=0)
                                                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                               >  1:1378002/4‑93 (MQ=0)
                                                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                               >  1:1000087/4‑93 (MQ=0)
                                                            gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                               >  2:5060495/4‑93 (MQ=0)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  2:3842000/1‑93 (MQ=2)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  1:89702/1‑93 (MQ=2)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  1:92400/1‑93 (MQ=2)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 <  2:620534/93‑1 (MQ=11)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  2:680571/1‑93 (MQ=2)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 <  2:1710688/93‑1 (MQ=11)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  1:602687/1‑93 (MQ=2)
                                                                          aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 >  1:5806432/1‑93 (MQ=2)
                                                                          aCATATTGATTGCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                 <  2:3368210/93‑1 (MQ=11)
                                                                           cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                >  2:4390664/1‑93 (MQ=32)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGTATATATc                >  1:787065/1‑93 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  2:2309651/1‑93 (MQ=2)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  2:847487/1‑93 (MQ=2)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                <  1:5619277/93‑1 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                <  1:4917677/93‑1 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  1:4363248/1‑93 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  1:3454639/1‑93 (MQ=2)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  1:2158758/1‑93 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                >  1:1337043/1‑93 (MQ=11)
                                                                           cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAAGACCGGGGCTGCTGAGGGAATATATc                >  1:4389072/1‑93 (MQ=2)
                                                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATca               >  1:2770265/4‑92 (MQ=2)
                                                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg               <  1:651362/90‑1 (MQ=11)
                                                                            atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg               <  2:250584/90‑1 (MQ=11)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa              >  2:3442629/3‑93 (MQ=2)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa              <  2:408470/91‑1 (MQ=11)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa              <  2:4273946/91‑1 (MQ=11)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa              >  1:64915/3‑93 (MQ=11)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa              >  1:6087248/3‑93 (MQ=11)
                                                                             tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGCGAATATATCGa              >  1:2589016/3‑93 (MQ=11)
                                                                             tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa              <  1:157302/93‑1 (MQ=34)
                                                                             tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa              >  1:364161/1‑93 (MQ=34)
                                                                              gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa             <  2:1915356/93‑1 (MQ=35)
                                                                              gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa             >  1:2022268/1‑93 (MQ=35)
                                                                               ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaaa            >  2:3277818/1‑92 (MQ=25)
                                                                               ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc            <  1:4119313/93‑1 (MQ=37)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            >  2:199999/1‑93 (MQ=11)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            <  1:1969596/93‑1 (MQ=11)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            <  2:3789592/93‑1 (MQ=11)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            <  2:5543956/93‑1 (MQ=11)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            >  1:4195661/1‑93 (MQ=11)
                                                                               ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc            <  2:4452808/93‑1 (MQ=11)
                                                                                tGATTTCCTCAAGAAGGAGATCATTAAAGCCATGAAAAAACGAGCTGCGCTGGATGAAATAATACGGGGGTTGCTGAGTGAATATATCGAACa           <  1:5325720/93‑1 (MQ=9)
                                                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa           <  2:2987389/93‑1 (MQ=14)
                                                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa           <  1:4595290/93‑1 (MQ=14)
                                                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa           >  2:5375347/1‑93 (MQ=14)
                                                                                tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGAAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGCGAATATGTCGAACa           >  1:387439/1‑93 (MQ=2)
                                                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg          >  1:1009252/1‑93 (MQ=39)
                                                                                 gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg          <  1:3222314/93‑1 (MQ=39)
                                                                                 gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg          >  1:4249001/1‑93 (MQ=14)
                                                                                 gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg          >  1:5132580/1‑93 (MQ=14)
                                                                                 gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg          <  1:1833045/93‑1 (MQ=14)
                                                                                 gATTTCCCGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCAGGACGAACTGATACCGGGGTTGCTGAGTGAATATAGCGAACAg          >  2:5471772/1‑93 (MQ=11)
                                                                                  aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGt         <  1:1053421/93‑1 (MQ=40)
                                                                                  aTTTCCTGAAACGGGATATCATAAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGt         <  1:4873411/93‑1 (MQ=255)
                                                                                  aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt         <  2:159048/93‑1 (MQ=14)
                                                                                  aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt         <  1:1746160/93‑1 (MQ=14)
                                                                                   tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc        >  1:3135364/1‑93 (MQ=255)
                                                                                    ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa       <  2:5327633/93‑1 (MQ=255)
                                                                                     tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg      >  2:2719441/1‑93 (MQ=255)
                                                                                     tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg      <  2:4238870/93‑1 (MQ=255)
                                                                                     tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg      >  1:1518965/1‑93 (MQ=255)
                                                                                      ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg     <  2:2738822/93‑1 (MQ=255)
                                                                                        tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt   >  1:3396879/1‑93 (MQ=255)
                                                                                        agAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt   <  1:1417995/92‑1 (MQ=255)
                                                                                         gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa  <  2:3626604/93‑1 (MQ=255)
                                                                                         |                                                                                            
GCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTA  >  NC_000913/580420‑580601

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