Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,564,939 G→T 20.5% pseudogene (1204/2937 nt) yjiV → pseudogene; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,564,9390GT20.5% 874.3 / 142.9 304pseudogene (1204/2937 nt)yjiVpseudogene; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (119/122);  new base T (34/29);  total (153/151)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

AGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGAC  >  CP009273/4564791‑4565088
                                                                                                                                                    |                                                                                                                                                     
 gCCGCGACTATCCTGGTATCCCGAACCTGGACGGATATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGaa                                                                                                                                                     <  1:812774/150‑1 (MQ=255)
   cgcgACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAAt                                                                                                                                                   >  1:1460064/1‑150 (MQ=255)
     cgACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAAtgt                                                                                                                                                 <  2:2704804/150‑1 (MQ=255)
      gACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAAtgtt                                                                                                                                                >  1:3743232/1‑150 (MQ=255)
      gACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAAtgtt                                                                                                                                                <  2:2218262/150‑1 (MQ=255)
         tATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAAtgttgtt                                                                                                                                             >  1:2038303/1‑150 (MQ=255)
         tATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCAATATCTTAAGCCTGAAGTGTTTAATGAAAACGTCGGCCACCTTAAGCTTTCTGAAAtgttgtt                                                                                                                                             >  1:3335168/1‑150 (MQ=255)
            ccTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATgttgttgtt                                                                                                                                          <  2:2539365/150‑1 (MQ=255)
            ccTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATgttgttgtt                                                                                                                                          >  1:491709/1‑150 (MQ=255)
               ggTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAAGGTTGTTGTTCCg                                                                                                                                       >  2:172067/1‑150 (MQ=255)
                 tATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCggg                                                                                                                                     <  2:3936117/150‑1 (MQ=255)
                 tATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCggg                                                                                                                                     <  1:3569566/150‑1 (MQ=255)
                 tATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCggg                                                                                                                                     >  2:1298973/1‑150 (MQ=255)
                  aTCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGa                                                                                                                                    >  1:2953414/1‑150 (MQ=255)
                   tCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGaa                                                                                                                                   <  1:454085/150‑1 (MQ=255)
                     ccGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAAtt                                                                                                                                 <  2:2177646/150‑1 (MQ=255)
                       gAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCa                                                                                                                               >  2:3407072/1‑150 (MQ=255)
                            tGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTTAGCTTTCTGTATTGTTGTTGTTCCGGGAATTCATGCGt                                                                                                                          >  2:3855686/1‑150 (MQ=255)
                             ggACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTc                                                                                                                         >  2:1244398/1‑150 (MQ=255)
                             ggACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTc                                                                                                                         >  1:585007/1‑150 (MQ=255)
                             ggACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTc                                                                                                                         >  2:4497006/1‑150 (MQ=255)
                             ggACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTc                                                                                                                         >  2:4039912/1‑150 (MQ=255)
                                  gagaTGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCt                                                                                                                    <  2:2867322/150‑1 (MQ=255)
                                  gagaTGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCt                                                                                                                    >  2:247591/1‑150 (MQ=255)
                                    taTGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGAAATATAACCTTTATCTTATGCCTGATGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTaa                                                                                                                  <  1:2460848/149‑1 (MQ=255)
                                      tGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAAc                                                                                                                <  1:593686/150‑1 (MQ=255)
                                      tGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAAc                                                                                                                <  1:3909686/150‑1 (MQ=255)
                                           aaCGATCTTAAAGAGAGTGCAGATATCCGCGGTCGGGTTCTGCATTATAAAAATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCGCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTAATGCGTCGCCCTAAACGGACt                                                                                                           <  2:2344158/150‑1 (MQ=255)
                                            aCGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTa                                                                                                          >  2:4025222/1‑150 (MQ=255)
                                             cGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTaa                                                                                                         <  2:1944395/150‑1 (MQ=255)
                                              gAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAAc                                                                                                        >  2:3554583/1‑150 (MQ=255)
                                                gCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAg                                                                                                      >  2:4197311/1‑150 (MQ=255)
                                                   tAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCt                                                                                                   <  1:4109418/150‑1 (MQ=255)
                                                   tAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCt                                                                                                   <  1:1867968/150‑1 (MQ=255)
                                                     aaGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTgg                                                                                                 >  1:3622798/1‑150 (MQ=255)
                                                       gagagagCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGaa                                                                                               >  1:316976/1‑150 (MQ=255)
                                                        agagagCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGaaa                                                                                              <  2:3360635/150‑1 (MQ=255)
                                                         gagagaCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGATGTTGTTCCGGGAATTCATGCGTCGCCCTAAACTGACTAACAGCCTGGAAac                                                                                             >  1:981316/1‑150 (MQ=255)
                                                           gagCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAacac                                                                                           <  1:173082/150‑1 (MQ=255)
                                                             gCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAg                                                                                         <  2:920245/150‑1 (MQ=255)
                                                              ccGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAgg                                                                                        <  1:1443908/150‑1 (MQ=255)
                                                              ccGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGGTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAgg                                                                                        >  2:1525282/1‑150 (MQ=255)
                                                                gATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAgggg                                                                                      >  1:2488412/1‑150 (MQ=255)
                                                                   atCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTg                                                                                   <  2:3042144/150‑1 (MQ=255)
                                                                      cgcgGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGtt                                                                                >  1:2362555/1‑150 (MQ=255)
                                                                            ccGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggtt                                                                          <  2:242204/150‑1 (MQ=255)
                                                                              ggTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttgg                                                                        >  1:4397780/1‑150 (MQ=255)
                                                                               gTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggt                                                                       >  1:3041196/1‑150 (MQ=255)
                                                                                ttCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggtt                                                                      <  2:4593022/150‑1 (MQ=255)
                                                                                   tGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTAcc                                                                   >  1:420964/1‑150 (MQ=255)
                                                                                   tGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTAcc                                                                   >  1:4633717/1‑150 (MQ=255)
                                                                                   tGCAATATAACAATTATATTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTAcc                                                                   <  1:2864379/150‑1 (MQ=255)
                                                                                     cAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAg                                                                 >  2:3276096/1‑150 (MQ=255)
                                                                                      aaTATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                                                <  1:2425423/150‑1 (MQ=255)
                                                                                       atatAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAggg                                                               <  2:3724293/150‑1 (MQ=255)
                                                                                       atatAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAggg                                                               >  2:262015/1‑150 (MQ=255)
                                                                                           aaCCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTg                                                           >  2:2115752/1‑150 (MQ=255)
                                                                                             ccATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGga                                                         <  1:3961289/150‑1 (MQ=255)
                                                                                                 tATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGaaa                                                     >  2:4280554/1‑150 (MQ=255)
                                                                                                   tCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAAt                                                   <  2:4151793/150‑1 (MQ=255)
                                                                                                    cTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATa                                                  <  2:754969/150‑1 (MQ=255)
                                                                                                     ttATGCCTGAAGTGTTTAATGAATACGGCGGCCCCATTAAGCTTTCTGAAATGTTGTTGTTGCGTGAATTCATGCGTCGCACTAACCGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGTCTGGAGTAAATAc                                                 <  2:31854/150‑1 (MQ=255)
                                                                                                       aaGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACAt                                               <  2:57087/150‑1 (MQ=255)
                                                                                                       aaGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACAt                                               <  1:3813183/150‑1 (MQ=255)
                                                                                                       aaGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACAt                                               >  1:3265578/1‑150 (MQ=255)
                                                                                                        aGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATa                                              <  2:3492422/150‑1 (MQ=255)
                                                                                                           cTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAaga                                           >  1:202447/1‑150 (MQ=255)
                                                                                                           cTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAaga                                           >  2:649044/1‑150 (MQ=255)
                                                                                                              aaGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGcc                                        <  1:3521623/150‑1 (MQ=255)
                                                                                                               aGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGccc                                       >  1:2104899/1‑150 (MQ=255)
                                                                                                                  gtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCtt                                    <  1:1044104/150‑1 (MQ=255)
                                                                                                                  gtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCtt                                    <  1:2334/150‑1 (MQ=255)
                                                                                                                  gtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCtt                                    >  2:4000602/1‑150 (MQ=255)
                                                                                                                  gtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCtt                                    <  2:2750598/150‑1 (MQ=255)
                                                                                                                  gtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCtt                                    <  1:447270/150‑1 (MQ=255)
                                                                                                                      aaTGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCAc                                <  2:3453439/150‑1 (MQ=255)
                                                                                                                       aTGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACt                               <  1:4057380/150‑1 (MQ=255)
                                                                                                                        tGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTg                              <  2:3017804/150‑1 (MQ=255)
                                                                                                                         gAAAACGGCGGCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTgg                             >  2:3662915/1‑150 (MQ=255)
                                                                                                                            aaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAg                          <  1:1133965/150‑1 (MQ=255)
                                                                                                                             aCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAgg                         >  2:2238340/1‑150 (MQ=255)
                                                                                                                              cggcggCCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGa                        >  2:4499702/1‑150 (MQ=255)
                                                                                                                               ggcggcCCCCTTAAGCTTTCTTAAATGTTGTTGTTCCTGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGaa                       <  1:1518043/150‑1 (MQ=255)
                                                                                                                                  ggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGaaaaa                    >  2:766301/1‑150 (MQ=255)
                                                                                                                                    cccccTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAgg                  >  2:1135897/1‑150 (MQ=255)
                                                                                                                                    cccccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAgg                  >  1:3359948/1‑150 (MQ=255)
                                                                                                                                    cccccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAgg                  >  1:688885/1‑150 (MQ=255)
                                                                                                                                     ccccTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGa                 >  1:314418/1‑150 (MQ=255)
                                                                                                                                       ccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGAtt               >  1:52142/1‑150 (MQ=255)
                                                                                                                                       ccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGAtt               >  1:3684646/1‑150 (MQ=255)
                                                                                                                                        cTTAAGCTTTCTTAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTa              >  1:3445276/1‑150 (MQ=255)
                                                                                                                                             gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCTCTGGCAGGAAAAAGGATTAACGCt         <  2:2042114/150‑1 (MQ=255)
                                                                                                                                             gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCt         <  2:188349/150‑1 (MQ=255)
                                                                                                                                             gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCt         <  2:3235783/150‑1 (MQ=255)
                                                                                                                                              cTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTg        >  2:2029739/1‑150 (MQ=255)
                                                                                                                                               tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTgg       >  1:790392/1‑150 (MQ=255)
                                                                                                                                               tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGACAGGAAAAAGGATTAACGCTgg       <  2:2874127/150‑1 (MQ=255)
                                                                                                                                                   tGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGa   <  2:456196/150‑1 (MQ=255)
                                                                                                                                                   tGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGa   <  2:1612849/150‑1 (MQ=255)
                                                                                                                                                   tGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGa   >  2:3777583/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCCGCGACTATCCTGGTATCCCGAACCTGGACGGAGATGGTTAACGAGCTTAAAGAGAGAGCCGATATCCGCGGGCCGGTTCTGCAATATAACCATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGAC  >  CP009273/4564791‑4565088

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: