Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,125,734 G→C 6.2% G58A (GGC→GCC)  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,7340GC 6.2% 283.7 / 8.3 81G58A (GGC→GCC) rimJribosomal‑protein‑S5‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base C (2/3);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

AAAAGTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATC  >  NC_000913/1125588‑1125880
                                                                                                                                                  |                                                                                                                                                  
aaaaGTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGc                                                                                                                                                   <  2:341111/148‑1 (MQ=255)
    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: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                                                      <  4:891904/143‑1 (MQ=255)
                                                                                                  tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGattatt                                                      >  3:891904/1‑143 (MQ=255)
                                                                                                  tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                >  3:188541/1‑149 (MQ=255)
                                                                                                  tCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGtt                                                >  4:232988/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)
                                                                                                      gCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGcct                                                                                                                                             <  4:414175/52‑3 (MQ=255)
                                                                                                      gCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGcct                                                                                                                                             >  3:414175/1‑50 (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                                 <  1:381155/137‑1 (MQ=255)
                                                                                                                             aaaGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAAt                                 >  2:381155/1‑137 (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                                                           >  3:97980/1‑105 (MQ=255)
                                                                                                                                   aCTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAgaga                                                           <  4:97980/105‑1 (MQ=255)
                                                                                                                                   aCTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCttt                >  1:204484/1‑148 (MQ=255)
                                                                                                                                      gTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCa             <  1:579718/148‑1 (MQ=255)
                                                                                                                                      gTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCAt            <  3:687355/149‑1 (MQ=255)
                                                                                                                                        tcatccccAGCCTGGCCTGCCAGCCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCtttt               <  3:888612/137‑1 (MQ=255)
                                                                                                                                        tcatccccAGCCTGGCCTGCCAGCCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCtttt               >  4:888612/8‑144 (MQ=255)
                                                                                                                                        tATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGaaaaa                                                               <  4:512374/96‑1 (MQ=255)
                                                                                                                                        tATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGaaaaa                                                               >  3:512374/1‑96 (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)
                                                                                                                                                  |                                                                                                                                                  
AAAAGTGCGCTTAACCACAGACCGACTGGTCGTGCGTCTGGTGCATGATCGTGATGCCTGGCGTCTTGCGGATTATTACGCAGAGAATCGCCATTTCCTCAAGCCCTGGGAGCCAGTGCGCGACGAAAGCCACTGTTATCCATCAGGCTGGCAGGCCAGGCTGGGGATGATTAACGAATTTCATAAACAAGGTTCAGCTTTCTACTTTGGCTTATTCGACCCGGACGAAAAAGAGATTATTGGCGTTGCCAATTTTTCCAATGTTGTTCGTGGCTCTTTTCATGCCTGCTATC  >  NC_000913/1125588‑1125880

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: