Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,368 C→G 8.2% A369A (GCC→GCG ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,3680CG 8.2% 698.5 / 41.9 196A369A (GCC→GCGygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (85/95);  new base G (7/9);  total (92/104)
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.04e-01

Only 100 of 203 total aligned reads displayed.

CCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCG  >  NC_000913/3007221‑3007516
                                                                                                                                                   |                                                                                                                                                    
ccGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATcgccg                                                                                                                                                     >  3:1034542/1‑149 (MQ=255)
ccGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATcgccg                                                                                                                                                     >  2:95267/1‑149 (MQ=255)
       ggCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCgatgat                                                                                                                                              >  4:588902/1‑149 (MQ=255)
        gCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgata                                                                                                                                             <  3:1362882/149‑2 (MQ=255)
          gAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                                                                             >  4:541357/1‑147 (MQ=255)
          gAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                                                                             <  3:541357/147‑1 (MQ=255)
            aGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                                                                             >  1:1781168/1‑145 (MQ=255)
            aGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                                                                             <  2:1781168/145‑1 (MQ=255)
               ccTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTccc                                                                                                                                      <  1:425100/149‑1 (MQ=255)
                cTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGt                                                                                                                                         <  4:175298/145‑1 (MQ=255)
                cTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGt                                                                                                                                         >  3:175298/1‑145 (MQ=255)
                  gTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAGGAAGCCAGCTGGAACCCTTATATCATCGCCGTGATGATCCTGTCCCGTa                                                                                                                                   >  4:1598765/1‑149 (MQ=255)
                  gTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                                                                                             >  4:858104/1‑139 (MQ=255)
                       tGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                                                                                              <  1:1355506/149‑1 (MQ=255)
                       tGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                                                                                              >  2:1355506/1‑149 (MQ=255)
                          aCTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgatataag                                                                                                                                         <  1:808235/135‑6 (MQ=255)
                          aCTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                                                                                           >  2:635548/1‑149 (MQ=255)
                          aCTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                                                                                           >  1:1327957/1‑149 (MQ=255)
                            tgcctgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                                                                                              >  2:17907/1‑144 (MQ=255)
                              cctgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                                                                                       >  2:1679351/1‑149 (MQ=255)
                              cctgccAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                                                                                       >  2:1650335/1‑149 (MQ=255)
                                   cAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGt                                                                                                                  <  4:755581/149‑1 (MQ=255)
                                    aGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTg                                                                                                                 >  3:1541304/1‑149 (MQ=255)
                                       tGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATa                                                                                                                               <  1:276431/132‑1 (MQ=255)
                                       tGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGc                                                                                                                            >  2:1213977/1‑135 (MQ=255)
                                       tGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCAc                                                                                                              <  2:1327957/149‑1 (MQ=255)
                                       tGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCAc                                                                                                              <  3:845738/149‑1 (MQ=255)
                                              aGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGc                                                                                                           >  4:1546439/1‑145 (MQ=255)
                                              aGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                                                                                       <  3:308372/149‑1 (MQ=255)
                                               gcggcgTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAg                                                                                                      <  1:1679351/149‑1 (MQ=255)
                                                   cgTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                                                                                                  <  4:564955/149‑1 (MQ=255)
                                                     tATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                                                                                                    >  3:636390/1‑129 (MQ=255)
                                                     tATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                                                                                                <  4:494096/149‑1 (MQ=255)
                                                      aTCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                                                                               <  4:519334/149‑1 (MQ=255)
                                                      aTCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATAGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                                                                               <  1:1610778/149‑1 (MQ=255)
                                                         cTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCaa                                                                                                   <  4:1387315/142‑1 (MQ=255)
                                                          tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgatataag                                                                                                                                         <  2:1114431/103‑6 (MQ=255)
                                                          tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgatataag                                                                                                                                         >  1:1114431/1‑98 (MQ=255)
                                                          tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCGGCgatgat                                                                                                                                              <  1:1532748/98‑1 (MQ=255)
                                                          tGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                                                                                                    <  3:576571/108‑1 (MQ=255)
                                                           gTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                                                                          <  3:1607481/149‑1 (MQ=255)
                                                            tAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                                                                                         >  3:97458/1‑117 (MQ=255)
                                                            tAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                                                                          <  3:1678122/148‑1 (MQ=255)
                                                            tAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                                                                                          >  4:1678122/1‑148 (MQ=255)
                                                            tAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaag                                                                                         >  2:1088066/1‑149 (MQ=255)
                                                              aagaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCaa                                                                                                   >  3:595888/1‑137 (MQ=255)
                                                               agaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaag                                                                                      <  1:635548/149‑1 (MQ=255)
                                                               agaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaa                                                                                       >  1:1447800/1‑148 (MQ=255)
                                                                gaagGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                                                                                     <  3:157173/149‑1 (MQ=255)
                                                                      gAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGa                                                                               >  2:1420787/1‑149 (MQ=255)
                                                                       aaGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                                                                                                 <  2:783884/98‑1 (MQ=255)
                                                                         gTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgctg                                                                                             >  1:1481836/1‑132 (MQ=255)
                                                                          tGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                                                                               <  1:1186426/129‑1 (MQ=255)
                                                                          tGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgc                                                                           <  1:1650335/149‑1 (MQ=255)
                                                                           gACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg                                                                          >  2:146849/1‑149 (MQ=255)
                                                                              tGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                                                                               >  2:363769/1‑125 (MQ=255)
                                                                              tGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGa                                                                       >  4:1333526/1‑149 (MQ=255)
                                                                                aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGaaa                                                                     <  1:1088066/149‑1 (MQ=255)
                                                                                        ttCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTc                                                                                                                                        <  3:87121/74‑1 (MQ=255)
                                                                                          cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                                                                                         <  4:1720867/103‑1 (MQ=255)
                                                                                          cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                                                                                         >  3:1720867/1‑103 (MQ=255)
                                                                                           gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCgat                                                                                                                                                 <  1:60428/62‑1 (MQ=255)
                                                                                           gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCgat                                                                                                                                                 >  2:60428/1‑62 (MQ=255)
                                                                                                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGa                                                     >  3:1876806/1‑149 (MQ=255)
                                                                                                    cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAg                                                 >  4:923486/1‑149 (MQ=255)
                                                                                                     gTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGc                                                >  2:157491/1‑149 (MQ=255)
                                                                                                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTc                                                            >  4:106926/1‑135 (MQ=255)
                                                                                                            tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCgg                                         <  1:1420787/149‑1 (MQ=255)
                                                                                                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                                                                           >  4:1712495/1‑96 (MQ=255)
                                                                                                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                                                                                           <  3:1712495/96‑1 (MQ=255)
                                                                                                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaa                                                                                       >  2:230653/1‑100 (MQ=255)
                                                                                                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaa                                                                                       <  1:230653/100‑1 (MQ=255)
                                                                                                                tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCct                                     <  4:1876806/149‑1 (MQ=255)
                                                                                                                 acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAcccc                                                              <  2:286954/123‑1 (MQ=255)
                                                                                                                   aaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc                                                                             <  1:287011/106‑1 (MQ=255)
                                                                                                                   aaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCggg                                                         >  4:1534933/1‑126 (MQ=255)
                                                                                                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAtt                                >  1:512336/1‑149 (MQ=255)
                                                                                                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAtt                                >  3:1277588/1‑149 (MQ=255)
                                                                                                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAtt                                >  1:972397/1‑149 (MQ=255)
                                                                                                                            aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCagag                         >  3:256915/1‑149 (MQ=255)
                                                                                                                               tGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGcc                      <  4:1463282/149‑1 (MQ=255)
                                                                                                                                 gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAt                                 <  3:1514567/136‑1 (MQ=255)
                                                                                                                                     ccTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAttt                               >  1:1249590/1‑134 (MQ=255)
                                                                                                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCgg                                                          >  4:633992/1‑106 (MQ=255)
                                                                                                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagga               <  4:256915/149‑1 (MQ=255)
                                                                                                                                       ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat              >  2:14831/1‑149 (MQ=255)
                                                                                                                                       ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGc                                                >  3:966449/1‑115 (MQ=255)
                                                                                                                                       ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGc                                                <  4:966449/115‑1 (MQ=255)
                                                                                                                                        tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata             <  4:1277588/149‑1 (MQ=255)
                                                                                                                                        tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata             <  2:512336/149‑1 (MQ=255)
                                                                                                                                        aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                  <  1:1141369/125‑1 (MQ=255)
                                                                                                                                        aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                  >  2:1141369/4‑128 (MQ=255)
                                                                                                                                        aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                  <  1:1741025/125‑1 (MQ=255)
                                                                                                                                        aggatcatcGCGGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                  >  2:1741025/4‑128 (MQ=255)
                                                                                                                                            tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggataggat         >  1:1066089/1‑149 (MQ=255)
                                                                                                                                               tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                                                                                                  <  1:1744154/57‑1 (MQ=255)
                                                                                                                                               tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                                                                                                  >  2:1744154/1‑57 (MQ=255)
                                                                                                                                               tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGc                                                <  3:145593/107‑1 (MQ=255)
                                                                                                                                                  ccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaa                    >  4:1231229/1‑132 (MQ=255)
                                                                                                                                                  ccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaa                    <  3:1231229/132‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCG  >  NC_000913/3007221‑3007516

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