Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,317 A→C 15.3% N21H (AAT→CAT)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3170AC15.3% 415.3 / 69.8 157N21H (AAT→CAT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (58/75);  new base C (10/14);  total (68/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 158 total aligned reads displayed.

AAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGT  >  NC_000913/2560219‑2560414
                                                                                                  |                                                                                                 
  gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtg                                                                                                >  1:1389407/1‑100 (MQ=255)
   tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACCCAAATAAAAGCCCCCGGAAtgc                                                                                               <  2:2289956/100‑1 (MQ=255)
    ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                              <  2:1590953/100‑1 (MQ=255)
      aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                            <  2:578625/100‑1 (MQ=255)
      aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                            <  1:2593678/100‑1 (MQ=255)
       aTAAATTATTCAATATCCAGAACTTAAATCCACGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                           <  2:3729505/100‑1 (MQ=255)
       aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                           <  1:2993640/100‑1 (MQ=255)
         aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                         <  2:1101850/100‑1 (MQ=255)
          aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                        <  2:3269946/100‑1 (MQ=255)
          aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                        <  2:3523055/100‑1 (MQ=255)
                 cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                 <  2:1917487/100‑1 (MQ=38)
                 cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                 <  1:33776/100‑1 (MQ=38)
                 cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                 >  1:3670762/1‑100 (MQ=255)
                 cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                 >  1:2416466/1‑100 (MQ=255)
                  aaTATCAAGAACTTAAATCCATGTTTGCAACGATAGATCAAGCGTTGGAAAATAAGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                 <  2:3861637/99‑1 (MQ=25)
                  aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCGGGAATGCTGCAACATGCCgggg                                                                                <  1:3605615/100‑1 (MQ=255)
                  aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                <  1:2639866/100‑1 (MQ=38)
                  aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                <  1:1860143/100‑1 (MQ=255)
                   atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                               <  1:378612/100‑1 (MQ=38)
                   atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                               <  1:983588/100‑1 (MQ=38)
                   atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                               <  1:1024969/100‑1 (MQ=38)
                   atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                >  2:2001917/1‑99 (MQ=255)
                    tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACGAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                              <  1:2410057/100‑1 (MQ=255)
                    tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                              <  1:1369428/100‑1 (MQ=255)
                      tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                            <  1:716053/100‑1 (MQ=255)
                      tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                            >  2:2594922/1‑100 (MQ=255)
                      tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                            <  1:984375/100‑1 (MQ=255)
                      tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                            <  2:3140455/100‑1 (MQ=255)
                       cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACCCAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                           <  1:41127/100‑1 (MQ=255)
                       cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttt                                                                           <  1:2311321/100‑4 (MQ=37)
                        aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                          <  1:1227198/100‑1 (MQ=255)
                        aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                          <  1:3731914/100‑1 (MQ=255)
                        aaGAACTTAAATCCAGGTTTGCAACGAAAGATCCAGAGGTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                               <  2:1188054/95‑1 (MQ=37)
                         aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                         <  1:1666129/100‑1 (MQ=255)
                          gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                        <  2:3161862/100‑1 (MQ=255)
                            aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                      <  1:3275234/100‑1 (MQ=255)
                            aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                      <  2:499744/100‑1 (MQ=255)
                             cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                     >  1:2491227/1‑100 (MQ=255)
                              ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                    >  1:3366512/1‑100 (MQ=255)
                              ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                    <  1:1894014/100‑1 (MQ=255)
                                aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                  <  1:1575340/100‑1 (MQ=255)
                                 aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                 >  2:402942/1‑100 (MQ=255)
                                    ccATGTTTGCAACGAAAGATCCGGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                              <  2:1840109/100‑1 (MQ=255)
                                     cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                <  2:3015755/97‑1 (MQ=255)
                                     cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                             >  2:3061641/1‑100 (MQ=255)
                                       tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCCGAATGCTGCATCATGCCGGGGGCGGTATCAGAAACAAGCtca                                                           <  1:2668034/100‑3 (MQ=255)
                                        gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                          <  1:2081198/100‑1 (MQ=255)
                                             cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                     <  1:1075270/100‑1 (MQ=255)
                                              aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                    >  1:1752548/1‑100 (MQ=255)
                                              aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                    <  1:1987665/100‑1 (MQ=255)
                                               aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCACCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                   >  1:1563092/1‑100 (MQ=255)
                                                cgatctaTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAtt                                                  <  1:3689134/94‑1 (MQ=255)
                                                ccaaaGATGGAGATTTCGAAAATATGATCAACACAAATAAAAGCCACCGGAATGCTGCAACAGGCCGGGGGCGGGATGAGACACAAGCTAAAGGAGAAtt                                                  <  2:2346255/98‑1 (MQ=255)
                                                   aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                               >  1:2329684/1‑100 (MQ=255)
                                                      aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                            <  2:3586621/100‑1 (MQ=255)
                                                      aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                            >  1:2025611/1‑100 (MQ=255)
                                                        ccAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGt                                          >  2:1293201/1‑100 (MQ=255)
                                                         cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                         >  1:3645464/1‑100 (MQ=255)
                                                             tctCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                     <  2:2209713/98‑1 (MQ=255)
                                                               tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                   >  2:3104232/1‑100 (MQ=255)
                                                                cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGataga                                  >  1:2209725/1‑98 (MQ=255)
                                                                cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                  >  1:293626/1‑100 (MQ=255)
                                                                 gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                 >  1:451194/1‑100 (MQ=255)
                                                                 gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                 >  2:3753222/1‑100 (MQ=255)
                                                                 gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                 >  2:1692997/1‑100 (MQ=255)
                                                                  aaaaTAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                <  1:697323/100‑1 (MQ=38)
                                                                  aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                <  1:125308/100‑1 (MQ=255)
                                                                  aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                <  2:3646205/100‑1 (MQ=255)
                                                                  aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                <  2:3645436/100‑1 (MQ=255)
                                                                  aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCGCAAGGAGAATTTCTCATGTCTGATATGCt                                >  2:1764869/1‑100 (MQ=255)
                                                                     aTAGCATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                             <  1:581157/100‑1 (MQ=255)
                                                                      tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                            <  2:1486679/100‑1 (MQ=255)
                                                                        ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                          <  1:400537/100‑1 (MQ=255)
                                                                          aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                        >  1:166953/1‑100 (MQ=255)
                                                                            cAACACAAATAAAAGCCCCCGGAATGCTGCAGCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                      <  2:3366531/100‑1 (MQ=255)
                                                                             aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                     >  2:84690/1‑100 (MQ=255)
                                                                             aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                     <  2:3875954/100‑1 (MQ=255)
                                                                             aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                     >  1:3496708/1‑100 (MQ=255)
                                                                              acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                    >  1:3513522/1‑100 (MQ=255)
                                                                              acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                    <  2:3473112/100‑1 (MQ=255)
                                                                               cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                   >  2:169276/1‑100 (MQ=255)
                                                                               cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                   >  1:2160237/1‑100 (MQ=255)
                                                                               cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                   >  1:2636632/1‑100 (MQ=255)
                                                                                acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                  >  2:1141540/1‑100 (MQ=255)
                                                                                acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                  >  2:1755998/1‑100 (MQ=255)
                                                                                acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                  <  1:1457959/100‑1 (MQ=255)
                                                                                   aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa               >  1:774049/1‑100 (MQ=255)
                                                                                    aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac              <  2:2179466/100‑1 (MQ=255)
                                                                                    aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac              <  2:3909579/100‑1 (MQ=255)
                                                                                     tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca             >  1:391231/1‑100 (MQ=255)
                                                                                        aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTGATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac          >  1:2769972/1‑100 (MQ=255)
                                                                                        aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac          <  2:1124405/100‑1 (MQ=255)
                                                                                         cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca         >  2:611947/2‑100 (MQ=38)
                                                                                          gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca         >  2:1360005/1‑99 (MQ=38)
                                                                                          gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa        >  1:848557/1‑100 (MQ=38)
                                                                                          gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa        >  1:2524621/1‑100 (MQ=38)
                                                                                          gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa        >  2:3334843/1‑100 (MQ=38)
                                                                                             cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga     <  2:2984592/100‑1 (MQ=255)
                                                                                              ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat    >  1:3742398/1‑100 (MQ=255)
                                                                                                ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt  <  1:3574586/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
AAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGT  >  NC_000913/2560219‑2560414

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