Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,177,447 C→T 5.2% G155D (GGT→GAT)  nudF ← ADP‑ribose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,4470CT 5.2% 248.6 / 2.1 77G155D (GGT→GAT) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (36/37);  new base T (2/2);  total (38/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 = 4.41e-01

ACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATC  >  NC_000913/3177301‑3177592
                                                                                                                                                  |                                                                                                                                                 
ccGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCg                                                                                                                                                 >  4:191748/2‑149 (MQ=255)
   cTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCtt                                                                                                                                              >  8:86181/1‑149 (MQ=255)
       atgatgCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCg                                                                                                                                           <  7:165316/148‑1 (MQ=255)
       atgatgCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCg                                                                                                                                           <  7:200747/148‑1 (MQ=255)
        tgatgCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATATCGCTTGCGGt                                                                                                                                         >  3:31984/1‑149 (MQ=255)
        tgatgCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATATCGCTTGCGGt                                                                                                                                         >  3:27692/1‑149 (MQ=255)
          atgCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTc                                                                                                                                        <  6:306366/148‑1 (MQ=255)
             cagctgcagcCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTg                                                                                                                                     <  4:153869/148‑1 (MQ=255)
              agctgcagcCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGc                                                                                                                                   >  1:11930/1‑149 (MQ=255)
                 tgcagcCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTc                                                                                                                                >  2:25557/1‑149 (MQ=255)
                     gcCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACt                                                                                                                            >  3:174644/1‑149 (MQ=255)
                     gcCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCACTCTTCTACCCACTGGTATGCATGTTCCCTGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGACAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACt                                                                                                                            >  7:239354/1‑149 (MQ=255)
                     gcCATAGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCACTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGAGAATACCGCTTGCGGTCGTGGCGTCCACt                                                                                                                            >  3:17995/1‑149 (MQ=255)
                      ccATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACtt                                                                                                                           >  7:297870/1‑149 (MQ=255)
                       cATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATATCGCTTGCGGTCGTGGCGTCCACTTc                                                                                                                          <  4:27692/149‑1 (MQ=255)
                       cATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATATCGCTTGCGGTCGTGGCGTCCACTTc                                                                                                                          <  4:31984/149‑1 (MQ=255)
                         ttGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGAGTTCGTGGCGTCCACTTCGc                                                                                                                        >  6:203449/1‑149 (MQ=255)
                         ttGCAAAGCGATGACCGACGCTGCGTCGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGc                                                                                                                        >  1:54288/1‑149 (MQ=255)
                            cAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCGTCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGccc                                                                                                                      <  8:214326/148‑1 (MQ=255)
                             aaaGCGATGACCGACGCTGCGTTGTCGATTATCCCCTCTTCTACCCACTGGTATGCCTGTTACCGGAATACCACATGAACGCGAATATCATCGTTTTCATCAGCCAGACAGGGAATACCGCTTGCGCTCGTTGCGTCCACTTAGccaa                                                                                                                     >  5:95897/1‑146 (MQ=255)
                               aGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCccacc                                                                                                                   >  4:97315/1‑148 (MQ=255)
                                  gATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATaa                                                                                                               >  2:199491/1‑149 (MQ=255)
                                      aCCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGa                                                                                                           <  3:97315/149‑1 (MQ=255)
                                      aCCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGa                                                                                                           <  1:44489/149‑1 (MQ=255)
                                          aCGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCccacca                                                                                                                  <  5:181811/138‑1 (MQ=255)
                                          aCGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCccacca                                                                                                                  >  6:181811/1‑138 (MQ=255)
                                                cGTTGTCGATTTTCCCCTCTTCTAGCCACTGGTATGCCTGTTCCCGGCTTTCCACATGAACGCGAATATCTTCGGTTTCATCTGGCAGACAGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCt                                                                                                   <  6:37680/147‑1 (MQ=255)
                                                cGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGCCCGTGAATACCGCTTGCTGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCt                                                                                                   <  6:74905/147‑1 (MQ=255)
                                                    gTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTg                                                                                              <  1:199491/148‑1 (MQ=255)
                                                     tCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGt                                                                                            >  8:87546/1‑149 (MQ=255)
                                                         ttttCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGccc                                                                                        <  7:86181/149‑1 (MQ=255)
                                                              cccTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCg                                                                                    <  2:11930/148‑1 (MQ=255)
                                                               ccTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCggg                                                                                  >  1:50308/1‑149 (MQ=255)
                                                                 tcttctACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCt                                                                                >  5:131468/1‑149 (MQ=255)
                                                                   ttctACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTg                                                                              >  2:262615/1‑149 (MQ=255)
                                                                      tACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGAAGAACGCTCACTGGTGCCCCACGGGCTTGc                                                                             >  6:183904/1‑147 (MQ=255)
                                                                        cctaCTGGTATGCCTGTTCCGGGCTTACCACATGAACGCGAATTTCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGTCTTGCCAg                                                                          <  6:95897/145‑1 (MQ=255)
                                                                        cccACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGcc                                                                            <  1:25557/146‑1 (MQ=255)
                                                                        cccACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCa                                                                           <  7:87546/147‑1 (MQ=255)
                                                                        cccACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTATCATCAGCCAGACCGTGAATACCGCTTCCGGTCGTGGCGTACACTTCGCCCACCATAATTGAAGAACGCt                                                                                                   >  6:217297/1‑123 (MQ=255)
                                                                         ccACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGa                                                                        <  3:191748/149‑1 (MQ=255)
                                                                           aCTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGaa                                                                       >  6:14088/1‑148 (MQ=255)
                                                                            cTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAAc                                                                     <  7:181809/149‑1 (MQ=255)
                                                                               gTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGACCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACtt                                                                   <  5:203449/148‑1 (MQ=255)
                                                                                 aTGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCa                                                                           >  4:37101/1‑138 (MQ=255)
                                                                                 aTGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCa                                                                           <  3:37101/138‑1 (MQ=255)
                                                                                   gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTg                                                                                                            >  4:174534/1‑103 (MQ=255)
                                                                                   gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTg                                                                                                            <  3:174534/103‑1 (MQ=255)
                                                                                   gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAacac                                                              >  1:152415/1‑149 (MQ=255)
                                                                                   gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAacac                                                              >  3:192504/1‑149 (MQ=255)
                                                                                       gTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                          >  1:242364/1‑149 (MQ=255)
                                                                                                  acgacaTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACAGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                <  4:17995/145‑1 (MQ=255)
                                                                                                  aCCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                <  4:174644/148‑1 (MQ=255)
                                                                                                   ccACATTCACGCGAATATCTTCGTTTTCATCAGCCAGACCGTTAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               <  8:239354/148‑1 (MQ=255)
                                                                                                   ccACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               <  3:94767/148‑1 (MQ=255)
                                                                                                   ccACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               <  1:262615/148‑1 (MQ=255)
                                                                                                    cacaTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                             >  2:279265/1‑149 (MQ=255)
                                                                                                    cacaTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                             <  1:279265/149‑1 (MQ=255)
                                                                                                    cacaTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                             <  4:192504/149‑1 (MQ=255)
                                                                                                      caTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                            >  5:235267/1‑148 (MQ=255)
                                                                                                        tGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                         <  2:54288/149‑1 (MQ=255)
                                                                                                           aCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCAACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAg                                      <  5:14088/149‑1 (MQ=255)
                                                                                                                    tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               <  4:75258/131‑1 (MQ=255)
                                                                                                                    tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               >  3:75258/1‑131 (MQ=255)
                                                                                                                      ttCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               >  7:99382/1‑129 (MQ=255)
                                                                                                                      ttCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                               <  8:99382/129‑1 (MQ=255)
                                                                                                                        cGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctc                             <  5:58487/145‑1 (MQ=255)
                                                                                                                        cGTTTTCATCAGCCAGACCCTGAATACCGCGTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                            <  6:120638/146‑1 (MQ=255)
                                                                                                                                tcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc                 >  6:54189/1‑149 (MQ=255)
                                                                                                                                tcaGCCAGACCGTGAATACCGCATGCTTTCGTGGCGTCCACTTCTAACACAATAAATGACGAACGCTCACAGCAACCCCCCGCGCATCACAGGAAACTTAACAAATGATTTCTCCGATAGACTATCAGTACCGCCTCTTCAAACCCTTc                 >  4:143133/1‑149 (MQ=255)
                                                                                                                                tcaGCCAGACCGTGAATACCGCAAGCGGACCTAGCGACCACTTAGCACACCATAAATGACGAACCCTCACAGGTACCCCCCGCGCTTCCCAGCAAACATAACAACCGAATGCTCCGTATCACTATCACTCCCTCCACATCAAACCCTTc                 >  4:138122/1‑149 (MQ=255)
                                                                                                                                tcaGCCAGACCGTGAATACAGCTAGAGGTCGTGGCGTCCACTTCGCCCAACATAATTGACGAACTCTCACTGGTGCCCCCCGGCCTTGCCAGGAAACTTAACACCGGTTTGGTCAGTTTGACTATCAGTAACGCCTCTTCAATCGCTTc                 >  2:321369/1‑149 (MQ=255)
                                                                                                                                       gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACgg          <  2:152415/149‑1 (MQ=255)
                                                                                                                                         ccGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGc        <  8:297870/149‑1 (MQ=255)
                                                                                                                                          cGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGc        >  2:348769/1‑148 (MQ=255)
                                                                                                                                              aaTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACAAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACa    >  2:327943/1‑148 (MQ=255)
                                                                                                                                               aTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATc  <  2:242364/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATC  >  NC_000913/3177301‑3177592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑