Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,885 G→T 6.2% Q36K (CAG→AAG)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,8850GT 6.2% 882.1 / 20.8 241Q36K (CAG→AAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (106/120);  new base T (9/6);  total (115/126)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 241 total aligned reads displayed.

CTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGC  >  CP009273/3958739‑3959030
                                                                                                                                                  |                                                                                                                                                   
 tGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGaa                                                                                                                                                   <  2:2971990/148‑1 (MQ=255)
    ttCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGACTgg                                                                                                                                                >  2:2863797/1‑148 (MQ=255)
    ttCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATgg                                                                                                                                                >  2:2316405/1‑148 (MQ=255)
     tCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATggg                                                                                                                                               <  2:2056512/148‑1 (MQ=255)
       aCCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAACCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCg                                                                                                                                             <  2:707342/148‑1 (MQ=255)
       aCCTGGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGGTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCg                                                                                                                                             <  2:2646660/148‑1 (MQ=255)
         cgtGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa                                                                                                                                           <  2:1855252/146‑1 (MQ=255)
         cTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa                                                                                                                                           >  1:850338/1‑148 (MQ=255)
           tGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACa                                                                                                                                         >  1:459444/1‑148 (MQ=255)
            gCGATCGGCAAGCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCTGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAg                                                                                                                                        >  1:1043916/1‑148 (MQ=255)
             cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGt                                                                                                                                       <  1:3105711/148‑1 (MQ=255)
               aTCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTa                                                                                                                                     >  1:1455846/1‑148 (MQ=255)
               aTCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTa                                                                                                                                     <  1:1851837/148‑1 (MQ=255)
                  ggCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGa                                                                                                                                  <  2:3008980/148‑1 (MQ=255)
                   gCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGaa                                                                                                                                 <  2:769379/148‑1 (MQ=255)
                      aTCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAcc                                                                                                                              >  1:2738763/1‑148 (MQ=255)
                        cGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑t                                                                                                                          >  2:3552734/1‑148 (MQ=255)
                        cGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑t                                                                                                                          >  1:2953677/1‑148 (MQ=255)
                            ggATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttt                                                                                                                      <  1:313936/148‑1 (MQ=255)
                             gATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑tttttt                                                                                                                     <  2:1875668/148‑1 (MQ=255)
                              aTGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttttt                                                                                                                    >  1:109125/1‑148 (MQ=255)
                                gagagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttttttt                                                                                                                  >  1:563916/1‑147 (MQ=255)
                                gagagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttttttt                                                                                                                  >  1:3531325/1‑147 (MQ=255)
                                 agagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCa                                                                                                                 >  1:1590459/1‑148 (MQ=255)
                                  gagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑tttttttttca                                                                                                                >  1:1010552/1‑145 (MQ=255)
                                  gagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCaa                                                                                                                >  1:283806/1‑148 (MQ=255)
                                   agagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCTGAAGGCCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑tttttttttcaa                                                                                                               <  2:1050982/148‑5 (MQ=255)
                                      gaagaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCTGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGc                                                                                                            <  2:779319/148‑1 (MQ=255)
                                         gaTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTc                                                                                                         >  2:2461210/1‑148 (MQ=255)
                                          aTAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCt                                                                                                        >  1:2997662/1‑148 (MQ=255)
                                           tAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCt                                                                                                        >  2:2846142/1‑148 (MQ=255)
                                              tGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACaa                                                                                                    >  2:2216602/1‑148 (MQ=255)
                                                aaaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAc                                                                                                   <  1:2495170/148‑1 (MQ=255)
                                                aaaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAAcc                                                                                                  >  2:3253796/1‑148 (MQ=255)
                                                 aaCGTTGCCGTAAGAAACGCCATCGTTTTCCCGGGACCGGGTTGCGCCGGCCCGGCTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTTCAAGGGCTTCTACAAc                                                                                                   <  1:1283730/148‑1 (MQ=255)
                                                 aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAAcc                                                                                                  >  2:3322786/1‑148 (MQ=255)
                                                 aaCGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCt                                                                                                 >  1:960925/1‑148 (MQ=255)
                                                   cGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCtt                                                                                                <  2:3102643/148‑1 (MQ=255)
                                                        aCGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                                                                          <  2:2968144/148‑1 (MQ=255)
                                                           tAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCGGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCa                                                                                        <  2:591774/148‑1 (MQ=255)
                                                            aaGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAg                                                                                       >  2:2531348/1‑148 (MQ=255)
                                                            aaGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAg                                                                                       >  2:1501585/1‑148 (MQ=255)
                                                             aGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                                                                     >  1:2702521/1‑148 (MQ=255)
                                                              gAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAggg                                                                                     >  1:1181562/1‑148 (MQ=255)
                                                                aaCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCg                                                                                   >  1:1373338/1‑148 (MQ=255)
                                                                 aCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa                                                                                  >  1:2018985/1‑148 (MQ=255)
                                                                 aCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                                                                 >  2:862051/1‑148 (MQ=255)
                                                                    ccATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                                                                               <  1:1084530/148‑1 (MQ=255)
                                                                    ccATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc                                                                              >  1:2025171/1‑148 (MQ=255)
                                                                      aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                                                                            >  1:3356863/1‑148 (MQ=255)
                                                                      aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                                                                            >  1:3356904/1‑148 (MQ=255)
                                                                        cGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGa                                                                           >  1:2461695/1‑148 (MQ=255)
                                                                          ttttCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                        >  2:420252/1‑148 (MQ=255)
                                                                           tttCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                        <  1:838582/148‑1 (MQ=255)
                                                                           tttCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaaaa                                                                        <  2:459444/148‑1 (MQ=255)
                                                                              cccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt                                                                     >  2:788402/1‑148 (MQ=255)
                                                                               ccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                                                                    <  2:160973/148‑1 (MQ=255)
                                                                                cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg                                                                  <  1:714928/148‑1 (MQ=255)
                                                                                 ggTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt                                                                 <  1:3394570/148‑1 (MQ=255)
                                                                                  gTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGACACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttgtgt                                                                 >  2:1914801/1‑148 (MQ=255)
                                                                                  gTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgat                                                                >  2:1919075/1‑146 (MQ=255)
                                                                                   tACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt                                                                >  1:1602348/1‑148 (MQ=255)
                                                                                    aCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct                                                              <  1:774494/148‑1 (MQ=255)
                                                                                     ccGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct                                                              >  2:3084809/1‑148 (MQ=255)
                                                                                      cGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg                                                             >  1:2148267/1‑148 (MQ=255)
                                                                                          ttGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaa                                                         >  1:2305733/1‑148 (MQ=255)
                                                                                           tGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa                                                        <  2:3174701/148‑1 (MQ=255)
                                                                                               ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgtg                                                    >  2:716958/1‑148 (MQ=255)
                                                                                                cTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt                                                   <  1:1960808/148‑1 (MQ=255)
                                                                                                 tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                  >  2:1210833/1‑148 (MQ=255)
                                                                                                 tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt                                                 <  2:1150746/148‑1 (MQ=255)
                                                                                                 ggCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt                                                  <  2:1788650/147‑1 (MQ=255)
                                                                                                    ccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTg                                               <  2:563916/148‑1 (MQ=255)
                                                                                                        cTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCa                                           >  1:3176624/1‑148 (MQ=255)
                                                                                                          aCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                        >  1:3118459/1‑148 (MQ=255)
                                                                                                           cGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                        <  2:850338/148‑1 (MQ=255)
                                                                                                           cGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAg                                        <  2:2119057/148‑1 (MQ=255)
                                                                                                             tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgtg                                      >  1:3129213/1‑148 (MQ=255)
                                                                                                                cGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                  <  1:2612984/148‑1 (MQ=255)
                                                                                                                cGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                  <  1:382206/148‑1 (MQ=255)
                                                                                                                 gACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                  <  1:2947048/148‑1 (MQ=255)
                                                                                                                 gACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt                                  <  1:2680639/148‑1 (MQ=255)
                                                                                                                    ccGCCAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGaa                              >  2:2932917/1‑148 (MQ=255)
                                                                                                                    ccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGa                               <  2:1266806/148‑1 (MQ=255)
                                                                                                                     cGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAat                             <  2:337887/148‑1 (MQ=255)
                                                                                                                       ccAGCGTCAGCGGAAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAatat                           <  2:2702521/148‑1 (MQ=255)
                                                                                                                         aGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATc                          >  2:1277619/1‑148 (MQ=255)
                                                                                                                           cGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGc                       >  2:1749750/1‑148 (MQ=255)
                                                                                                                               aGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTa                    <  1:1709361/148‑1 (MQ=255)
                                                                                                                                gCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                  <  1:1464904/148‑1 (MQ=255)
                                                                                                                                   gtaGGGCCAGTGACTGAATGGGCGTACAGTTATTAAACCCTTTTTTTTTCAAGGGATTCTACAACCTTCTGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGATCATAGTGTGGTAGAATATCAGCTTACTAt                <  1:51079/146‑1 (MQ=255)
                                                                                                                                   gAAGGGCCAGTGCCTGAATGGTCGTACAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTAt                <  2:2948551/148‑1 (MQ=255)
                                                                                                                                     agtgCCAGTGCATGAATTGTCGTACATTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTg              <  2:2503399/145‑1 (MQ=255)
                                                                                                                                     aGGGCCTGTGCCTGAATGGGCGTAAAGTTATGAAACCCTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTg              <  2:2013920/148‑1 (MQ=255)
                                                                                                                                      gggCCAGTGCCTTAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCt            >  2:1471489/1‑148 (MQ=255)
                                                                                                                                       ggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCtt           >  2:2416736/1‑148 (MQ=255)
                                                                                                                                        gCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCttt          >  1:234977/1‑148 (MQ=255)
                                                                                                                                          cAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTAc        >  2:1391013/1‑148 (MQ=255)
                                                                                                                                              gCCTGAATGGGCTTACAGTTATGAAACCCTTTTTTTTTCAAGGGATTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGaa     <  2:3183846/148‑1 (MQ=255)
                                                                                                                                                cTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGc  <  1:1281486/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGC  >  CP009273/3958739‑3959030

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: