Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,317 A→C 18.0% 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,3170AC18.0% 442.3 / 94.3 192N21H (AAT→CAT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (69/88);  new base C (13/22);  total (82/110)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 195 total aligned reads displayed.

TAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTC  >  NC_000913/2560218‑2560416
                                                                                                   |                                                                                                   
taagccgAATAAATTATTAAAAAACAAAAACTTCAATCAATATTAGCAACCAAAGATCCATACTCCGAAAATAGGATCACTACACATAAAAGCCCCCCGa                                                                                                     <  1:3730198/94‑1 (MQ=255)
 aaGTGGAATAAATTATTCAATATCAAAAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGaa                                                                                                    <  1:1409694/100‑1 (MQ=255)
  aGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAt                                                                                                   >  1:2456332/1‑100 (MQ=255)
   gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCAtg                                                                                                  <  1:3906911/100‑1 (MQ=255)
   gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtg                                                                                                  <  2:2069352/100‑1 (MQ=255)
    tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                 <  1:232756/100‑1 (MQ=255)
       aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                              >  1:573495/1‑100 (MQ=255)
        aTAAATTATTCAATATCAAGAACTTAAATCCATTTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                             <  2:894430/100‑1 (MQ=255)
        aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                             <  1:3602104/100‑1 (MQ=255)
        aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                             >  2:692337/1‑100 (MQ=255)
        aTAAATTAATCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGAACCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                             <  1:1097640/100‑1 (MQ=255)
          aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                           >  2:2126230/1‑100 (MQ=255)
           aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                          >  1:2877187/1‑100 (MQ=255)
           aaTTATTCAATATAAAGAACTTAAATACATGTTTGCAACGAAAGATACAGAGTTCAAAAATAGGATAAAAACAAATAAAAGCGAGCGGAAtgcaaaaaca                                                                                          >  1:4026738/1‑93 (MQ=255)
            attattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                          <  1:893461/99‑1 (MQ=255)
             ttattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATg                                                                                        >  1:2162785/1‑100 (MQ=255)
             ttattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATg                                                                                        >  2:3712955/1‑100 (MQ=255)
               attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                      >  2:2183968/1‑100 (MQ=255)
                 tCAATATCAAGAACTTAAGTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgg                                                                                    <  2:374555/100‑1 (MQ=38)
                  cAATATCAAGAGCTTATATCCTTTTTTGCAACGATAGATCCAGAGTTCGAAAATAGGATCAACACAAATATAAGCCGCCGGCATGTTGCAGCATTCCggg                                                                                   <  1:187340/100‑1 (MQ=255)
                  cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                   <  2:3288197/100‑1 (MQ=38)
                   atcttCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                  <  2:4110383/96‑1 (MQ=38)
                   aaTATCAAAAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                  <  2:3013288/100‑1 (MQ=38)
                    atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                 <  2:3265192/100‑1 (MQ=38)
                    atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                 <  2:2738834/100‑1 (MQ=38)
                    atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                 >  1:3995884/1‑100 (MQ=255)
                     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                <  1:1314270/100‑1 (MQ=38)
                     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                <  2:1186149/100‑1 (MQ=38)
                     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                <  2:4305471/100‑1 (MQ=38)
                     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                <  2:1824027/100‑1 (MQ=38)
                      ctCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                               <  2:444231/99‑1 (MQ=255)
                      atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGct                                                                               <  2:750113/100‑2 (MQ=38)
                      atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                <  1:4265533/99‑1 (MQ=38)
                      atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                               >  1:2825709/1‑100 (MQ=255)
                       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                              >  1:2490302/1‑100 (MQ=255)
                       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                              <  2:4158322/100‑1 (MQ=255)
                       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                              <  1:3882836/100‑1 (MQ=255)
                       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGCTAGAAAATAGGATCAACACAAGTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                              <  1:1052170/100‑1 (MQ=255)
                        cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                             >  2:3384288/1‑100 (MQ=255)
                         aaGATCTTAAATCTATGTTTGCAACGAAAGATCCAGACTTCGAATATAGGATCAACCCAAATAAATGTCCCCGGAATGTAGCAACATGCTGGGGTCgg                                                                              <  1:3777194/98‑1 (MQ=255)
                           gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                          <  2:830915/100‑1 (MQ=255)
                            aaCTTAAATCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAATAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                          <  1:3717810/99‑1 (MQ=255)
                            aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                         <  1:933251/100‑1 (MQ=255)
                             aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                        <  2:2825725/100‑1 (MQ=255)
                             aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                         <  1:1842651/99‑1 (MQ=255)
                              cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                       <  2:3139065/100‑1 (MQ=255)
                                tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                     >  2:1949426/1‑100 (MQ=255)
                                 aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                    <  1:159219/100‑1 (MQ=255)
                                 aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                    <  2:2373382/100‑1 (MQ=255)
                                   aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                  <  2:545150/100‑1 (MQ=255)
                                   aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                  <  2:2710819/100‑1 (MQ=255)
                                   aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                  >  2:3836363/1‑100 (MQ=255)
                                    tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                 <  2:3796917/100‑1 (MQ=255)
                                    tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                 <  1:4054315/100‑1 (MQ=255)
                                     ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                <  1:4333270/100‑1 (MQ=255)
                                     ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                <  1:2699958/100‑1 (MQ=255)
                                      cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                               >  1:3347186/1‑100 (MQ=255)
                                         gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                            <  2:2936457/100‑1 (MQ=255)
                                             gCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATTCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                        <  2:4179326/100‑1 (MQ=255)
                                             gCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                        <  1:752098/100‑1 (MQ=255)
                                                aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                     >  2:3873663/1‑100 (MQ=255)
                                                  gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                   >  2:358108/1‑100 (MQ=255)
                                                  gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                   >  2:249632/1‑100 (MQ=255)
                                                   aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                  >  2:3299076/1‑100 (MQ=255)
                                                    aaaaTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                 <  1:3068107/97‑1 (MQ=255)
                                                    aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAGGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                 <  1:873371/100‑1 (MQ=255)
                                                     aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                <  1:4260048/100‑1 (MQ=255)
                                                      gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                               >  1:1659539/1‑100 (MQ=255)
                                                      gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                               >  2:1106174/1‑100 (MQ=255)
                                                       aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                              <  2:1028429/100‑1 (MQ=255)
                                                       aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                              >  2:2650047/1‑100 (MQ=255)
                                                            gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                         <  1:3014555/100‑1 (MQ=255)
                                                             agTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                        <  2:2069253/100‑1 (MQ=255)
                                                              gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                       >  1:1243830/1‑100 (MQ=255)
                                                              gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                       >  1:3490651/1‑100 (MQ=255)
                                                                 cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACCAGCTAAAGGAGAATTTCTCATGTCTGATATg                                    <  1:4102585/100‑1 (MQ=255)
                                                                 cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                    >  1:2385482/1‑100 (MQ=255)
                                                                 cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGACCAACAAGCTAAAGGAGACTTTCTCATGTCTGCTATg                                    >  1:4137623/1‑100 (MQ=255)
                                                                     aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                <  1:2944729/100‑1 (MQ=255)
                                                                          gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                           <  2:957005/100‑1 (MQ=255)
                                                                          gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                           <  2:755703/100‑1 (MQ=255)
                                                                           aTCAACACAAATAAAAGCCCCCTGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACCAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                            <  1:2511650/98‑1 (MQ=255)
                                                                             cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                          <  1:448555/98‑1 (MQ=255)
                                                                               acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                      <  2:2407139/100‑1 (MQ=255)
                                                                                cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                     >  1:49950/1‑100 (MQ=255)
                                                                                cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                     >  2:2571759/1‑100 (MQ=255)
                                                                                 acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                    >  1:15386/1‑100 (MQ=255)
                                                                                        aaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa             <  2:1395941/100‑1 (MQ=255)
                                                                                          cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAGTTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca           >  1:2816525/2‑100 (MQ=37)
                                                                                          aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca           >  2:1880962/1‑100 (MQ=255)
                                                                                            cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat         >  2:149418/1‑100 (MQ=38)
                                                                                             ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg        >  1:2100066/1‑100 (MQ=38)
                                                                                             ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg        >  2:3446544/1‑100 (MQ=38)
                                                                                             ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATCTCTGACATGCTAGCAGTAGAACAAGAAACAAACAatg        >  1:2960068/1‑100 (MQ=38)
                                                                                             ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg        >  2:1744347/1‑100 (MQ=255)
                                                                                             acccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGAAGAATTTCTAATGTATGATATGCTAACAGTAGAACAAGAAacaaaaaa          >  1:2659441/2‑98 (MQ=25)
                                                                                              cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga       >  1:2061351/1‑100 (MQ=38)
                                                                                               ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat      >  2:2898815/1‑100 (MQ=255)
                                                                                               ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga       >  2:2527150/1‑99 (MQ=255)
                                                                                                   aaTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc  <  2:1637103/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTC  >  NC_000913/2560218‑2560416

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