Predicted mutation
evidence seq id position mutation annotation gene description
RA CP043523 1,885,813 2 bp→CC coding (39‑40/351 nt) FTN78_chr1931 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0435231,885,8130TC99.6% 499.9 / NA 266I13I (ATT→ATCFTN78_chr1931hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (133/132);  minor base A (0/1);  total (133/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*CP0435231,885,8140AC99.2% 503.6 / NA 266R14R (AGG→CGG) FTN78_chr1931hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (133/131);  total (133/133)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

AAGTAGTTTATCCGAGGGATCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGAAAAAGATGTCAGTTGACCGTTCATATAGTCCCTACGAGGTTATTAGGATGTATCACGATAGAGGAATGATGAAATGGGGGGCATTTGCGACTGGGGAATTAACTGAAGCCCAAAATACTTTTGAAAAAGAAAAAAAAGAAGATAAGATAATTCAAGCATTAC  >  CP043523/1885689‑1885931
                                                                                                                             ||                                                                                                                      
  gtagtTTATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGatga                                                                                               >  2:4354096/1‑148 (MQ=14)
     gtTTATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGAtgatga                                                                                            >  1:4094771/1‑148 (MQ=9)
      tttATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGaa                                                                                           <  1:2160498/148‑1 (MQ=9)
      tttATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGaa                                                                                           <  1:3016658/148‑1 (MQ=9)
      tttATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGaa                                                                                           <  2:2239899/148‑1 (MQ=14)
      tttATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGaa                                                                                           >  1:530145/1‑148 (MQ=14)
       ttATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGaaa                                                                                          <  1:1696591/148‑1 (MQ=9)
        tATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAAt                                                                                         >  1:112008/1‑148 (MQ=9)
        tATCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAAt                                                                                         >  1:2548492/1‑148 (MQ=9)
         aTCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATg                                                                                        <  1:3465910/148‑1 (MQ=9)
         aTCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATg                                                                                        <  1:2633136/148‑1 (MQ=9)
         aTCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATg                                                                                        >  2:3957817/1‑148 (MQ=14)
          tCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATgg                                                                                       >  2:411461/1‑148 (MQ=9)
          tCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATgg                                                                                       <  1:1212382/148‑1 (MQ=9)
          tCCGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATgg                                                                                       >  1:2920411/1‑148 (MQ=9)
            cgaggattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATgggg                                                                                     <  1:449889/141‑1 (MQ=9)
            cGAGGATTCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATgggg                                                                                     <  1:206999/148‑1 (MQ=14)
             gaggattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATggggc                                                                                    >  2:1231262/7‑147 (MQ=9)
             gaggattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATggggc                                                                                    >  1:2313519/7‑147 (MQ=14)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 <  1:1062907/145‑1 (MQ=9)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 <  1:4459021/145‑1 (MQ=9)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 >  2:243568/4‑148 (MQ=14)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 <  1:1576746/145‑1 (MQ=9)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 >  1:1504845/4‑148 (MQ=14)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 >  2:2743521/4‑148 (MQ=14)
                gattCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCa                                                                                 >  1:1035787/4‑148 (MQ=14)
                 attCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCAt                                                                                >  2:1154676/3‑148 (MQ=14)
                  ttCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCAtt                                                                               <  1:4353908/147‑1 (MQ=14)
                  ttCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCAtt                                                                               <  1:1761469/147‑1 (MQ=14)
                  ttCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCAtt                                                                               >  1:3158717/2‑148 (MQ=14)
                  ttCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCAtt                                                                               >  1:4345688/2‑148 (MQ=14)
                    cACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTg                                                                             <  1:752226/148‑1 (MQ=14)
                     aCGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGc                                                                            >  2:3479348/1‑148 (MQ=14)
                      cGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCg                                                                           >  2:4221398/1‑148 (MQ=14)
                       gtgtCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGa                                                                          >  2:1415097/1‑148 (MQ=14)
                       gtgtCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGa                                                                          <  2:4062288/148‑1 (MQ=14)
                        tgtCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGAc                                                                         >  2:2986703/1‑148 (MQ=14)
                           cATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTgg                                                                      >  1:964130/1‑148 (MQ=14)
                            aTTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTgga                                                                     <  2:1714070/148‑2 (MQ=14)
                            aTTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTgga                                                                     <  1:4044225/148‑2 (MQ=14)
                            aTTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAATCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTTCGACTgga                                                                     >  2:808044/1‑147 (MQ=14)
                             tttATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTggag                                                                    <  1:2322851/148‑3 (MQ=9)
                             tttATCGGAACAAACTCATTGGAGGACATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTggag                                                                    >  2:3702487/1‑146 (MQ=255)
                                 tCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaact                                                                >  1:495335/1‑146 (MQ=14)
                                  cGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactc                                                               <  2:4275548/148‑4 (MQ=9)
                                  cGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactc                                                               >  2:3189738/1‑145 (MQ=14)
                                   ggAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactca                                                              >  2:2870570/1‑144 (MQ=14)
                                     aaCAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactcacg                                                            <  1:35277/148‑7 (MQ=9)
                                     aaCAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactcacg                                                            <  2:3553693/148‑7 (MQ=9)
                                     aaCAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGaactcacg                                                            >  2:3764657/1‑142 (MQ=9)
                                                    ggTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACt                                             >  1:1363456/1‑148 (MQ=2)
                                                    ggTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACt                                             <  2:1993585/148‑1 (MQ=255)
                                                    ggTCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACt                                             <  1:4088061/148‑1 (MQ=255)
                                                      tCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACttt                                           <  1:4681264/148‑1 (MQ=2)
                                                      tCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACttt                                           >  1:243556/1‑148 (MQ=2)
                                                      tCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACttt                                           <  1:3719186/148‑1 (MQ=2)
                                                      tCATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACttt                                           >  2:2107114/1‑148 (MQ=2)
                                                       cATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACtttt                                          >  1:1121449/1‑148 (MQ=2)
                                                       cATGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACtttt                                          >  1:3131464/1‑148 (MQ=9)
                                                        aTGCAGCAAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGGCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTg                                         <  1:3350338/148‑1 (MQ=255)
                                                          gcagcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaa                                       <  2:961564/148‑1 (MQ=9)
                                                          gcagcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaa                                       >  2:4071376/1‑148 (MQ=9)
                                                          gcagcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaa                                       >  1:3268670/1‑148 (MQ=9)
                                                          gcagcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaa                                       <  2:3030589/148‑1 (MQ=9)
                                                           cagcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaa                                      <  1:3052059/148‑1 (MQ=9)
                                                            agcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaa                                     <  1:1631205/148‑1 (MQ=255)
                                                            agcaAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaa                                     >  1:4285257/1‑148 (MQ=9)
                                                              caAGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaag                                   >  1:707391/1‑148 (MQ=9)
                                                               aaGTCAAGACAAGGAGTTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaaga                                  <  2:1526055/148‑1 (MQ=9)
                                                                aGTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaa                                 >  2:857055/1‑148 (MQ=9)
                                                                 gTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaa                                <  1:1417187/148‑1 (MQ=9)
                                                                 gTCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaa                                <  1:1314753/148‑1 (MQ=9)
                                                                  tCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTTCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaaa                               >  1:700408/1‑148 (MQ=9)
                                                                  tCAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaaa                               <  1:387898/148‑1 (MQ=9)
                                                                   cAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCATGCTTATAACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaaaaaaaa                              <  2:2577637/148‑1 (MQ=9)
                                                                   cAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaaaa                              <  2:2049074/148‑1 (MQ=9)
                                                                   cAAGACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTgaaaaagaaaaa                              >  1:3617556/1‑148 (MQ=9)
                                                                      gACAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGaaaaaaaa                           >  2:4666855/1‑148 (MQ=9)
                                                                       aCAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAaag                          >  2:1881852/1‑148 (MQ=9)
                                                                        cAAGGAGGTAAAAGATGTCAGTTGACCGTTCATACAGGCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGTATGTTGAAATGGGGCGCAGTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAaaga                         <  1:4157984/148‑1 (MQ=255)
                                                                           ggaggtAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAagatga                      <  1:1884460/148‑1 (MQ=12)
                                                                            gaggtAAAAGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAgatgat                     <  1:3827059/148‑1 (MQ=255)
                                                                              ggtaaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataa                   <  1:2410604/145‑1 (MQ=9)
                                                                              ggtaaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataa                   >  2:2323879/4‑148 (MQ=9)
                                                                               gtaaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataag                  >  1:2293771/3‑148 (MQ=9)
                                                                               gtaaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataag                  >  1:839984/3‑148 (MQ=9)
                                                                                taaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataaga                 >  2:1468276/2‑148 (MQ=9)
                                                                                taaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAAAAAAAAAAGatgataaga                 <  2:1372082/147‑1 (MQ=255)
                                                                                taaaaGATGTCAGTTGACCGTTAATACAGCCAATATGAAGTAATACGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAAAAAAAAAAGatgataaga                 <  2:618916/147‑1 (MQ=255)
                                                                                taaaaGATGGCAGTTGACCCTTCATACAGCCCCTATGAAGACATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTTGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataaga                 <  2:2794461/147‑1 (MQ=255)
                                                                                 aaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataagaa                <  1:1981747/148‑2 (MQ=9)
                                                                                 aaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataagaa                <  1:2648924/148‑2 (MQ=9)
                                                                                 aaaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGACATCCGTGCTAATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAAAAAAAAAAGatgataagaa                <  2:4549686/148‑2 (MQ=9)
                                                                                  aaaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGatgataagaaa               <  2:1031804/148‑3 (MQ=9)
                                                                                   aaGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATgataagaaaa              <  1:2859944/148‑1 (MQ=9)
                                                                                    aGATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATGATAAGAAAAt             >  1:2720011/1‑148 (MQ=9)
                                                                                     gATGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATGATAAGAAAAtt            >  1:959565/1‑148 (MQ=255)
                                                                                      aTGTCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATGATAAGAAAATTc           >  2:461221/1‑148 (MQ=255)
                                                                                         tCAGTTGACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATGATAAGAAAATTCAAg        <  2:1772092/148‑1 (MQ=255)
                                                                                               gACCGTTCATACAGCCCCTATGAAGTCATCCGTGCTTATCACGATAGAGGGATGATGAAATGGGGCGCATTTGCGACTGGAGAACTCACGGAGGCTCAAAATACTTTTGAAAAAGAAAAAAAAGATGATAAGAAAATTCAAGCTTTAc  <  2:2389290/148‑1 (MQ=255)
                                                                                                                             ||                                                                                                                      
AAGTAGTTTATCCGAGGGATCACGTGTCATTTATCGGAACAAACTCATTGGAGGTCATGCAGCAAGTCAAGACAAGGAGGAAAAAGATGTCAGTTGACCGTTCATATAGTCCCTACGAGGTTATTAGGATGTATCACGATAGAGGAATGATGAAATGGGGGGCATTTGCGACTGGGGAATTAACTGAAGCCCAAAATACTTTTGAAAAAGAAAAAAAAGAAGATAAGATAATTCAAGCATTAC  >  CP043523/1885689‑1885931

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