Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,655 T→A 5.3% intergenic (‑32/+152) yggN ← / ← yggL DUF2884 family putative periplasmic protein/DUF469 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,6550TA 5.3% 327.2 / 7.2 94intergenic (‑32/+152)yggN/yggLDUF2884 family putative periplasmic protein/DUF469 family protein
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base A (2/3);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

ACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGC  >  NC_000913/3101508‑3101802
                                                                                                                                                   |                                                                                                                                                   
aCGGTTTGCGGGCTGACAATCACATCGTCCCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                                                                                                    >  4:589131/1‑149 (MQ=255)
    tttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                                                                                                <  2:394979/149‑1 (MQ=255)
    tttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                                                                                                >  1:482496/1‑149 (MQ=255)
     ttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                                                                                               <  1:734746/149‑1 (MQ=255)
      tGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                                                                                                              <  3:222414/149‑1 (MQ=255)
      tGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                                                                                                              <  3:801784/149‑1 (MQ=255)
        cGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCCTAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcc                                                                                                                                            <  2:13326/149‑1 (MQ=255)
           gCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                                                                                                 >  3:293571/1‑141 (MQ=255)
           gCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                                                                                                 <  4:293571/141‑1 (MQ=255)
           gCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctga                                                                                                                                         <  1:116357/149‑4 (MQ=255)
              gACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                                                                                                      >  1:848778/1‑149 (MQ=255)
              gACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                                                                                                      >  1:835207/1‑149 (MQ=255)
               aCAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                                                                                                                                     >  2:475888/1‑149 (MQ=255)
                 aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                                                                                                   <  3:35787/149‑1 (MQ=255)
                 aaTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                                                                                                   >  1:504667/1‑149 (MQ=255)
                        tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                                                                                                            <  1:492220/149‑1 (MQ=255)
                        tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                                                                                                            <  2:322109/149‑1 (MQ=255)
                        gcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                                                                                                            >  2:558299/2‑149 (MQ=255)
                          gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                                                                                                                          >  2:539464/1‑149 (MQ=255)
                           tCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                                                                                                         <  4:275710/149‑1 (MQ=255)
                             aCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                                                                                                       >  2:567859/1‑149 (MQ=255)
                               gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                                                                                                                     >  4:859460/1‑149 (MQ=255)
                                  gcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                                                                                                  >  3:860734/1‑149 (MQ=255)
                                  gcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                                                                                                  <  4:472078/149‑1 (MQ=255)
                                   cgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtttt                                                                                                                 <  1:333654/149‑1 (MQ=255)
                                    gTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                                                                                                >  3:390341/1‑149 (MQ=255)
                                            gCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCAGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCgg                                                                                                        >  3:230999/1‑149 (MQ=255)
                                             cACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa                                                                                                       <  2:355103/149‑1 (MQ=255)
                                               cTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                                                                                                     >  2:794052/1‑149 (MQ=255)
                                                tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc                                                                                                    >  1:676601/1‑149 (MQ=255)
                                                tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgat                                                                                                                                        >  2:246465/1‑109 (MQ=255)
                                                tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgat                                                                                                                                        <  1:246465/113‑5 (MQ=255)
                                                     gTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGt                                                                                               <  4:604232/149‑1 (MQ=255)
                                                      tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                                                                                              >  4:490973/1‑149 (MQ=255)
                                                      tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                                                                                              >  3:332649/1‑149 (MQ=255)
                                                           gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac                                                                                         <  4:390341/149‑1 (MQ=255)
                                                                cGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCTATGCCGGATGCGGCGTAAACACCTTAt                                                                                    <  3:590751/149‑1 (MQ=255)
                                                                cGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt                                                                                    <  4:87697/149‑1 (MQ=255)
                                                                   tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATAGATAACAAAGTACGAGTTTTCCATGTCGGATGCGGCGTAAACACCATATccc                                                                                 <  4:659988/149‑1 (MQ=255)
                                                                     aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                           >  1:74591/1‑105 (MQ=255)
                                                                     aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                           <  2:74591/105‑1 (MQ=255)
                                                                     aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGc                                                                               >  3:652746/1‑149 (MQ=255)
                                                                          cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                                                                                                >  1:220549/1‑95 (MQ=255)
                                                                          cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                                                                                                <  2:220549/95‑1 (MQ=255)
                                                                             cACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGa                                                                       >  4:661944/1‑149 (MQ=255)
                                                                               cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATc                                                                     >  2:366358/1‑149 (MQ=255)
                                                                                tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                                                                                                  >  2:94782/1‑103 (MQ=255)
                                                                                tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                                                                                                  <  1:94782/103‑1 (MQ=255)
                                                                                tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCg                                                                    >  3:845401/1‑149 (MQ=255)
                                                                                gcaaaaTGCTGCTGGCGGCAGCACTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCg                                                                    >  4:536140/3‑149 (MQ=255)
                                                                                 gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                                   <  4:860734/149‑1 (MQ=255)
                                                                                 gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGt                                                                   >  1:888785/1‑149 (MQ=255)
                                                                                  aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                  <  3:402352/149‑1 (MQ=255)
                                                                                  aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCACTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                  >  4:52102/1‑149 (MQ=255)
                                                                                           ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGTATCTTTcc                                                         <  1:425861/149‑1 (MQ=255)
                                                                                           ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTcc                                                         >  1:670580/1‑149 (MQ=255)
                                                                                            cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCa                                                        <  1:333956/149‑1 (MQ=255)
                                                                                            cgccAGCAGCATTTGGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACCCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGGGCATCTTTCCa                                                        <  4:867206/149‑1 (MQ=255)
                                                                                              gcggcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAgg                                                      >  3:126244/4‑149 (MQ=255)
                                                                                                  cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTg                                                  >  4:374652/1‑149 (MQ=255)
                                                                                                     caTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATa                                               >  2:82365/1‑149 (MQ=255)
                                                                                                       ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAg                                             <  3:517031/149‑1 (MQ=255)
                                                                                                       ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAg                                             >  4:127724/1‑149 (MQ=255)
                                                                                                       ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAg                                             <  2:835207/149‑1 (MQ=255)
                                                                                                              catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAgcg                                      <  1:706031/149‑1 (MQ=255)
                                                                                                                 catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCt                                                   >  4:485830/1‑133 (MQ=255)
                                                                                                                 catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCt                                                   <  3:485830/133‑1 (MQ=255)
                                                                                                                   tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCa                                 <  2:848778/149‑1 (MQ=255)
                                                                                                                      tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGc                              >  1:571105/1‑149 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAg                                        <  1:51390/137‑1 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAg                                        >  2:51390/1‑137 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAg                            >  3:306219/1‑149 (MQ=255)
                                                                                                                        ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAg                            >  4:586570/1‑149 (MQ=255)
                                                                                                                          cccTCACGAAAAAATCAGGCCAGCATTTGCTGGCCCAAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTCTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                          <  3:589131/149‑1 (MQ=255)
                                                                                                                          cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                          <  3:661944/149‑1 (MQ=255)
                                                                                                                          cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                          >  2:267936/1‑149 (MQ=255)
                                                                                                                            cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                        >  2:699401/1‑149 (MQ=255)
                                                                                                                            cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAAAAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                        <  3:490973/149‑1 (MQ=255)
                                                                                                                             tCAAGAAAAAATTAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac                                                                                         >  4:380143/1‑83 (MQ=255)
                                                                                                                             tCAAGAAAAAATTAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac                                                                                         <  3:380143/83‑1 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAttt                  <  4:513284/149‑1 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAttt                  <  2:571105/149‑1 (MQ=255)
                                                                                                                                  aaaaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAttt                  <  3:127724/149‑1 (MQ=255)
                                                                                                                                    aaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGt                <  1:82365/149‑1 (MQ=255)
                                                                                                                                      aaTCAGGCCAGCATTTGCAGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCa              <  4:230999/149‑1 (MQ=255)
                                                                                                                                       aTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCAt             >  1:325759/1‑149 (MQ=255)
                                                                                                                                        tggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATc            >  2:279282/4‑149 (MQ=255)
                                                                                                                                        tCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATc            >  4:353133/1‑149 (MQ=255)
                                                                                                                                        tCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATc            >  4:623496/1‑149 (MQ=255)
                                                                                                                                            gCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTCTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGccc        <  3:682097/149‑1 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt      >  1:284644/1‑149 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt      >  4:519250/1‑149 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt      <  3:859460/149‑1 (MQ=255)
                                                                                                                                              cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt      >  2:890167/1‑149 (MQ=255)
                                                                                                                                                  aTTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGc  <  3:594856/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGC  >  NC_000913/3101508‑3101802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑