Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,407,134 G→T 5.6% G161C (GGC→TGC)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,1340GT 5.6% 262.4 / 2.5 71G161C (GGC→TGC) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base G (26/41);  new base T (2/2);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGT  >  NC_000913/4406999‑4407273
                                                                                                                                       |                                                                                                                                           
aaaaCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGt                                                                                                                                         <  1:1472501/140‑1 (MQ=255)
       gCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCttg                                                                                                                                  >  2:294731/1‑139 (MQ=255)
           tCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTgg                                                                                                                              <  2:928871/140‑1 (MQ=255)
             aTGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGtt                                                                                                                            >  1:428597/1‑140 (MQ=255)
              tGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGATGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTc                                                                                                                           <  1:1514017/140‑1 (MQ=255)
                       ggTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACa                                                                                                                  <  1:328520/140‑1 (MQ=255)
                         tGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGc                                                                                                                <  2:802322/140‑1 (MQ=255)
                           cTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCg                                                                                                              >  1:823949/1‑140 (MQ=255)
                           cTGGATCAGCCGCGGGGTGCTGACCGTAAAGTTAGTATTGTAGCGCGTATTGTCCTCGTACTGTTGCCAAAAAACAGCCTGATTGTTGGTCGCTTCTTTACCGTAGCGGGCGTCGGCTTTGTGGTTCCTGACGTCAGCCg                                                                                                              >  1:1011901/1‑140 (MQ=255)
                                 cAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAg                                                                                                        <  1:518115/140‑1 (MQ=255)
                                     cGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTc                                                                                                    <  2:768611/140‑1 (MQ=255)
                                       cTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGa                                                                                                  <  1:105728/140‑1 (MQ=255)
                                       cTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGa                                                                                                  <  2:1132895/140‑1 (MQ=255)
                                       cTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGa                                                                                                  <  2:1132975/140‑1 (MQ=255)
                                       cTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGa                                                                                                  <  2:1772051/140‑1 (MQ=255)
                                       cTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGa                                                                                                  <  1:108855/140‑1 (MQ=255)
                                              cTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTa                                                                                           >  1:1474971/1‑140 (MQ=255)
                                                    gTAAAGGTCGTCGTGAAGCGCGTAGTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATGGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTGTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATCTCTTAATcccg                                                                                     >  2:171443/1‑140 (MQ=255)
                                                      aaaGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGGTTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCGTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATcccgcc                                                                                   >  2:197708/1‑140 (MQ=255)
                                                      aaaGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATcccgcc                                                                                   >  2:1848669/1‑140 (MQ=255)
                                                        aGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATcccgcccg                                                                                 <  2:1055424/140‑1 (MQ=255)
                                                          gtcgtcGTTATGAGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGATACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  2:172480/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTTGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  1:867345/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  1:421313/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  2:172641/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  2:172625/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  2:172406/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  1:408492/140‑1 (MQ=255)
                                                          gtcgtcGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTGGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgat                                                                               <  1:421068/140‑1 (MQ=255)
                                                             gtcgtGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCTTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTACGATATCTTAATCCCGCCCgatcag                                                                            >  1:1862712/1‑140 (MQ=255)
                                                             gtcgtGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgatcag                                                                            >  1:1507091/1‑140 (MQ=255)
                                                             gtcgtGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgatcag                                                                            >  1:1507012/1‑140 (MQ=255)
                                                             gtcgtGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgatcag                                                                            >  2:1313892/1‑140 (MQ=255)
                                                                   aaGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCCGATCATg                                                                      >  1:133316/1‑140 (MQ=255)
                                                                     gcgcGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATggg                                                                    >  2:73983/1‑140 (MQ=255)
                                                                       gcgTATTGTCCGCGTACTGGTGCCAAAAACCCGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCATAATCCCGCCCGATCAGATCATGGgcg                                                                  <  2:984778/140‑1 (MQ=255)
                                                                       gcgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTACTCCCGCCCGATCAGATCATGGgcg                                                                  <  1:1446582/140‑1 (MQ=255)
                                                                       gcgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  <  1:1739998/140‑1 (MQ=255)
                                                                       gcgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  <  2:131518/140‑1 (MQ=255)
                                                                       gcgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  <  2:984779/140‑1 (MQ=255)
                                                                        cgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGTGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  >  2:682454/1‑139 (MQ=255)
                                                                        cgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGTGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  >  2:682415/1‑139 (MQ=255)
                                                                        cgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGTGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  <  1:682454/139‑1 (MQ=255)
                                                                        cgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGTGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcg                                                                  <  1:682415/139‑1 (MQ=255)
                                                                        cgTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcgc                                                                 >  1:1572697/1‑140 (MQ=255)
                                                                         gTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcgcg                                                                <  2:1754252/140‑1 (MQ=255)
                                                                         gTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGgcgcg                                                                <  2:1756590/140‑1 (MQ=255)
                                                                              gTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGAt                                                           <  1:841516/140‑1 (MQ=255)
                                                                              gTCCGCATACTGGTGCCAAACACCAGCCAGATTGTTGGTCGCTACTTTACCGCAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGAt                                                           <  2:828941/140‑1 (MQ=255)
                                                                                   cgTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGTTCAGATCATGGGCGCGCGGATGGGCt                                                      >  2:899374/1‑140 (MQ=255)
                                                                                   cgTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCt                                                      >  1:503499/1‑140 (MQ=255)
                                                                                      aCTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTtg                                                   <  2:992130/140‑1 (MQ=255)
                                                                                      aCTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTtg                                                   <  1:690976/140‑1 (MQ=255)
                                                                                      aCTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTtg                                                   <  2:904377/140‑1 (MQ=255)
                                                                                      aCTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTtg                                                   <  1:1731031/140‑1 (MQ=255)
                                                                                         ggTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTgg                                                >  1:457/1‑140 (MQ=255)
                                                                                          gTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGt                                               <  1:417737/140‑1 (MQ=255)
                                                                                             ccAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGt                                            >  2:46862/1‑140 (MQ=255)
                                                                                                    ccagccagATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGaa                                     <  2:240342/140‑1 (MQ=255)
                                                                                                    ccagccagATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGaa                                     <  2:1643582/140‑1 (MQ=255)
                                                                                                          agATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACt                               <  1:73983/140‑1 (MQ=255)
                                                                                                                       cTACTTTACCGAAGCAGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTc                  >  1:1789175/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  2:112699/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  2:649775/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  2:649843/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  2:1538464/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  1:118742/1‑140 (MQ=255)
                                                                                                                                  aaGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCaaa       >  1:55882/1‑140 (MQ=255)
                                                                                                                                     cGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCg    <  1:860049/140‑1 (MQ=255)
                                                                                                                                     cGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCg    <  2:1708916/140‑1 (MQ=255)
                                                                                                                                     cGGGCGTCGGCATTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCATAATCCCGCCCGAACAGATCATGGGCGCGCGGATGGGCATTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCg    <  2:248349/140‑1 (MQ=255)
                                                                                                                                       ggCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCTATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCTGt  >  1:1274144/1‑140 (MQ=255)
                                                                                                                                       ggCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCTATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCTGt  >  1:1274007/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                           
AAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTGACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCTTTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGT  >  NC_000913/4406999‑4407273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: