Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 523,413 G→T 19.4% pseudogene (6/708 nt) rhsH → pseudogene, Rhs familly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273523,4130GT19.4% 915.2 / 156.9 324pseudogene (6/708 nt)rhsHpseudogene, Rhs familly protein
Reads supporting (aligned to +/- strand):  ref base G (160/101);  new base T (28/35);  total (188/136)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

AAAACCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACC  >  CP009273/523264‑523561
                                                                                                                                                     |                                                                                                                                                    
  aaCCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGctct                                                                                                                                                    >  1:3958064/1‑150 (MQ=255)
     cATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTgat                                                                                                                                                 <  2:4899810/150‑1 (MQ=255)
      atatTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTgatg                                                                                                                                                >  2:1143640/1‑150 (MQ=255)
      atatTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTgatg                                                                                                                                                <  2:4018776/150‑1 (MQ=255)
      atatTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTgatg                                                                                                                                                >  2:4280556/1‑150 (MQ=255)
      atatTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGAGGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTgatg                                                                                                                                                >  2:3567464/1‑150 (MQ=255)
          ttAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTgatggatg                                                                                                                                            >  1:3771218/1‑150 (MQ=255)
           tAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGc                                                                                                                                           >  1:2742323/1‑150 (MQ=255)
            aaTTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCg                                                                                                                                          >  1:4680522/1‑150 (MQ=255)
             aTTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCgg                                                                                                                                         >  1:2179946/1‑150 (MQ=255)
               tACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGAt                                                                                                                                       >  2:523403/1‑150 (MQ=255)
                 cGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATgg                                                                                                                                     >  1:1584288/1‑150 (MQ=255)
                  ggAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGa                                                                                                                                    >  1:592776/1‑150 (MQ=255)
                       tACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACAt                                                                                                                               >  1:2966712/1‑150 (MQ=255)
                         ccAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTg                                                                                                                             <  2:2328492/150‑1 (MQ=255)
                               aGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGa                                                                                                                       >  1:2155057/1‑150 (MQ=255)
                                gTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAg                                                                                                                      >  1:3480583/1‑150 (MQ=255)
                                  tGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCg                                                                                                                    <  1:4403347/150‑1 (MQ=255)
                                           gAGGGAGCACAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTatg                                                                                                           <  2:3840406/150‑1 (MQ=255)
                                             gggAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGTATGCGGATGGAAACATTGAGTGGAGCGGGGAGTatgat                                                                                                         >  2:2454178/1‑150 (MQ=255)
                                              ggAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTatgatg                                                                                                        >  1:4887796/1‑150 (MQ=255)
                                               gAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTAtgatga                                                                                                       >  2:3674584/1‑150 (MQ=255)
                                                aGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAg                                                                                                      <  1:1272125/150‑1 (MQ=255)
                                                 gCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGt                                                                                                     >  1:5096165/1‑150 (MQ=255)
                                                  cAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTg                                                                                                    >  1:2044826/1‑150 (MQ=255)
                                                  cAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTg                                                                                                    >  1:1773982/1‑150 (MQ=255)
                                                   aaaGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTgg                                                                                                   >  1:3580945/1‑150 (MQ=255)
                                                    aaGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCATACATAGTAAATAATTCTGGCGCTGATGTATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTggg                                                                                                  >  1:3385926/1‑150 (MQ=255)
                                                     aGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTgggg                                                                                                 <  2:4957251/150‑1 (MQ=255)
                                                         aCCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAAc                                                                                             <  2:2694895/150‑1 (MQ=255)
                                                           cTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCa                                                                                           >  1:1196921/1‑150 (MQ=255)
                                                              aTGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCt                                                                                        <  2:3207493/150‑1 (MQ=255)
                                                              aTGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCt                                                                                        >  2:4589194/1‑150 (MQ=255)
                                                               tGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtg                                                                                       >  2:1665178/1‑150 (MQ=255)
                                                               tGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtg                                                                                       <  1:4544280/150‑1 (MQ=255)
                                                                  aGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaat                                                                                    <  1:288799/150‑1 (MQ=255)
                                                                  aGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaat                                                                                    >  2:1734226/1‑150 (MQ=255)
                                                                  aGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaat                                                                                    <  2:2581645/150‑1 (MQ=255)
                                                                   gTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGAGGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaatg                                                                                   >  2:620450/1‑150 (MQ=255)
                                                                        gggTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAagag                                                                              >  2:3492488/1‑150 (MQ=255)
                                                                        gggTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGAAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGAATGATGAGTGGGGCAACCAGCAGAATGAagag                                                                              >  2:3713273/1‑150 (MQ=255)
                                                                         ggtggtGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAgaga                                                                             >  1:1561494/1‑150 (MQ=255)
                                                                           tggtgTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAAc                                                                           >  1:2836056/1‑150 (MQ=255)
                                                                            ggtgTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAAcc                                                                          >  2:263117/1‑150 (MQ=255)
                                                                             gtgtATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAAccc                                                                         <  1:3610693/150‑1 (MQ=255)
                                                                              tgtATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCg                                                                        <  2:3342440/150‑1 (MQ=255)
                                                                                 atttattATAGAAATTAATAAGAAAAATGGAGGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCAt                                                                     <  2:3057802/150‑1 (MQ=255)
                                                                                     attatAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCAcc                                                                 <  1:1610239/150‑1 (MQ=255)
                                                                                      ttatAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCt                                                                >  2:912962/1‑150 (MQ=255)
                                                                                         taGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCa                                                             >  1:703071/1‑150 (MQ=255)
                                                                                            aaaTTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCa                                                          >  1:2752832/1‑150 (MQ=255)
                                                                                              aTTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGc                                                        >  2:127215/1‑150 (MQ=255)
                                                                                                 aataaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGt                                                     >  2:4671696/1‑150 (MQ=255)
                                                                                                 aataaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGt                                                     >  2:238874/1‑150 (MQ=255)
                                                                                                  ataaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTa                                                    >  1:758018/1‑150 (MQ=255)
                                                                                                  ataaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTa                                                    <  2:4149780/150‑1 (MQ=255)
                                                                                                   taaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAc                                                   >  2:772344/1‑150 (MQ=255)
                                                                                                    aaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAcc                                                  <  1:4602848/150‑1 (MQ=255)
                                                                                                    aaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGACCCCGCATCACCTGCACCAGCCGTAcc                                                  >  1:2246515/1‑150 (MQ=255)
                                                                                                    aaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAcc                                                  >  1:407351/1‑150 (MQ=255)
                                                                                                     aGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                 <  2:1907620/150‑1 (MQ=255)
                                                                                                     aGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                 <  1:863149/150‑1 (MQ=255)
                                                                                                     aGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                 <  1:4808676/150‑1 (MQ=255)
                                                                                                      gAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCgg                                                >  1:3055073/1‑150 (MQ=255)
                                                                                                      gAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCgg                                                >  1:5003727/1‑150 (MQ=255)
                                                                                                        aaaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCt                                              >  2:1401978/1‑150 (MQ=255)
                                                                                                        aaaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCt                                              <  2:2408267/150‑1 (MQ=255)
                                                                                                        aaaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCt                                              <  2:2502116/150‑1 (MQ=255)
                                                                                                         aaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTg                                             <  2:568873/150‑1 (MQ=255)
                                                                                                           aTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGcc                                           <  1:1354842/150‑1 (MQ=255)
                                                                                                              gATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCggg                                        >  2:3913955/1‑150 (MQ=255)
                                                                                                               aTGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCgggg                                       >  1:1797986/1‑150 (MQ=255)
                                                                                                               aTGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCgggg                                       <  2:2069630/150‑1 (MQ=255)
                                                                                                                tgtgTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgc                                      <  1:1950995/150‑1 (MQ=255)
                                                                                                                tgtgTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgc                                      <  1:2344003/150‑1 (MQ=255)
                                                                                                                  tgtTTTGAATTTCCTACCTAGTAAATAATGCTGGCGCTGCTGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcag                                    <  2:1453475/150‑1 (MQ=255)
                                                                                                                  tgtTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcag                                    <  1:1558513/150‑1 (MQ=255)
                                                                                                                  tgtTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcag                                    <  1:4718306/150‑1 (MQ=255)
                                                                                                                     tttGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGcagcag                                 <  2:2436906/150‑1 (MQ=255)
                                                                                                                        gAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTAt                              >  1:1373294/1‑150 (MQ=255)
                                                                                                                         aaTTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATg                             <  2:1984239/150‑1 (MQ=255)
                                                                                                                         aaTTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATg                             <  1:2639913/150‑1 (MQ=255)
                                                                                                                         aaTTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATg                             >  1:2570617/1‑150 (MQ=255)
                                                                                                                             tCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATaa                         >  2:4644373/1‑150 (MQ=255)
                                                                                                                             tCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATaa                         >  1:978637/1‑150 (MQ=255)
                                                                                                                               cTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAgg                       >  1:3175130/1‑150 (MQ=255)
                                                                                                                                tACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGa                      >  2:2144586/1‑150 (MQ=255)
                                                                                                                                tACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGa                      >  1:4113093/1‑150 (MQ=255)
                                                                                                                                 aCATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAg                     <  1:4236987/150‑1 (MQ=255)
                                                                                                                                 aCATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAg                     >  2:2239519/1‑150 (MQ=255)
                                                                                                                                   aTAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTc                   >  2:85537/1‑150 (MQ=255)
                                                                                                                                     aGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCgg                 >  2:4981078/1‑150 (MQ=255)
                                                                                                                                        aaaTAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGc              >  2:1836617/1‑150 (MQ=255)
                                                                                                                                         aataatGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCt             >  2:3899655/1‑150 (MQ=255)
                                                                                                                                          ataatGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTg            >  2:2785645/1‑150 (MQ=255)
                                                                                                                                             atGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtac         >  1:1997020/1‑150 (MQ=255)
                                                                                                                                             atGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtac         >  1:3345706/1‑150 (MQ=255)
                                                                                                                                              tGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtact        <  2:3596782/150‑1 (MQ=255)
                                                                                                                                               gcgggCGCTGATGGATGCGGCTGGACACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtacta       <  1:3864619/147‑1 (MQ=255)
                                                                                                                                                cTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtactac      >  2:3574762/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACC  >  CP009273/523264‑523561

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: