Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,768,689 C→A 100% S215* (TCA→TAA)  ydiB → quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,768,6890CA99.2% 823.3 / ‑2.2 255S215* (TCA→TAA) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (115/138);  total (116/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 258 total aligned reads displayed.

TTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGG  >  CP009273/1768544‑1768833
                                                                                                                                                 |                                                                                                                                                
  cgcgCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTg                                                                                                                                              <  1:740536/148‑1 (MQ=255)
   gcgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTgg                                                                                                                                             <  1:2745338/148‑1 (MQ=255)
    cgcAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACAGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGt                                                                                                                                            <  1:3254260/148‑1 (MQ=255)
      cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTa                                                                                                                                          <  1:355037/148‑1 (MQ=255)
      cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTa                                                                                                                                          >  1:2932712/1‑148 (MQ=255)
      cAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTa                                                                                                                                          >  2:1248418/1‑148 (MQ=255)
        gCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                        >  2:3421612/1‑148 (MQ=255)
        gCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                        >  1:1738262/1‑148 (MQ=255)
        gCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAAt                                                                                                                                        >  1:1511057/1‑148 (MQ=255)
         cGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATg                                                                                                                                       >  1:1603842/1‑148 (MQ=255)
            gTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGata                                                                                                                                    >  1:882376/1‑148 (MQ=255)
             ttAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGatat                                                                                                                                   <  2:1405156/148‑1 (MQ=255)
              tAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATc                                                                                                                                  >  2:2564050/1‑148 (MQ=255)
               aaTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCa                                                                                                                                 <  1:1396257/148‑1 (MQ=255)
                aTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                <  1:870487/148‑1 (MQ=255)
                aTGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAg                                                                                                                                >  2:846665/1‑148 (MQ=255)
                  gAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTc                                                                                                                              <  1:2526144/148‑1 (MQ=255)
                       cacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTa                                                                                                                         >  1:546994/1‑148 (MQ=255)
                       cacCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTa                                                                                                                         >  1:52054/1‑148 (MQ=255)
                           gATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATc                                                                                                                     >  2:786507/1‑148 (MQ=255)
                              tgtgTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCgg                                                                                                                  >  2:2281713/1‑148 (MQ=255)
                                 gtcgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGAc                                                                                                               <  2:80746/148‑1 (MQ=255)
                                   cgtcACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACtt                                                                                                             >  1:3356879/1‑148 (MQ=255)
                                       aCGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTgg                                                                                                         >  2:2754057/1‑148 (MQ=255)
                                        cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                        <  1:184466/148‑1 (MQ=255)
                                        cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                        >  2:3202521/1‑148 (MQ=255)
                                        cGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGt                                                                                                        <  1:1350725/148‑1 (MQ=255)
                                          gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCa                                                                                                      <  1:3052475/148‑1 (MQ=255)
                                          gTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCa                                                                                                      <  2:1158292/148‑1 (MQ=255)
                                              ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGa                                                                                                  <  2:169492/148‑1 (MQ=255)
                                              ccGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGa                                                                                                  <  2:169660/148‑1 (MQ=255)
                                               cGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGaa                                                                                                 <  1:677996/148‑1 (MQ=255)
                                                gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAAt                                                                                                <  1:1975808/148‑1 (MQ=255)
                                                gATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAAt                                                                                                <  1:1184674/148‑1 (MQ=255)
                                                  tctcGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                              >  1:1006959/1‑148 (MQ=255)
                                                  tctcGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGc                                                                                              <  1:122637/148‑1 (MQ=255)
                                                      gCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCgtgt                                                                                          >  2:4634/1‑148 (MQ=255)
                                                      gCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCgtgt                                                                                          <  2:1255283/148‑1 (MQ=255)
                                                      gCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCgtgt                                                                                          <  1:1422295/148‑1 (MQ=255)
                                                          aTCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                      >  1:996221/1‑148 (MQ=255)
                                                          aTCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATaa                                                                                      >  1:413950/1‑148 (MQ=255)
                                                                aaGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCa                                                                                >  2:2489037/1‑148 (MQ=255)
                                                                aaGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCa                                                                                >  1:2162446/1‑148 (MQ=255)
                                                                aaGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCa                                                                                <  1:3367134/148‑1 (MQ=255)
                                                                 aGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCat                                                                               <  2:3339961/148‑1 (MQ=255)
                                                                 aGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCat                                                                               <  2:1598155/148‑1 (MQ=255)
                                                                           gAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGt                                                                     >  1:595164/1‑148 (MQ=255)
                                                                            aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                    <  2:2380323/148‑1 (MQ=255)
                                                                            aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                    <  2:1887068/148‑1 (MQ=255)
                                                                            aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                    >  2:272267/1‑148 (MQ=255)
                                                                            aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                    >  1:351105/1‑148 (MQ=255)
                                                                            aaGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGtt                                                                    >  1:3480791/1‑148 (MQ=255)
                                                                              gCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAt                                                                  <  2:3146896/148‑1 (MQ=255)
                                                                               cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                 >  1:1384522/1‑148 (MQ=255)
                                                                               cccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTAtt                                                                 <  2:1981573/148‑1 (MQ=255)
                                                                                ccTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTg                                                                <  2:545440/148‑1 (MQ=255)
                                                                                 cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                               <  1:1130364/148‑1 (MQ=255)
                                                                                 cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                               <  1:2644138/148‑1 (MQ=255)
                                                                                 cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                               >  2:2653056/1‑148 (MQ=255)
                                                                                 cTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgc                                                               <  2:2112827/148‑1 (MQ=255)
                                                                                    gCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTgcagc                                                            >  2:579349/1‑148 (MQ=255)
                                                                                       tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                         >  2:3486369/1‑148 (MQ=255)
                                                                                       tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                         <  1:2387690/148‑1 (MQ=255)
                                                                                       tCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAgg                                                         >  2:3109386/1‑148 (MQ=255)
                                                                                         cgccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgcg                                                       >  1:1811765/1‑148 (MQ=255)
                                                                                         cgccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGgcg                                                       <  2:1738158/148‑1 (MQ=255)
                                                                                           ccgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGca                                                     <  2:1094488/148‑1 (MQ=255)
                                                                                            cgACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAAAAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaa                                                    >  1:1129855/1‑148 (MQ=255)
                                                                                             gACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaac                                                   <  1:2473197/148‑1 (MQ=255)
                                                                                               cATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGcaacaa                                                 >  1:1498089/1‑148 (MQ=255)
                                                                                                aTTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAg                                                <  1:3425724/148‑1 (MQ=255)
                                                                                                  tttAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCt                                              <  1:2698740/148‑1 (MQ=255)
                                                                                                    tAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTgg                                            <  2:52054/148‑1 (MQ=255)
                                                                                                        cAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGc                                        >  2:1567759/1‑148 (MQ=255)
                                                                                                        cAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGc                                        >  2:417704/1‑148 (MQ=255)
                                                                                                           tGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaa                                     <  1:1428386/148‑1 (MQ=255)
                                                                                                            ggCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaaa                                    <  2:1806993/148‑1 (MQ=255)
                                                                                                            ggCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCaaaa                                    <  1:3271310/148‑1 (MQ=255)
                                                                                                                caAAAGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAt                                <  2:725131/148‑1 (MQ=255)
                                                                                                                 aaaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGAtt                               >  1:2060302/1‑148 (MQ=255)
                                                                                                                  aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                              <  1:126672/148‑1 (MQ=255)
                                                                                                                  aaaGTGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTg                              <  1:1923161/148‑1 (MQ=255)
                                                                                                                      tGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATgg                          >  2:2900147/1‑148 (MQ=255)
                                                                                                                      tGGGTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATgg                          <  2:595164/148‑1 (MQ=255)
                                                                                                                        ggTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGAt                        <  2:1333566/148‑1 (MQ=255)
                                                                                                                         gTATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATa                       >  2:2643795/1‑148 (MQ=255)
                                                                                                                          tATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATAc                      <  1:3023257/148‑1 (MQ=255)
                                                                                                                          tATGAAACCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATAc                      <  1:759107/148‑1 (MQ=255)
                                                                                                                               aaCCCCTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAt                 >  2:3041949/1‑148 (MQ=255)
                                                                                                                                  cccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgtt              >  1:2594191/1‑148 (MQ=255)
                                                                                                                                   ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg             <  1:550646/148‑1 (MQ=255)
                                                                                                                                   ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg             >  2:1053710/1‑148 (MQ=255)
                                                                                                                                   ccTTGAGAATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCAtgttg             <  2:882376/148‑1 (MQ=255)
                                                                                                                                        agaATGAATAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTtgtg        >  1:3067479/1‑148 (MQ=255)
                                                                                                                                         gaatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg       >  1:2734441/1‑148 (MQ=255)
                                                                                                                                         gaatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTgg       >  1:979653/1‑148 (MQ=255)
                                                                                                                                          aatgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCACCAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGc      >  1:812887/1‑148 (MQ=255)
                                                                                                                                            tgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa    <  2:3290440/148‑1 (MQ=255)
                                                                                                                                            tgaatAATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCaa    <  2:1511057/148‑1 (MQ=255)
                                                                                                                                              aataatTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAgg  <  2:941277/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTTGAGAATGAATCATTGGTTAATGATATCAGTCTGTTACATCCGGGACTTCTGGTCACTGAATGCGTGTATAACCCGCATATGACGAAGTTATTGCAGCAGGCGCAACAAGCTGGTTGCAAAACGATTGATGGATACGGCATGTTGTTGTGGCAAGG  >  CP009273/1768544‑1768833

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