Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 326,821 C→T 5.2% G143S (GGC→AGC)  betA ← choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913326,8210CT 5.2% 256.8 / 2.7 80G143S (GGC→AGC) betAcholine dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base T (2/2);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGG  >  NC_000913/326683‑326957
                                                                                                                                          |                                                                                                                                        
gATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAccgcc                                                                                                                                 <  8:22183/148‑1 (MQ=255)
 aTCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAccgccg                                                                                                                                >  2:69345/1‑148 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCTATCAc                                                                                                                                     >  7:417622/1‑142 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCTATCAc                                                                                                                                     <  8:417622/142‑1 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCaa                                                                                                                                     <  6:424619/142‑2 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAc                                                                                                                                     <  4:25998/142‑1 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAc                                                                                                                                     >  3:25998/1‑142 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAc                                                                                                                                     >  5:424619/1‑142 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGt                                                                                                                               >  3:276837/1‑148 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGt                                                                                                                               >  1:323446/1‑148 (MQ=255)
  tCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGt                                                                                                                               >  1:317462/1‑148 (MQ=255)
   ccATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGt                                                                                                                               >  5:231265/1‑147 (MQ=255)
    cATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTg                                                                                                                              <  6:231265/147‑1 (MQ=255)
    cATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTg                                                                                                                              <  8:57791/147‑1 (MQ=255)
     aTCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGa                                                                                                                             >  3:273341/1‑147 (MQ=255)
      tCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGa                                                                                                                             >  5:173446/1‑146 (MQ=255)
        ggACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGt                                                                                                                         >  7:239337/1‑148 (MQ=255)
          aCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCg                                                                                                                       >  8:247455/1‑148 (MQ=255)
          aCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCg                                                                                                                       <  2:390568/148‑1 (MQ=255)
             aaaaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCTATCACCGCCGTGATAGTCGttt                                                                                                                    >  6:405277/1‑148 (MQ=255)
             aaaaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGttt                                                                                                                    >  8:309115/1‑148 (MQ=255)
              aaaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGtttt                                                                                                                   >  4:351872/1‑148 (MQ=255)
               aaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGtttt                                                                                                                   >  3:240590/1‑147 (MQ=255)
                  cTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCAcc                                                                                                               <  1:123029/148‑1 (MQ=255)
                         ctgATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCatat                                                                                                          <  4:240590/146‑1 (MQ=255)
                            aTAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATcgc                                                                                                       >  7:366263/1‑146 (MQ=255)
                               aCCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtc                                                                                                  >  2:54987/1‑148 (MQ=255)
                                         tCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc                                                                                        <  1:364405/148‑1 (MQ=255)
                                         tCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc                                                                                        <  7:73887/148‑1 (MQ=255)
                                         tCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc                                                                                        <  7:309115/148‑1 (MQ=255)
                                           gTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCg                                                                                       <  1:69345/147‑1 (MQ=255)
                                            tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCTATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgt                                                                                     <  5:405277/148‑1 (MQ=255)
                                            tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgt                                                                                     <  8:366263/148‑1 (MQ=255)
                                            tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgt                                                                                     <  3:351872/148‑1 (MQ=255)
                                             ccttGCGCGGGGAGCCCGCCTGCACGCCCGCTTCAATCATCGCTGCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCGCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTGGCGgta                                                                                    <  6:331385/145‑1 (MQ=255)
                                                 gcgcGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGtagtag                                                                                <  8:239337/148‑1 (MQ=255)
                                                  cgcgGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGtagtag                                                                                <  7:247455/147‑1 (MQ=255)
                                                    cgGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggc                                                                             <  1:54987/148‑1 (MQ=255)
                                                     gggTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggca                                                                            <  4:44620/148‑1 (MQ=255)
                                                      ggTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggcag                                                                           <  4:276837/148‑1 (MQ=255)
                                                       gctacccGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgt                                                                                     <  5:161717/133‑1 (MQ=255)
                                                       gctacccGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgt                                                                                     >  6:161717/5‑137 (MQ=255)
                                                       gTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATcgcg                                                                                                      <  8:286375/120‑1 (MQ=255)
                                                       gTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATcgcg                                                                                                      >  7:286375/1‑120 (MQ=255)
                                                        tAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggcaggc                                                                         <  4:341326/148‑1 (MQ=255)
                                                             cGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCg                                                                    <  6:316888/148‑1 (MQ=255)
                                                              gCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTc                                                                     >  8:46684/1‑146 (MQ=255)
                                                                   cACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCg                                                                                                                       <  8:158176/91‑1 (MQ=255)
                                                                   cACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCg                                                                                                                       >  7:158176/1‑91 (MQ=255)
                                                                         cGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAg                                                        >  3:110233/1‑148 (MQ=255)
                                                                          gCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGt                                                       <  4:273341/148‑1 (MQ=255)
                                                                             tCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTtct                                                    >  2:123979/1‑148 (MQ=255)
                                                                               aaTCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTcc                                                  <  1:357338/148‑1 (MQ=255)
                                                                                 tcatcGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAg                                                <  2:59387/148‑1 (MQ=255)
                                                                                  catcGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGa                                               >  8:389633/1‑148 (MQ=255)
                                                                                  catcGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGAGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCAACAGTTATCCAgc                                               >  4:306279/1‑147 (MQ=255)
                                                                                    tcGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTg                                                                                                                              >  5:90179/1‑67 (MQ=255)
                                                                                    tcGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTg                                                                                                                              <  6:90179/67‑1 (MQ=255)
                                                                                       cTTCAAACAGCGTATTGACTCCGTGTTTGGAGGTAGTGACGCTCACCGGGCGATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCggg                                          <  8:108836/148‑1 (MQ=255)
                                                                                          cAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc                                                                                        >  7:221121/1‑99 (MQ=255)
                                                                                          cAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc                                                                                        <  8:221121/99‑1 (MQ=255)
                                                                                           aaaCAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCt                                      <  1:141723/148‑1 (MQ=255)
                                                                                             aCAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGAGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTg                                    <  2:317462/148‑1 (MQ=255)
                                                                                             aCAGCGGATTGACGCCGGGGTTGGAGGTCGTGACGCTCACCGGGCCATCCCCGCCGGGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTg                                    <  2:323446/148‑1 (MQ=255)
                                                                                              cAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgtagta                                                                                 <  7:34135/102‑1 (MQ=255)
                                                                                              cAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgtagta                                                                                 >  8:34135/1‑102 (MQ=255)
                                                                                                gCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCtt                                     <  7:202346/144‑1 (MQ=255)
                                                                                                gCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCtt                                     >  8:202346/1‑144 (MQ=255)
                                                                                                      tGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGt                            >  3:124374/1‑147 (MQ=255)
                                                                                                        aCGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTAt                         >  5:225175/1‑148 (MQ=255)
                                                                                                             gggTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTtctc                                                   <  3:197092/117‑1 (MQ=255)
                                                                                                                    gAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAgc              <  6:173446/147‑1 (MQ=255)
                                                                                                                      ggTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCGCCAGCCCGGGTTCTTGCGCCCAGTTATCGAGATCCCGCGCa           >  4:301848/1‑148 (MQ=255)
                                                                                                                           tGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGcc      <  4:110233/148‑1 (MQ=255)
                                                                                                                            gACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTcc                                                  >  3:125350/1‑103 (MQ=255)
                                                                                                                            gACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTcc                                                  <  4:125350/103‑1 (MQ=255)
                                                                                                                             gcgcgcgccGGGCCCTCACCGGCGTGATAGGCGGTGTCCGGCATTTCGCGGGGCTCGGCGTTGCGGTAGTAGGGCAGGCAGTCGCGGTGGCGCCAGGGCTCCAGACCGGGTGCTTGCGCCCAGTGAGCGCGATCCAGCGCTTTGCCAc    >  4:30783/8‑148 (MQ=255)
                                                                                                                               gCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACgg  <  7:389633/148‑1 (MQ=255)
                                                                                                                                   aCCGGGCCATCTCCGCCGTGATAGTTGTTTTCACCCATATGGCGAGTCTCGGCCTTGCGGTAGTAGGGCTGGCAGTCGAGGTAGCTCCTGTTCTCCAGACCGGGtt                                        <  2:401068/106‑1 (MQ=255)
                                                                                                                                    ccGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTc                                       >  1:401068/1‑106 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGG  >  NC_000913/326683‑326957

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