Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,125,740 A→C 6.2% Q60P (CAG→CCG) ‡ rimJ → ribosomal‑protein‑S5‑alanine N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,125,7400AC 6.2% 282.2 / 9.4 80Q60P (CAG→CCG) ‡rimJribosomal‑protein‑S5‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base C (2/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGT  >  NC_000913/1125592‑1125885
                                                                                                                                                    |                                                                                                                                                 
gTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTACTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGca                                                                                                                                                   >  4:608319/1‑149 (MQ=255)
 tGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGcag                                                                                                                                                  <  3:428556/149‑1 (MQ=255)
         aCCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCt                                                                                                                                          >  4:741342/1‑149 (MQ=255)
           cacaGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGCCTGGCCTGCCAGCCTga                                                                                                                                        <  4:381723/149‑2 (MQ=255)
            acaGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGCCTGGCCTGCCAGCCTgat                                                                                                                                       <  3:813821/149‑3 (MQ=255)
             caGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTggg                                                                                                                                       >  1:165194/1‑148 (MQ=255)
              aGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGga                                                                                                                                     >  3:196127/1‑149 (MQ=255)
               gACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGga                                                                                                                                     >  2:274118/1‑148 (MQ=255)
                aCCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGgatg                                                                                                                                   >  3:660285/1‑149 (MQ=255)
                 ccGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGgatg                                                                                                                                   <  2:777598/148‑1 (MQ=255)
                      tGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAAc                                                                                                                             <  3:156788/149‑1 (MQ=255)
                              cGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTc                                                                                                                      <  2:845135/148‑1 (MQ=255)
                                tCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATa                                                                                                                   >  3:471213/1‑149 (MQ=255)
                                  tGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATaaa                                                                                                                 <  4:789170/149‑1 (MQ=255)
                                      gCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAg                                                                                                             <  4:681286/149‑1 (MQ=255)
                                      gCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAg                                                                                                             <  4:81600/149‑1 (MQ=255)
                                      gCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAg                                                                                                             >  4:69924/1‑149 (MQ=255)
                                       cATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAgg                                                                                                            <  3:172644/149‑1 (MQ=255)
                                         tGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGt                                                                                                           >  1:36297/1‑148 (MQ=255)
                                         tGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCCCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGt                                                                                                           >  1:423695/1‑148 (MQ=255)
                                            tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCt                                                                                                                                          <  4:235404/114‑1 (MQ=255)
                                            tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCt                                                                                                                                          >  3:235404/1‑114 (MQ=255)
                                            tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAg                                                                                                       >  3:488803/1‑149 (MQ=255)
                                            tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCAGTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCa                                                                                                        <  2:380404/148‑1 (MQ=255)
                                              gTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCGTTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGc                                                                                                      <  2:538500/148‑1 (MQ=255)
                                                     cTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACt                                                                                              >  4:797172/1‑149 (MQ=255)
                                                      tGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACtt                                                                                             <  3:69924/149‑1 (MQ=255)
                                                       ggCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACttt                                                                                            <  4:660285/149‑1 (MQ=255)
                                                        gCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACttt                                                                                            <  1:832301/148‑1 (MQ=255)
                                                          gTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGc                                                                                         <  4:330658/149‑1 (MQ=255)
                                                              tGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTa                                                                                      >  1:316956/1‑148 (MQ=255)
                                                               gCGGATTATTACGCAGAGAGTCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTAt                                                                                     <  2:803056/148‑1 (MQ=255)
                                                                cGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTc                                                                                   <  4:241614/149‑1 (MQ=255)
                                                                cGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTc                                                                                   >  4:632290/1‑149 (MQ=255)
                                                                    ttattaCGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGAc                                                                                >  2:579718/1‑148 (MQ=255)
                                                                     tattaCGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGAcc                                                                               >  1:27505/1‑148 (MQ=255)
                                                                      attaCCCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGAccc                                                                              >  2:446483/1‑148 (MQ=255)
                                                                       ttaCGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCgg                                                                            >  4:57786/1‑149 (MQ=255)
                                                                            cAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGa                                                                        >  2:530620/1‑148 (MQ=255)
                                                                               agaATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGaaaaa                                                                    <  3:741342/149‑1 (MQ=255)
                                                                                gaATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGACGATTAACGAATTTCATAAACAAGGtt                                                                                                          <  1:342310/110‑1 (MQ=255)
                                                                                gaATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGACGATTAACGAATTTCATAAACAAGGtt                                                                                                          >  2:342310/1‑110 (MQ=255)
                                                                                   tCGCCATTTACTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAgaga                                                                <  3:608319/149‑1 (MQ=255)
                                                                                           tCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                         <  1:274118/148‑1 (MQ=255)
                                                                                           tCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                         <  2:431567/148‑1 (MQ=255)
                                                                                              tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGattatt                                                           <  4:891904/143‑1 (MQ=255)
                                                                                              tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGattatt                                                           >  3:891904/1‑143 (MQ=255)
                                                                                              tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                     >  4:232988/1‑149 (MQ=255)
                                                                                              tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                     >  3:188541/1‑149 (MQ=255)
                                                                                              tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGt                                                      >  2:300894/1‑148 (MQ=255)
                                                                                               cAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                         <  2:813169/144‑1 (MQ=255)
                                                                                               cAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                         >  1:813169/1‑144 (MQ=255)
                                                                                               cAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTg                                                    >  4:712783/1‑149 (MQ=255)
                                                                                                aaGCCCTGGGAGCCAGTGCACGACGAAAGCCACTGTTATCCATCAGGCTGGCATGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGc                                                   >  4:624026/1‑149 (MQ=255)
                                                                                                                     gacgaAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTcc                                         >  3:851390/1‑138 (MQ=255)
                                                                                                                     gacgaAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTcc                                         <  4:851390/138‑1 (MQ=255)
                                                                                                                      acgaAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGt                             >  4:714740/1‑149 (MQ=255)
                                                                                                                       cgaAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTg                            <  3:632290/149‑1 (MQ=255)
                                                                                                                         aaaGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAAt                                      >  2:381155/1‑137 (MQ=255)
                                                                                                                         aaaGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAAt                                      <  1:381155/137‑1 (MQ=255)
                                                                                                                          aaGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGct                         >  4:756767/1‑149 (MQ=255)
                                                                                                                              cACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGCTTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGATAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCttt                     >  3:797165/1‑149 (MQ=255)
                                                                                                                              cACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGt                                                      >  3:647817/1‑116 (MQ=255)
                                                                                                                              cACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGt                                                      <  4:647817/116‑1 (MQ=255)
                                                                                                                               aCTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAgaga                                                                <  4:97980/105‑1 (MQ=255)
                                                                                                                               aCTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAgaga                                                                >  3:97980/1‑105 (MQ=255)
                                                                                                                               aCTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCttt                     >  1:204484/1‑148 (MQ=255)
                                                                                                                                  gTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCa                  <  1:579718/148‑1 (MQ=255)
                                                                                                                                  gTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCAt                 <  3:687355/149‑1 (MQ=255)
                                                                                                                                    tcatccccAGCCTGGCCTGCCAGCCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCtttt                    >  4:888612/8‑144 (MQ=255)
                                                                                                                                    tcatccccAGCCTGGCCTGCCAGCCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCtttt                    <  3:888612/137‑1 (MQ=255)
                                                                                                                                    tATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGaaaaa                                                                    >  3:512374/1‑96 (MQ=255)
                                                                                                                                    tATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGaaaaa                                                                    <  4:512374/96‑1 (MQ=255)
                                                                                                                                    tATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATg                <  2:325959/148‑1 (MQ=255)
                                                                                                                                         aTCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGc           <  1:42377/148‑1 (MQ=255)
                                                                                                                                          tCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGATTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTg            <  2:727092/146‑1 (MQ=255)
                                                                                                                                          tCAGGCTGGCAGGCCAGGCAGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTg            >  1:727092/1‑146 (MQ=255)
                                                                                                                                          ccAGCCTGGCCTGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTa         >  3:332533/2‑149 (MQ=255)
                                                                                                                                            aagCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTAtc       >  4:650449/3‑149 (MQ=255)
                                                                                                                                               cTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCg    <  4:188541/149‑1 (MQ=255)
                                                                                                                                                tGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCgg   <  4:592855/149‑1 (MQ=255)
                                                                                                                                                 ggCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGt  >  4:329125/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGT  >  NC_000913/1125592‑1125885

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