Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,146,744 C→T 6.5% noncoding (156/169 nt) psrD → novel Novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,7440CT 6.5% 181.6 / 2.4 65noncoding (156/169 nt)psrDnovel Novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (37/22);  new base T (2/2);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGTTGAG  >  NC_000913/1146598‑1146889
                                                                                                                                                  |                                                                                                                                                 
ttttttCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTTCACGAAAATAAAGCAGTAACACGCAATGTGTGCGAAGTATTGGCAAAAGGCAACCACAGGcgg                                                                                                                                                 >  2:378530/1‑147 (MQ=255)
      ccATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCtttt                                                                                                                                           <  1:628294/149‑1 (MQ=255)
        aTCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAGCCTGTGGTTGCCttttgc                                                                                                                                         <  1:944687/149‑3 (MQ=255)
        aTCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                                                                                                                         >  2:995238/1‑149 (MQ=255)
          cAGATATAGCGTATTGAAGATAGCGATTTTAAACTATGCGCTGCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCGGCCTTTTTCtt                                                                                                                                       <  1:4050/149‑1 (MQ=255)
           aGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                                                                                                                      <  1:24356/149‑1 (MQ=255)
            gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                                     >  2:244614/1‑149 (MQ=255)
            gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                                                                                                                     >  1:2865/1‑149 (MQ=255)
               ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                                                                                                                                  <  2:1102915/149‑1 (MQ=255)
               ataGCGTATTGATGATAGCCATTTTAAACTATGCGATTCGTTTTGCAGGTTGCTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGATTTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                                                                                                                                  <  1:349062/149‑1 (MQ=255)
                   cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                                                                                                                              <  2:915990/149‑1 (MQ=255)
                   cGTATTGATGATAGCCATTTTAAACTAAGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                                                                                                                              <  2:1100055/149‑1 (MQ=255)
                          atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                                                                                                                       <  1:880334/149‑1 (MQ=255)
                             atAGCCATTTTAAACTATGCGCGTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTATCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTGCTTTGACTCTATGACGGTACaa                                                                                                                    >  2:901814/1‑149 (MQ=255)
                                 ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAgat                                                                                                                >  1:116047/1‑147 (MQ=255)
                                 ccATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt                                                                                                                >  2:143529/1‑149 (MQ=255)
                                  cATTTTAAACAATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa                                                                                                               >  2:226770/1‑149 (MQ=255)
                                            tctGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAACTAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgcgc                                                                                                     <  2:512812/147‑1 (MQ=255)
                                                cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTa                                                                                                 >  1:341261/1‑149 (MQ=255)
                                                 ccTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAATCAGTAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGTCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTAt                                                                                                <  1:114239/148‑1 (MQ=255)
                                                   ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTATTCTTAGACTCTATGACGTTACAAAGTTAATAAGCGCGCCCTATGc                                                                                              >  1:78199/1‑149 (MQ=255)
                                                               gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGGTGCCTTTTTCTTTGCCTCTATGACGTTACAAAGTTAAAATGCTCGCCCTATGAAAAAGGTAAAAt                                                                                  >  1:815124/1‑149 (MQ=255)
                                                               gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCAGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTGTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAAt                                                                                  >  1:296079/1‑149 (MQ=255)
                                                                  gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAc                                                                               <  2:845534/149‑1 (MQ=255)
                                                                   aTGTTTGTTATCAGCACTGAACGAAAATACAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTAcc                                                                              <  1:244614/149‑1 (MQ=255)
                                                                     gtttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAAGTTACAAAGTTAATATGAGCGCCCTATGCAAAAGGTAAAATTAcccc                                                                            >  2:474851/1‑149 (MQ=255)
                                                                      tttgttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCATAAGGTAAAATTACCCCt                                                                           >  1:578024/1‑149 (MQ=255)
                                                                         gttATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGCCTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGGAAATTTACCCCTGaa                                                                        >  1:958472/1‑148 (MQ=255)
                                                                            aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                                     <  1:55057/149‑1 (MQ=255)
                                                                            aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                                     <  2:135777/149‑1 (MQ=255)
                                                                            aTCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                                     >  1:490293/1‑149 (MQ=255)
                                                                            aTCAGCACTGAACCAAAATACAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGTCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTGACAAAGTTAATATGTGCGCCCTATGCAAAAGGTAAAATTACCCCTGActct                                                                     <  2:1063580/149‑1 (MQ=255)
                                                                               aGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGa                                                                  <  1:143529/149‑1 (MQ=255)
                                                                                  aCTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATcc                                                               >  1:1064116/1‑149 (MQ=255)
                                                                                    tGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                             <  1:19773/149‑1 (MQ=255)
                                                                                     gAACGAAAAGAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAATGGCAACCACAGGCTGCATTTTTCTTTGACTCTATGACGTTACAAAGATATTATGCGCGCGCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                            >  2:877080/1‑149 (MQ=255)
                                                                                      aaCGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGtt                                                           >  1:395135/1‑149 (MQ=255)
                                                                                      aaCGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGtt                                                           >  1:639297/1‑149 (MQ=255)
                                                                                        cGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCg                                                         >  2:290640/1‑149 (MQ=255)
                                                                                          aaaaTAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCTGACTCTCGATCCGGTTCGTAc                                                      >  1:418245/1‑149 (MQ=255)
                                                                                              tAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCATATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCCGATCGTACgga                                                   >  2:1116823/1‑148 (MQ=255)
                                                                                                aaGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCGTTGACTCTATGCCGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCGGGTTCGTACGGCTc                                                 >  2:727165/1‑149 (MQ=255)
                                                                                                     gTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGACTTTTTCTTTGACTCTATGACTTTACAAAGTTAATATGCTCGCCCTATGCAAAAGGTAAAATTACCACTGACTCTCGATACGGTTCGTACAGCTCaaaaa                                            >  1:815871/1‑149 (MQ=255)
                                                                                                      tAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTAAATATGCGCGCCCTATGCAAAACGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                                           <  2:116047/149‑1 (MQ=255)
                                                                                                         cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTATATGACGTTACAAAGTTAAAATGCGCGCCCTATTCAAAAGGAAAAATTACCACTCACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                                        >  2:774490/1‑149 (MQ=255)
                                                                                                          atgCATTGTGTTCTAATTATTGGCAAAAGGCAAGCACAGGGTTTCTTTTTCTTTGACTCTATGACGTTACTAAGTTAATATTCGCGACATATGCATAAGGTAATATTAGGCGTGACTCTCGATGAGGTTCGTTGGGCTCAAAAACGcat                                       <  2:465540/147‑3 (MQ=255)
                                                                                                                tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTa                                 <  1:995238/149‑1 (MQ=255)
                                                                                                                    tgCGAATTATTGTCTAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAg                             <  2:296079/149‑1 (MQ=255)
                                                                                                                     gCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAgg                            <  1:304904/149‑1 (MQ=255)
                                                                                                                      cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGCCTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAggg                           >  1:594643/1‑149 (MQ=255)
                                                                                                                      cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCGTGTTTACCAggg                           >  1:768193/1‑149 (MQ=255)
                                                                                                                         attattGGCAAACGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTaa                        >  2:978206/1‑148 (MQ=255)
                                                                                                                              aagaaaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAATTTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata                   >  1:174692/5‑149 (MQ=255)
                                                                                                                              aagaaaaaGGCAGCCTGTGGTTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata                   >  2:497923/5‑149 (MQ=255)
                                                                                                                                gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc                 >  1:580761/1‑149 (MQ=255)
                                                                                                                                gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc                 >  2:434473/1‑149 (MQ=255)
                                                                                                                                gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATAACCAGGGTATCTATAcc                 >  1:543321/1‑149 (MQ=255)
                                                                                                                                gCAAAAGGCAACCACAGGCTGCCGTTTTCATTGGCTCTATGACGTTACACAGTTAATATGCGCCCCCTATGCAAAAGGTAAAATTACCCGTGACTCACGCTACGGTGCGTACGGCTCAAAAACCCGATTCTTACCAGGGTATCTATAcc                 >  2:1060957/1‑149 (MQ=255)
                                                                                                                                 cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGCTCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAccc                >  2:1009255/1‑149 (MQ=255)
                                                                                                                                 aaaaaGGCAAACACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAATTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAccc                <  2:641378/148‑1 (MQ=255)
                                                                                                                                  aaaaGGCAACCACAGTCTTCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCGTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCTAAAACGCCTTGATTACCAGGGTATCTATAcccc               <  2:78199/149‑1 (MQ=255)
                                                                                                                                  aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAACCGCCTTGATTACCAGGGTATCTATAcccc               >  2:729649/1‑149 (MQ=255)
                                                                                                                                    aaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGAATCTCGATCCGGTTAGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATCCCCCTg             >  1:541635/1‑149 (MQ=255)
                                                                                                                                       gCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATc          >  1:537698/1‑149 (MQ=255)
                                                                                                                                         aaCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAg        >  1:18142/1‑149 (MQ=255)
                                                                                                                                               aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAATGTAAAATTACCCCTGACTCTAGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGTTGAg  >  2:56360/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTTTTCCATCAGATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCTGATCAGGTTGAG  >  NC_000913/1146598‑1146889

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