Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,201 T→C 6.3% pseudogene (810/2518 nt) ydbA → putative outer membrane protein, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,2010TC 6.3% 189.7 / 2.6 80pseudogene (810/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (19/52);  new base C (4/2);  total (25/55)
Fisher's exact test for biased strand distribution p-value = 6.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCA  >  NC_000913/1466062‑1466338
                                                                                                                                           |                                                                                                                                          
ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATCCCACCGGTACGCCCCACATCATGAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGAc                                                                                                                                         <  1:878861/143‑1 (MQ=17)
 cacaGTAAACAACAATGGCAAAACCACCGTTGACGAGAAAGAATCCACCGGGACGGAAATCAACGGGAATAACGGGAAAGAGATTCAGGACGGCGATCTAGATGTCAGCGACGGCGGTCACGGTATTGATATCACCGGTGACa                                                                                                                                        >  2:458850/1‑143 (MQ=17)
  acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAg                                                                                                                                       <  2:861938/143‑1 (MQ=255)
   caGTAAACAAAAATGGCATAACCACCGTTGACGGTAAAGAGTACACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgc                                                                                                                                      <  2:123301/143‑1 (MQ=255)
        aacaacAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTAAGCGGAGGCAGTCACGGTATTGAGATAACCGGTGACAGcgcgtc                                                                                                                                 >  1:345377/1‑141 (MQ=14)
           aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTCCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGTATGTCAGCGGCGGCGATCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                                                                              >  1:699625/1‑143 (MQ=255)
           aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTACTAACGGGAAAGTGATTAAGAACAGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                                                                              >  1:244723/1‑143 (MQ=255)
           aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGATGACAGCGCGACGGt                                                                                                                              >  1:97807/1‑143 (MQ=255)
             cagtGGTAATACCACCGTTGATGGTAAAGATTCCCCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                                                                            <  2:181851/140‑1 (MQ=14)
              aactgCAAAACCACCGTTGCCGGTAACGCTTCCACCGGTACGGCAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGCCAGTGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGa                                                                                                                           <  1:272988/139‑1 (MQ=14)
                   cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTACTGATATCACCGGCGACAGCGCGACGGGGGATAACa                                                                                                                      >  2:27396/1‑143 (MQ=255)
                     aaaCCACCGTTGATCGTACATATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                    <  2:302669/143‑1 (MQ=255)
                     aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                    <  1:836137/143‑1 (MQ=255)
                     aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                    <  2:547596/143‑1 (MQ=255)
                     aaaCCACCGTGGACGGTAAAGAGTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                    <  1:203020/143‑1 (MQ=255)
                      aaCCGCCGGTGACGGTAAAGATGCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                                                                   <  1:437818/143‑1 (MQ=255)
                      aaCCACGGTTGACTGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                                                                   <  2:79176/143‑1 (MQ=255)
                      aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                                                                   <  1:950654/143‑1 (MQ=255)
                         caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCa                                                                                                                >  1:68324/1‑143 (MQ=255)
                           ccGTTGACGGTAAAGATTCCACGGTTACGTAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                                                                             <  2:238605/143‑1 (MQ=255)
                             gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCAt                                                                                                           <  1:899348/143‑1 (MQ=255)
                               tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCActttga                                                                                                         <  2:754105/143‑6 (MQ=255)
                                gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAc                                                                                                        <  2:60015/143‑1 (MQ=255)
                                 aCGGTAATGATTCCACGGGTACGGAAGTCAACGGTAATAACGGGAAAGTGATTGAGGACGGCGATCTGGATGTGAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAcc                                                                                                       <  2:269264/143‑1 (MQ=14)
                                 aCGGTAAAGATCCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAcc                                                                                                       <  1:751500/143‑1 (MQ=255)
                                    gTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACAGGGGACAGCGCGACGGCGGATCACAAGGGCACCATGACCGTc                                                                                                    >  2:531391/1‑143 (MQ=255)
                                     tAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACACGGGCACCATGACCGTCa                                                                                                   >  1:203562/1‑143 (MQ=255)
                                     tAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                   >  1:310433/1‑143 (MQ=21)
                                      cagtcttCCACCGGTACGGAAATCAACGGTAAGAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATACCCCCGGTGACAGCGCGATCGAGAAGAACAAGGGCACCATGACCGTCAc                                                                                                  <  2:863882/138‑1 (MQ=14)
                                             ccACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcg                                                                                                                                     <  2:571446/102‑1 (MQ=17)
                                             ccACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcg                                                                                                                                     >  1:571446/1‑102 (MQ=18)
                                               gccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                         <  2:504358/142‑1 (MQ=255)
                                               gccGGTACGCAAATCAACGGTAATAACGGGAAAGTGATGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                         <  2:505045/142‑1 (MQ=255)
                                               aCCGGTACGGAAGGCAACGGTAAGAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACTGTATTGGTATCACCGGTGACAGCGCGGCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                         <  2:118843/143‑1 (MQ=14)
                                               aCCGGTACGGAAATCAACGGTAATAACGGGAATGTGTTTTTTGACGGCGATCTGGATGTCAGCGGCGGCGGCCACGGTAGTGATAGCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                         <  1:468838/143‑1 (MQ=14)
                                                 tggTACGTGAATCAACGGTAGTAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCAGGGTATTGATATCTCCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                       <  2:218442/142‑1 (MQ=14)
                                                     aCGTAAATCGACGGTAGTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAt                                                                                   <  2:767893/143‑1 (MQ=255)
                                                     aCGGAAATCAACGGTAATAACGGGAAGGTGAGTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAt                                                                                   <  1:255177/143‑1 (MQ=255)
                                                     aCGGAAACCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGAATAGAGGGCATCCATGACCGTCACCGATCCGGAGTCCAt                                                                                   <  2:72556/143‑1 (MQ=14)
                                                      cGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATg                                                                                  <  2:370532/143‑1 (MQ=40)
                                                          aaTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTa                                                                              <  2:514510/143‑1 (MQ=255)
                                                           aTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTac                                                                             >  1:357881/1‑142 (MQ=255)
                                                              aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCa                                                                          >  2:917210/1‑143 (MQ=18)
                                                               aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGTTATTGACATCCCCGCGGGCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                                                                         >  1:982742/1‑143 (MQ=11)
                                                                cGGTAACAACGGGAAAGTGATTCCGGACGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGa                                                                        <  1:690558/143‑1 (MQ=255)
                                                                 ggTAATAACGGGAAAGTGATGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGAt                                                                       <  1:256760/143‑1 (MQ=255)
                                                                    gatgactggAAAGTGATTAGGGACGGCGAGCTGGATGTCTGCGGCGGCGGGCACGGTATTGGTGTATGCGTTGACACAACGACGGGCCAGTACAAGGGCGCCACGGCCGTCACCGATCCGGAGTCCATGGGTATCCAGACCGa                                                                    <  2:961207/136‑1 (MQ=9)
                                                                     ataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAggg                                                                                                                  <  1:28680/97‑1 (MQ=255)
                                                                     ataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAggg                                                                                                                  >  2:28680/1‑97 (MQ=255)
                                                                       agcGGGACAGTGAGTCAGGACGGCGCTCGGGATGCCAGGGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATCACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACgg                                                                 <  2:241356/141‑1 (MQ=14)
                                                                       aaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGGCCATGGGTATCCAGATCGACgg                                                                 >  2:769610/1‑143 (MQ=255)
                                                                        ccGGGGAAGTGGTTCAGGACGGCGATCTTGGCGTCAGGCGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATACCAAGCTCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt                                                                <  1:199646/142‑1 (MQ=14)
                                                                          gggAAAGTGATTCAGCACGGCGATCCGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGa                                                              <  1:863991/143‑1 (MQ=255)
                                                                              gaGTGATTCAGGACGGCGAGCGGGTTGTCAGCGGCGGCGGTCACGGTATGGATATCACCGGTGACAGCGCGGCGGTGGATAACAAGGGCACCATGACCGTCGCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAg                                                          <  2:490055/142‑1 (MQ=14)
                                                                              aaGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAg                                                          <  1:443163/143‑1 (MQ=255)
                                                                                    ttCAGGACGGCGATCGGGGTGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc                                                    <  1:771252/143‑1 (MQ=255)
                                                                                    gtCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc                                                    <  1:487034/142‑1 (MQ=255)
                                                                                          aCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAAt                                              <  1:889325/143‑1 (MQ=255)
                                                                                            ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGCTGACAGCGACAAGGCGGTCGTCagcaa                                            >  2:831037/1‑139 (MQ=11)
                                                                                               gATATGGATGGCAGCGGCGGCGGACATGGTATTGATGTCACCGGTGCGAGCGCGACGGTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGa                                         <  2:56432/143‑1 (MQ=14)
                                                                                                gtgtGGATGGCAGCGGCGGCGGTCACGGTATTGGTATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCATATCGACGGTGACAAGGCCATCGTCAATAACGaa                                        <  1:55373/140‑1 (MQ=255)
                                                                                                gtcgggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAGAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGaa                                        <  2:594448/139‑1 (MQ=255)
                                                                                                 tCTGGATGCGAGCGGCGGCGTCCACGGGACATATAGCAGCACTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAg                                       <  1:764689/143‑1 (MQ=14)
                                                                                                  cTGTATGTCCGCGGCGCAGTGCACGGTATTGATATCACCGGTGACAGCGCGACGGTGTATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAgg                                      <  1:494896/143‑1 (MQ=17)
                                                                                                    ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCg                                    >  2:962326/1‑143 (MQ=255)
                                                                                                           aGCGGCGGCGGTCACGGTATTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                             <  1:821369/143‑1 (MQ=255)
                                                                                                            tcggtggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAcc                            <  1:833738/142‑1 (MQ=255)
                                                                                                               gcggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATc                         <  1:654165/143‑1 (MQ=255)
                                                                                                                     gTCACGGTATTGATATCACCGGTGACAGCGCGACGGAGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGGCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCCCCAAc                   >  1:638624/1‑143 (MQ=255)
                                                                                                                        aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                           >  2:498307/1‑84 (MQ=25)
                                                                                                                         cGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                              >  2:770454/1‑64 (MQ=255)
                                                                                                                          ggTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTgg              >  1:117089/1‑143 (MQ=255)
                                                                                                                          ggTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATTGACGGTGACAAGGCGATCGTCAATAACGAAGGCGACAGCACCGCCACCGACGgggg              >  1:227279/1‑143 (MQ=14)
                                                                                                                             aTTGATATCACCGGTGACAGCGCGACGGTGGATGACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac           <  1:658601/143‑1 (MQ=255)
                                                                                                                             aTTGAGATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac           <  1:657480/143‑1 (MQ=255)
                                                                                                                                gATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCTCCGATCCCGACTCCATGGGTATCCAGATCGACGGTCACAAGGCCATCGTCCATAACGCAGGCGAGAGCACCATCACCAACGGTggcaccgg        <  2:697360/143‑1 (MQ=17)
                                                                                                                                gATATCACCGGGGACAGCGCGACGGTGGGTTGCCGCTGCACCAGGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccgg        <  2:433348/143‑1 (MQ=17)
                                                                                                                                      aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCa  >  1:730843/1‑143 (MQ=255)
                                                                                                                                        cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                      <  2:772632/57‑1 (MQ=255)
                                                                                                                                        cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                      >  1:772632/1‑57 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCA  >  NC_000913/1466062‑1466338

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