Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,567 C→A 5.6% A885S (GCA→TCA)  ygcB ← Cascade complex anti‑viral R‑loop helicase‑annealase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,5670CA 5.6% 255.0 / 4.7 72A885S (GCA→TCA) ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

AATATTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGT  >  NC_000913/2884419‑2884714
                                                                                                                                                    |                                                                                                                                                     
aaTATTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGc                                                                                                                                                       <  4:283904/149‑1 (MQ=255)
 atatTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGGTTTGCa                                                                                                                                                      <  4:467898/149‑1 (MQ=255)
  tatTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAg                                                                                                                                                     <  3:424911/149‑1 (MQ=255)
   atTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAgg                                                                                                                                                    <  2:251201/149‑1 (MQ=255)
   atTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAgg                                                                                                                                                    <  3:560325/149‑1 (MQ=255)
    ttCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAggg                                                                                                                                                   <  3:481110/149‑1 (MQ=255)
     tCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGACGATGTACATTGTGCACCTTCCCTACTTAAGTAAGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGCTTTGCAGGGa                                                                                                                                                  <  4:380302/149‑1 (MQ=255)
      cATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGAt                                                                                                                                                 <  1:300511/149‑1 (MQ=255)
       aTACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATg                                                                                                                                                >  3:279075/1‑149 (MQ=255)
             tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                          <  3:364962/149‑1 (MQ=255)
               aaTGTAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGt                                                                                                                                        >  3:551704/1‑149 (MQ=255)
                atatAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                       <  4:279075/149‑1 (MQ=255)
                 tataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCa                                                                                                                                      >  1:130641/1‑149 (MQ=255)
                     aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATccc                                                                                                                                  <  3:434319/149‑1 (MQ=255)
                            aaacgtagtttttctgaagtagtagatgaagatgggttactttggctggaagggaaacagaatctggatggatgggtctggcagggtaacagtattgttattacctatacagGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑Tcc                                                                                                                         >  4:425841‑M1/113‑149 (MQ=255)
                               ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCt                                                                                                                      <  2:102531/149‑1 (MQ=255)
                                cGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTg                                                                                                                     >  4:32050/1‑149 (MQ=255)
                                 gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGt                                                                                                                    >  2:518510/1‑149 (MQ=255)
                                       ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGt                                                                                                              <  1:546069/149‑1 (MQ=255)
                                       ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGACAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGt                                                                                                              >  3:192174/1‑149 (MQ=255)
                                            agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGACAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGtaataa                                                                                                         <  4:192174/149‑1 (MQ=255)
                                                 aTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGTGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATa                                                                                                    >  4:550448/1‑149 (MQ=255)
                                                         tgtgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTAcc                                                                                            <  2:18389/149‑1 (MQ=255)
                                                         agatgggttactttggctggaagggaaacagaatctggatggatgggtctggcagggtaacagtattgttattacctatacagGGGATTTGCAg                                                                                                                                                     <  2:232121‑M1/11‑1 (MQ=255)
                                                         agatgggttactttggctggaagggaaacagaatctggatggatgggtctggcagggtaacagtattgttattacctatacagGGGATTTGCAg                                                                                                                                                     >  1:232121‑M1/84‑94 (MQ=255)
                                                           tgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCt                                                                                          >  3:273201/1‑149 (MQ=255)
                                                             cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGc                                                                                        >  1:149804/1‑149 (MQ=255)
                                                              aCCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                       >  3:35346/1‑149 (MQ=255)
                                                               ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                       >  3:297721/1‑102 (MQ=255)
                                                               ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                       <  4:297721/102‑1 (MQ=255)
                                                                 ttCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGa                                                                                    <  4:273201/149‑1 (MQ=255)
                                                                  tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                   >  3:121545/1‑149 (MQ=255)
                                                                     cTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                >  3:65692/1‑149 (MQ=255)
                                                                         ttAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACaa                                                                                                      >  3:566216/1‑123 (MQ=255)
                                                                         ttAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACaa                                                                                                      <  4:566216/123‑1 (MQ=255)
                                                                             gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACt                                                                                                  >  1:25567/1‑123 (MQ=255)
                                                                             gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACt                                                                                                  <  2:25567/123‑1 (MQ=255)
                                                                               aGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                      >  2:521650/1‑149 (MQ=255)
                                                                               aGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                      >  3:126439/1‑149 (MQ=255)
                                                                                   aacagaatctggatggatgggtctggcagggtaacagtattgttattacctatacagGGGAT‑‑GAAGGGATGACTCTGGTCATccc                                                                                                                                  <  2:496943‑M1/28‑1 (MQ=255)
                                                                                   aacagaatctggatggatgggtctggcagggtaacagtattgttattacctatacagGGGAT‑‑GAAGGGATGACTCTGGTCATccc                                                                                                                                  >  1:496943‑M1/58‑85 (MQ=255)
                                                                                   aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCt                                                                  >  4:101303/1‑149 (MQ=255)
                                                                                   aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCt                                                                  >  3:505791/1‑149 (MQ=255)
                                                                                     aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGt                                                                                                              <  1:156442/103‑1 (MQ=255)
                                                                                     aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                <  1:161906/149‑1 (MQ=255)
                                                                                     aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                <  2:174345/149‑1 (MQ=255)
                                                                                     aaaCCGTTATCGGTCTTATTATCGACATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGt                                                                                                              >  2:156442/1‑103 (MQ=255)
                                                                                       aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCGCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                              >  2:339555/1‑149 (MQ=255)
                                                                                        ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                             <  4:505791/149‑1 (MQ=255)
                                                                                          gTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTccc                                                           <  4:223358/149‑1 (MQ=255)
                                                                                               tGGTCTTATTATCGTCATTGAGAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                      <  4:551704/149‑1 (MQ=255)
                                                                                                 gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGc                                                    >  1:423978/1‑149 (MQ=255)
                                                                                                    ttattaTCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaa                                                 <  4:65692/149‑1 (MQ=255)
                                                                                                         atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                            <  2:130641/149‑1 (MQ=255)
                                                                                                          tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAc                                           >  3:224877/1‑149 (MQ=255)
                                                                                                             tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                        >  1:266913/1‑149 (MQ=255)
                                                                                                             tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                        >  2:262893/1‑149 (MQ=255)
                                                                                                                   aTAACAATTATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                  >  3:567310/1‑149 (MQ=255)
                                                                                                                    acagtattgttattacctatacagGGGAT‑‑GAAGGGATGACCAGAGTCATCCCTGCAAATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCt                                                          >  1:259859‑M1/25‑124 (MQ=255)
                                                                                                                    acagtattgttattacctatacagGGGAT‑‑GAAGGGATGACCAGAGTCATCCCTGCAAATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCt                                                          <  2:259859‑M1/100‑1 (MQ=255)
                                                                                                                        aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                             >  1:211565/1‑149 (MQ=255)
                                                                                                                        aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                             >  2:323616/1‑149 (MQ=255)
                                                                                                                        aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                             >  1:120685/1‑149 (MQ=255)
                                                                                                                          tCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                           <  4:121545/149‑1 (MQ=255)
                                                                                                                          tCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                           >  1:65711/1‑149 (MQ=255)
                                                                                                                            aTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTactact                         >  2:259460/1‑149 (MQ=255)
                                                                                                                              tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                       <  3:496308/149‑1 (MQ=255)
                                                                                                                                ccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAg                     >  4:50864/1‑149 (MQ=255)
                                                                                                                                        aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTa             <  1:400010/149‑1 (MQ=255)
                                                                                                                                           tGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGt          <  2:423978/149‑1 (MQ=255)
                                                                                                                                               ATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTcc      >  2:266913‑M1/1‑149 (MQ=255)
                                                                                                                                                 ttGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg    <  4:155217/149‑1 (MQ=255)
                                                                                                                                                   gCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGt  <  4:35346/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATATTCATACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCA‑‑TCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGT  >  NC_000913/2884419‑2884714

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: