Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 27,887 G→T 5.9% D199Y (GAC→TAC)  rihC → ribonucleoside hydrolase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091327,8870GT 5.9% 239.2 / 2.5 68D199Y (GAC→TAC) rihCribonucleoside hydrolase 3
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGG  >  NC_000913/27750‑28025
                                                                                                                                         |                                                                                                                                          
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          <  2:1146701/140‑1 (MQ=255)
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          <  1:1066924/140‑1 (MQ=255)
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          <  1:1658212/140‑1 (MQ=255)
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          >  1:1502407/1‑140 (MQ=255)
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          >  1:1551000/1‑140 (MQ=255)
ccGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATAGTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGAc                                                                                                                                          <  1:365413/140‑1 (MQ=255)
  ggACGCGGCAACTGTACGCCACACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGATATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTa                                                                                                                                        <  2:1352878/140‑1 (MQ=255)
  ggACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTa                                                                                                                                        <  1:985049/140‑1 (MQ=255)
   gACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTAt                                                                                                                                       <  1:1633157/140‑1 (MQ=255)
       gggCAACGGTACGCCACACGCCGAGTTTAATATTGCTGCCGATCCAGTAGCTGCTGCCTGTGTCTTCCGCAGTGGAATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTAtctct                                                                                                                                   <  2:1685532/139‑1 (MQ=255)
        ggCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTAtctctc                                                                                                                                  >  2:1545504/1‑140 (MQ=255)
                tACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGc                                                                                                                          <  1:1211492/140‑1 (MQ=255)
                tACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGc                                                                                                                          <  1:294977/140‑1 (MQ=255)
                  cGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCg                                                                                                                        >  1:455070/1‑140 (MQ=255)
                        aCGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTa                                                                                                                  >  2:694991/1‑140 (MQ=255)
                                tttAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTAc                                                                                                          <  2:823169/140‑1 (MQ=255)
                                       ttGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGaaa                                                                                                   >  1:1398754/1‑140 (MQ=255)
                                          ctgcCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATg                                                                                                >  2:1485601/1‑140 (MQ=255)
                                          ctgcCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATg                                                                                                >  1:1546807/1‑140 (MQ=255)
                                          ctgcCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATg                                                                                                >  1:1546509/1‑140 (MQ=255)
                                                  ccAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGc                                                                                        <  1:597702/140‑1 (MQ=255)
                                                  ccAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGc                                                                                        <  2:1370091/140‑1 (MQ=255)
                                                   cAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAATGCTTCACGccc                                                                                      <  1:1542409/140‑1 (MQ=255)
                                                    aGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                                                      <  1:594209/140‑1 (MQ=255)
                                                    aGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                                                      <  1:666345/140‑1 (MQ=255)
                                                    aGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                                                      <  2:303176/140‑1 (MQ=255)
                                                           gctgcCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTAGGATGTCACCAAACAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTAACGGGAAAATTCTTCACGCCCTGTTTAg                                                                               >  2:716181/1‑140 (MQ=255)
                                                            ctgcCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGc                                                                              <  2:519162/140‑1 (MQ=255)
                                                            ctgcCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGc                                                                              >  1:835899/1‑140 (MQ=255)
                                                            ctgcCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGc                                                                              <  2:1100026/140‑1 (MQ=255)
                                                                       ttCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAg                                                                   <  2:608205/140‑1 (MQ=255)
                                                                         ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTTACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCg                                                                 >  2:560839/1‑140 (MQ=255)
                                                                         ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTTACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCg                                                                 >  2:561068/1‑140 (MQ=255)
                                                                            cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                                                              <  2:1435556/140‑1 (MQ=255)
                                                                            cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                                                              <  1:1421318/140‑1 (MQ=255)
                                                                            cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                                                              <  2:989087/140‑1 (MQ=255)
                                                                                 gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                                                         >  2:1069345/1‑140 (MQ=255)
                                                                                 gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                                                         >  2:1069080/1‑140 (MQ=255)
                                                                                 gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                                                         >  2:1069007/1‑140 (MQ=255)
                                                                                      gAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAg                                                    >  2:1524599/1‑140 (MQ=255)
                                                                                      gAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAg                                                    >  2:1412912/1‑140 (MQ=255)
                                                                                         aTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                                                 >  2:467916/1‑140 (MQ=255)
                                                                                            gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt                                              >  1:1028290/1‑140 (MQ=255)
                                                                                            gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt                                              >  1:637807/1‑140 (MQ=255)
                                                                                            gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt                                              >  2:969627/1‑140 (MQ=255)
                                                                                            gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt                                              >  2:969855/1‑140 (MQ=255)
                                                                                             tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCt                                                                                              <  2:1654977/91‑1 (MQ=255)
                                                                                             tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCt                                                                                              >  1:1654977/1‑91 (MQ=255)
                                                                                               aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTTACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg                                           <  1:561068/140‑1 (MQ=255)
                                                                                               aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTTACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg                                           <  1:560839/140‑1 (MQ=255)
                                                                                                      gTTTGGCTGTCACCTCTCAGGCAATATTAACTCCTGACTAGCTCTCTACACTGCCGCAGTTAAACCGTTCCGGGAAAATGCTTCAGGCCCTTTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCAc                                    <  1:838070/140‑1 (MQ=255)
                                                                                                      gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCAc                                    <  2:1514207/140‑1 (MQ=255)
                                                                                                      gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCAc                                    <  2:20846/140‑1 (MQ=255)
                                                                                                           gATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAATATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtct                               >  2:373279/1‑140 (MQ=255)
                                                                                                           gATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtct                               >  2:1083867/1‑140 (MQ=255)
                                                                                                           gATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtct                               >  2:373341/1‑140 (MQ=255)
                                                                                                           gATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtct                               >  2:1150897/1‑140 (MQ=255)
                                                                                                           gATGTCACCAATCAGGCAATAATAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtct                               >  2:373340/1‑140 (MQ=255)
                                                                                                            ttGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGGCGCAGTTAAACCGTTCCGGGAAAATGCGTCACGGCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGAtctc                              <  2:1652409/139‑1 (MQ=255)
                                                                                                                 aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTgcgc                         >  1:349964/1‑140 (MQ=255)
                                                                                                                 aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTgcgc                         >  1:1042836/1‑140 (MQ=255)
                                                                                                                  ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGcc                        <  2:175707/140‑1 (MQ=255)
                                                                                                                           cAATCTTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCctgg               >  2:276142/1‑140 (MQ=255)
                                                                                                                                ttAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGt          <  2:1025221/140‑1 (MQ=255)
                                                                                                                                   aCTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTgcg       >  1:1223461/1‑140 (MQ=255)
                                                                                                                                    cTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTgcgc      >  1:1560762/1‑140 (MQ=255)
                                                                                                                                        tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCgg  >  2:433365/1‑140 (MQ=255)
                                                                                                                                        tGACTATCTCTCTACACTGCCGCAGATAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCgg  >  1:969001/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CCGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGG  >  NC_000913/27750‑28025

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