Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,125,741 G→T 7.2% Q60H (CAG→CAT) ‡ 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,7410GT 7.2% 279.1 / 8.1 80Q60H (CAG→CAT) ‡rimJribosomal‑protein‑S5‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (3/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGTTATT  >  NC_000913/1125593‑1125889
                                                                                                                                                    |                                                                                                                                                    
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:69924/1‑149 (MQ=255)
                                     gCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAg                                                                                                                 <  4:681286/149‑1 (MQ=255)
                                     gCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAg                                                                                                                 <  4:81600/149‑1 (MQ=255)
                                      cATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAgg                                                                                                                <  3:172644/149‑1 (MQ=255)
                                        tGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGt                                                                                                               >  1:36297/1‑148 (MQ=255)
                                        tGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCCCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGt                                                                                                               >  1:423695/1‑148 (MQ=255)
                                           tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCt                                                                                                                                              >  3:235404/1‑114 (MQ=255)
                                           tCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCt                                                                                                                                              <  4:235404/114‑1 (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                                                                                                              >  2:342310/1‑110 (MQ=255)
                                                                               gaATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGACGATTAACGAATTTCATAAACAAGGtt                                                                                                              <  1:342310/110‑1 (MQ=255)
                                                                                  tCGCCATTTACTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAgaga                                                                    <  3:608319/149‑1 (MQ=255)
                                                                                          tCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                             <  1:274118/148‑1 (MQ=255)
                                                                                          tCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                             <  2:431567/148‑1 (MQ=255)
                                                                                             tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGattatt                                                               >  3:891904/1‑143 (MQ=255)
                                                                                             tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGattatt                                                               <  4:891904/143‑1 (MQ=255)
                                                                                             tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                         >  4:232988/1‑149 (MQ=255)
                                                                                             tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                         >  3:188541/1‑149 (MQ=255)
                                                                                             tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGt                                                          >  2:300894/1‑148 (MQ=255)
                                                                                              cAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                             >  1:813169/1‑144 (MQ=255)
                                                                                              cAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTgg                                                             <  2:813169/144‑1 (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                                                          <  4:647817/116‑1 (MQ=255)
                                                                                                                             cACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGt                                                          >  3:647817/1‑116 (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)
                                                                                                                                                    ggccaggcTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGTTAtt  <  3:797172/149‑1 (MQ=255)
                                                                                                                                                    ggccaggcTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGTTAt   <  2:36297/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATCTCGGTTATT  >  NC_000913/1125593‑1125889

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: