Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,311,107 C→G 5.5% G21A (GGC→GCC)  pnp ← polynucleotide phosphorylase/polyadenylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,311,1070CG 5.5% 254.4 / 5.8 73G21A (GGC→GCC) pnppolynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (2/2);  total (38/35)
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.79e-01

CGCTCCTGATAGTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAATGACA  >  NC_000913/3310972‑3311254
                                                                                                                                       |                                                                                                                                                   
cGCTCCTGATAGTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCagag                                                                                                                                         <  5:65806/148‑1 (MQ=255)
      tGATAGTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACgg                                                                                                                                   >  7:251382/1‑148 (MQ=255)
          aGTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGgtgtg                                                                                                                               >  7:164005/1‑148 (MQ=255)
           gTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGgtgtgt                                                                                                                              <  3:6250/148‑1 (MQ=255)
            atAACGGTCAGTGGGAAGAAGGCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCTGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGtt                                                                                                                             <  4:41815/147‑1 (MQ=255)
             tAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTcc                                                                                                                                             >  4:268118/1‑131 (MQ=255)
             tAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTcc                                                                                                                                             <  3:268118/131‑1 (MQ=255)
              aaCGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTgg                                                                                                                           >  8:240952/1‑148 (MQ=255)
                  gTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCg                                                                                                                        >  4:67089/1‑147 (MQ=255)
                    cAGTGGGAAGAAGTCCTGACCTCGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTAc                                                                                                                     >  1:214489/1‑148 (MQ=255)
                                    tGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGa                                                                                                                 >  4:282127/1‑136 (MQ=255)
                                    tGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGa                                                                                                                 <  3:282127/136‑1 (MQ=255)
                                     gACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGAt                                                                                                    >  5:45586/1‑148 (MQ=255)
                                     gACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGAt                                                                                                    >  8:52875/1‑148 (MQ=255)
                                      aCCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATc                                                                                                   <  8:164005/148‑1 (MQ=255)
                                         tGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGa                                                                                                >  7:88590/1‑148 (MQ=255)
                                            tttGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGAtt                                                                                              >  7:174434/1‑147 (MQ=255)
                                            tttGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTa                                                                                             <  5:11363/148‑1 (MQ=255)
                                              tGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTaa                                                                                            <  6:143545/147‑1 (MQ=255)
                                                 cTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCaa                                                                                        <  2:202473/148‑1 (MQ=255)
                                                 cTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCaa                                                                                        >  8:68251/1‑148 (MQ=255)
                                                      ttCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTa                                                                                    <  7:240952/147‑1 (MQ=255)
                                                              aacaacGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGACCGATCGGATTAAGCAATGTAATATCcttt                                                                           >  6:98274/1‑148 (MQ=255)
                                                              aacaacGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGt                                                                                     >  7:261679/1‑138 (MQ=255)
                                                              aacaacGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGt                                                                                     <  8:261679/138‑1 (MQ=255)
                                                                caacGGTAACGAATACCGCGGTGTGATCGATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCTTGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTAAGAACGATCGGAt                                                                                               >  5:52884/1‑126 (MQ=255)
                                                                  acGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctct                                                                       >  6:260595/1‑148 (MQ=255)
                                                                   cGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctt                                                                      <  5:260595/148‑1 (MQ=255)
                                                                     gTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    <  7:229596/148‑1 (MQ=255)
                                                                       aaCGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctt                                                                            <  6:228892/138‑1 (MQ=255)
                                                                       aaCGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctt                                                                            >  5:228892/1‑138 (MQ=255)
                                                                           aaTACCGGGGTGTCATCCATGCTAACCATAACAGCGGGAGTAGGCTGACGAGGCATCATGCCGGTTTCGAGAGTGTCGGTGTGTTGGCCGTACTGGAATTTTCGAACGATCGGTTTAAGCAATGTAATATCCTTTCTCTTTCTTAGAc                                                              <  7:283460/148‑1 (MQ=255)
                                                                                cgcgGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTAc                                                         <  8:174434/148‑1 (MQ=255)
                                                                                cgcgGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTAc                                                         >  3:57515/1‑148 (MQ=255)
                                                                                 gcgGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAg                                                                <  8:397412/140‑1 (MQ=255)
                                                                                 gcgGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAg                                                                >  7:397412/1‑140 (MQ=255)
                                                                                     tgtCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTAc                                                    <  6:363544/148‑1 (MQ=255)
                                                                                        cATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAatat                                                                                <  7:221039/117‑1 (MQ=255)
                                                                                        cATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAatat                                                                                >  8:221039/1‑117 (MQ=255)
                                                                                          tCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCAc                                               >  7:142996/1‑148 (MQ=255)
                                                                                           ccATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACt                                              >  7:413690/1‑148 (MQ=255)
                                                                                            cATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATTTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTg                                             <  2:214489/148‑1 (MQ=255)
                                                                                               gCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAg                                                                >  8:166861/1‑126 (MQ=255)
                                                                                               gCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAg                                                                <  7:166861/126‑1 (MQ=255)
                                                                                               gCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGgtg                                          >  1:394991/1‑148 (MQ=255)
                                                                                               gCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGgtg                                          >  1:394995/1‑148 (MQ=255)
                                                                                                cTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGgtgt                                         <  2:87210/148‑1 (MQ=255)
                                                                                                  aaCCATAACAGCGGCAGTACCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGAACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGtt                                        >  3:181539/1‑147 (MQ=255)
                                                                                                   aCCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       <  1:230525/147‑1 (MQ=255)
                                                                                                   aCCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       <  3:67089/147‑1 (MQ=255)
                                                                                                   aCCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       <  5:105412/147‑1 (MQ=255)
                                                                                                    ccATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATc                                                                               >  2:361701/1‑106 (MQ=255)
                                                                                                    ccATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATc                                                                               <  1:361701/106‑1 (MQ=255)
                                                                                                      aTAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       <  8:26553/144‑1 (MQ=255)
                                                                                                      aTAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       >  7:26553/1‑144 (MQ=255)
                                                                                                      aTAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATAc                                   <  8:88590/148‑1 (MQ=255)
                                                                                                        aaCAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       >  8:360673/1‑142 (MQ=255)
                                                                                                        aaCAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTa                                       <  7:360673/142‑1 (MQ=255)
                                                                                                               gCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGAtcttc                          <  5:323588/148‑1 (MQ=255)
                                                                                                                    aGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTgcg                      >  1:305234/1‑147 (MQ=255)
                                                                                                                     gCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTgcgc                     <  4:57515/147‑1 (MQ=255)
                                                                                                                      ccTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    <  1:373521/99‑1 (MQ=255)
                                                                                                                      ccTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    >  2:373521/1‑99 (MQ=255)
                                                                                                                       cTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATc                  >  4:347464/1‑148 (MQ=255)
                                                                                                                          aCGAGCCATCATGGAGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    <  5:176944/95‑1 (MQ=255)
                                                                                                                           cGAGCCATCATGGCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACTAAAGCTCGGATTAAGCAATGTAATATCctttctctttat                                                                   >  6:176944/1‑93 (MQ=255)
                                                                                                                           cGAGCCATCATGGCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    <  5:182722/94‑1 (MQ=255)
                                                                                                                           cGAGCCATCATGGCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCctttctctttc                                                                    >  6:182722/1‑94 (MQ=255)
                                                                                                                             aGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTcgcg            <  8:233181/148‑1 (MQ=255)
                                                                                                                             aGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTcgc             >  7:233181/1‑147 (MQ=255)
                                                                                                                                   catGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAAt      >  7:146863/1‑148 (MQ=255)
                                                                                                                                   catGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAAt      >  2:213896/1‑148 (MQ=255)
                                                                                                                                      gCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAATGa    >  8:62738/1‑147 (MQ=255)
                                                                                                                                       ccGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAATGACa  <  7:68251/148‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
CGCTCCTGATAGTTAACGGTCAGTGGGAAGAAGTCCTGACCTGGTTTGGCTTTTTTCTGGCCAACAACGGTAACGAATACCGCGGTGTCATCCATGCTAACCATAACAGCGGCAGTAGCCTGACGAGCCATCATGCCGGTTTCCAGAGTCACGGTGTGTTGGCCGTACTGGAATTTACGAACGATCGGATTAAGCAATGTAATATCCTTTCTCTTTCTTAGACAGTACCTTACGGCACTGGTGTTAATACCCGATCTTCTGCGCATCCTCGCGACTAATGACA  >  NC_000913/3310972‑3311254

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