Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,077,212 A→T 90.6% N326Y (AAT→TAT)  sbcB → exodeoxyribonuclease I; exonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,077,2120AT90.6% 669.5 / 49.6 223N326Y (AAT→TAT) sbcBexodeoxyribonuclease I; exonuclease I
Reads supporting (aligned to +/- strand):  ref base A (9/12);  new base T (96/106);  total (105/118)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 225 total aligned reads displayed.

GACACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTG  >  CP009273/2077065‑2077356
                                                                                                                                                   |                                                                                                                                                
gACACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTa                                                                                                                                                  <  1:2737836/148‑1 (MQ=255)
  cacaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTAt                                                                                                                                                >  2:667072/1‑148 (MQ=255)
   acaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATc                                                                                                                                               >  1:254527/1‑148 (MQ=255)
   acaTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATc                                                                                                                                               >  2:610577/1‑148 (MQ=255)
    caTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCg                                                                                                                                              >  2:1319415/1‑148 (MQ=255)
     aTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                             <  2:1256583/148‑1 (MQ=255)
     aTTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGACCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGt                                                                                                                                             >  2:1108436/1‑148 (MQ=255)
      ttGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTc                                                                                                                                            <  1:307914/148‑1 (MQ=255)
       tGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCa                                                                                                                                           >  2:1079298/1‑148 (MQ=255)
          gcgAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCa                                                                                                                                        <  2:2932726/148‑1 (MQ=255)
           cgAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCAt                                                                                                                                       <  1:1325045/148‑1 (MQ=255)
           cgAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCAt                                                                                                                                       <  1:1814268/148‑1 (MQ=255)
            gAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCAtt                                                                                                                                      <  2:1803305/148‑1 (MQ=255)
             aGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTg                                                                                                                                     >  1:1242792/1‑148 (MQ=255)
              gCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGCTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGc                                                                                                                                    <  2:2522081/148‑1 (MQ=255)
               cGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGcc                                                                                                                                   <  2:18374/148‑1 (MQ=255)
               cGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGcc                                                                                                                                   >  2:917684/1‑148 (MQ=255)
                 tttATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTc                                                                                                                                 >  2:1477756/1‑148 (MQ=255)
                      ataCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGATGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATaa                                                                                                                            <  1:619903/148‑1 (MQ=255)
                       taCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAAc                                                                                                                           <  2:1354473/148‑1 (MQ=255)
                       taCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAAc                                                                                                                           <  2:3125591/148‑1 (MQ=255)
                       taCCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAAc                                                                                                                           <  2:1354453/148‑1 (MQ=255)
                         ccGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCt                                                                                                                         >  1:2317778/1‑148 (MQ=255)
                           gCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGa                                                                                                                       >  1:1065504/1‑148 (MQ=255)
                               aaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                   >  2:2270890/1‑148 (MQ=255)
                               aaaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAAt                                                                                                                   <  2:2060815/148‑1 (MQ=255)
                                aaCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAAtt                                                                                                                  <  1:768620/148‑1 (MQ=255)
                                 aCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTc                                                                                                                 >  2:2393605/1‑148 (MQ=255)
                                 aCCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTc                                                                                                                 >  2:887829/1‑148 (MQ=255)
                                  ccGATCTTGGCGATAACGCCGCCGTTCCGGTTTAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATTCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCt                                                                                                                <  2:445313/148‑1 (MQ=255)
                                  ccGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCt                                                                                                                >  1:2398061/1‑148 (MQ=255)
                                   cGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTg                                                                                                               <  1:2950265/148‑1 (MQ=255)
                                     aTCTTGGCGATAACGCCGCCGTTCCGGGTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATAGTCAGCAATGCCTCGATAACCTTAAAAttattag                                                                                                             >  2:3278206/1‑143 (MQ=255)
                                       cTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTg                                                                                                           >  2:2044225/1‑148 (MQ=255)
                                             gATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATc                                                                                                     <  2:3096706/148‑1 (MQ=255)
                                              aTAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATcc                                                                                                    >  2:1625503/1‑148 (MQ=255)
                                              aTAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATcc                                                                                                    >  2:2828472/1‑148 (MQ=255)
                                               tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  1:686544/1‑148 (MQ=255)
                                               tAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCg                                                                                                   >  2:848651/1‑148 (MQ=255)
                                                 aCGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCa                                                                                                 >  1:487224/1‑148 (MQ=255)
                                                  cgccgcCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCaa                                                                                                <  2:2001073/148‑1 (MQ=255)
                                                  cgccgcCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCaa                                                                                                <  2:203877/148‑1 (MQ=255)
                                                   gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAg                                                                                               <  1:1926419/148‑1 (MQ=255)
                                                   gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAg                                                                                               <  1:1319415/148‑1 (MQ=255)
                                                   gccgccGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAg                                                                                               <  1:1926381/148‑1 (MQ=255)
                                                     cgccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTg                                                                                             <  2:2107556/148‑1 (MQ=255)
                                                       ccgTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTgcg                                                                                           >  1:3090509/1‑148 (MQ=255)
                                                               gTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAg                                                                                   <  2:342715/148‑1 (MQ=255)
                                                                 tAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtg                                                                                 >  1:80658/1‑148 (MQ=255)
                                                                 tAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtg                                                                                 >  2:2173552/1‑148 (MQ=255)
                                                                  aaGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAgtgg                                                                                <  2:2646716/148‑1 (MQ=255)
                                                                       ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           <  2:640/148‑1 (MQ=255)
                                                                       ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           <  2:3265918/148‑1 (MQ=255)
                                                                       ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           <  2:679/148‑1 (MQ=255)
                                                                       ggTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCg                                                                           >  2:2480044/1‑148 (MQ=255)
                                                                        gTGCATATCAATAAAGGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGa                                                                          <  1:3416698/148‑1 (MQ=255)
                                                                           cATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGatat                                                                       <  1:1734846/148‑1 (MQ=255)
                                                                            atatCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATAtt                                                                      >  1:2139794/1‑148 (MQ=255)
                                                                              atCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTcg                                                                    >  2:2418376/1‑148 (MQ=255)
                                                                                  aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                <  2:2648409/148‑1 (MQ=255)
                                                                                  aTAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGa                                                                <  1:2047240/148‑1 (MQ=255)
                                                                                    aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAg                                                              <  1:2746293/148‑1 (MQ=255)
                                                                                    aaaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAg                                                              <  2:228965/148‑1 (MQ=255)
                                                                                     aaTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGc                                                             >  2:3088049/1‑148 (MQ=255)
                                                                                      aTGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGcc                                                            >  1:316522/1‑148 (MQ=255)
                                                                                             gTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGt                                                     <  2:3204913/148‑1 (MQ=255)
                                                                                              tGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGtt                                                    <  1:1367764/148‑1 (MQ=255)
                                                                                              tGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGtt                                                    <  1:981801/148‑1 (MQ=255)
                                                                                               gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   >  2:2452654/1‑148 (MQ=255)
                                                                                               gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   <  1:1597139/148‑1 (MQ=255)
                                                                                               gCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGttt                                                   <  1:2270890/148‑1 (MQ=255)
                                                                                                cTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTa                                                  >  1:1710804/1‑148 (MQ=255)
                                                                                                cTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTa                                                  >  1:2365744/1‑148 (MQ=255)
                                                                                                   gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  1:314551/1‑148 (MQ=255)
                                                                                                   gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  2:1882269/1‑148 (MQ=255)
                                                                                                   gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  2:2781743/1‑148 (MQ=255)
                                                                                                   gCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGc                                               >  1:61434/1‑148 (MQ=255)
                                                                                                    cccAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                              <  1:1611404/148‑1 (MQ=255)
                                                                                                    cccAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                              <  2:1699647/148‑1 (MQ=255)
                                                                                                    cccAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGcc                                              <  1:673633/148‑1 (MQ=255)
                                                                                                       aGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTc                                           <  2:1921792/148‑1 (MQ=255)
                                                                                                       aGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTc                                           <  2:2608426/148‑1 (MQ=255)
                                                                                                                  ctacgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTgg                                <  1:1625503/148‑1 (MQ=255)
                                                                                                                     cgcCCGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATg                             <  1:1129876/148‑1 (MQ=255)
                                                                                                                        ccGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGcac                          <  2:1394007/148‑1 (MQ=255)
                                                                                                                        ccGGAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGcac                          <  1:1435504/148‑1 (MQ=255)
                                                                                                                           gAAGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGc                       >  1:1840286/1‑148 (MQ=255)
                                                                                                                            aaGATGCCGACCGACTTGGAATTTATCGTCAGCATTGCCTAGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTAGCTGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCt                      <  2:3005910/148‑1 (MQ=255)
                                                                                                                            aaGATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCt                      >  1:3499843/1‑148 (MQ=255)
                                                                                                                              gATGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCttt                    <  1:3539616/148‑1 (MQ=255)
                                                                                                                               aTGCCGATCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTta                   <  2:720836/148‑1 (MQ=255)
                                                                                                                               aTGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTta                   <  1:3027874/148‑1 (MQ=255)
                                                                                                                                tGCCGACCGACTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTtat                  <  1:1086927/148‑1 (MQ=255)
                                                                                                                                  ccgaccgaCTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATaa                >  2:1228896/1‑148 (MQ=255)
                                                                                                                                     accgacTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACgg             >  1:2005505/1‑148 (MQ=255)
                                                                                                                                      ccgacTGGGAATTTATCGTCAGCAGTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGc            <  2:1438586/148‑1 (MQ=255)
                                                                                                                                          cTGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCtttt        <  2:1656958/148‑1 (MQ=255)
                                                                                                                                           tGGGAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCttttt       >  1:3341172/1‑148 (MQ=255)
                                                                                                                                            gggAATTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTc      <  1:1101348/148‑1 (MQ=255)
                                                                                                                                                aTTTATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTg  <  1:3470543/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GACACATTGCGCGAGCGTTTATATACCGCAAAAACCGATCTTGGCGATAACGCCGCCGTTCCGGTTAAGCTGGTGCATATCAATAAATGTCCGGTGCTGGCCCAGGCGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTG  >  CP009273/2077065‑2077356

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: