Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,720 G→T 86.4% E295D (GAG→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273337,7200GT86.4% 668.1 / 111.5 285E295D (GAG→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base G (28/11);  new base T (135/109);  total (164/121)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 287 total aligned reads displayed.

CCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337574‑337865
                                                                                                                                                  |                                                                                                                                                 
ccGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaa                                                                                                                                                 <  2:541542/149‑1 (MQ=255)
  gACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaaaa                                                                                                                                               >  1:3243742/1‑149 (MQ=255)
   aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                              >  1:1974121/1‑149 (MQ=255)
    ccAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACg                                                                                                                                             >  2:3212259/1‑149 (MQ=255)
      aGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGcc                                                                                                                                           >  1:2835109/1‑149 (MQ=255)
       gtgtTGTTGGCGCGTTGGAAAAATCCCTCGACACCTTGTTCCAGATTGCGCTGGAATACGACAATGGCGTCGATATTCACCTGCACGAAACCACTCCGGCTGGCGTGGCAGCCATCATTTATATGGTTGAAACGGTAGAGAAAACGCca                                                                                                                                          <  1:2973623/149‑1 (MQ=255)
       gtgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCca                                                                                                                                          >  2:2363099/1‑149 (MQ=255)
        tgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTCGATAAAACGCcac                                                                                                                                         >  2:3914485/1‑149 (MQ=255)
         gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCCGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCcaca                                                                                                                                        >  1:2130379/1‑149 (MQ=255)
         gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCcaca                                                                                                                                        <  1:3021992/149‑1 (MQ=255)
         gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCcaca                                                                                                                                        >  2:1171294/1‑149 (MQ=255)
          ttGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAg                                                                                                                                       >  1:92238/1‑149 (MQ=255)
           tGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGc                                                                                                                                      <  2:3259816/149‑1 (MQ=255)
           tGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGc                                                                                                                                      <  2:4760204/149‑1 (MQ=255)
             aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTg                                                                                                                                    <  1:4687433/149‑1 (MQ=255)
             aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTg                                                                                                                                    >  1:1352562/1‑149 (MQ=255)
             aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTg                                                                                                                                    >  1:3741879/1‑149 (MQ=255)
                     aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCaa                                                                                                                            >  2:234276/1‑149 (MQ=255)
                      tGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                           >  2:1554357/1‑149 (MQ=255)
                      tGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                           >  1:1489575/1‑149 (MQ=255)
                          aaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGa                                                                                                                       >  2:310979/1‑149 (MQ=255)
                           aaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGAc                                                                                                                      <  1:1425721/149‑1 (MQ=255)
                             aTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCa                                                                                                                    >  2:1809851/1‑149 (MQ=255)
                                ccTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCa                                                                                                                 <  2:592659/149‑1 (MQ=255)
                                  tCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                               <  2:1296267/149‑1 (MQ=255)
                                  tCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                               >  1:4635403/1‑149 (MQ=255)
                                   cGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTc                                                                                                              <  1:347077/149‑1 (MQ=255)
                                     acacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCAc                                                                                                            >  1:1534149/1‑149 (MQ=255)
                                     acacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCAc                                                                                                            <  1:1328016/149‑1 (MQ=255)
                                     acacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCAc                                                                                                            >  2:213225/1‑149 (MQ=255)
                                        ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGcc                                                                                                         <  1:2574319/149‑1 (MQ=255)
                                         cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        <  1:1564205/149‑1 (MQ=255)
                                         cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        <  1:2251835/149‑1 (MQ=255)
                                           tGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCttt                                                                                                      >  1:3864897/1‑149 (MQ=255)
                                                cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCt                                                                                                 <  2:4682377/149‑1 (MQ=255)
                                                 aGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTg                                                                                                >  1:93760/1‑149 (MQ=255)
                                                  gATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTgg                                                                                               >  1:1126168/1‑149 (MQ=255)
                                                   atagcgcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGc                                                                                              >  2:3026415/4‑149 (MQ=255)
                                                   aTTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGc                                                                                              <  2:3969581/149‑1 (MQ=255)
                                                     ggcgcGGGACTACGATAAGGGCGTGGATATGCACCTGCACGAAACCACTCCGGCAGGGGTGGCAGCCATCAATTATTTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCaa                                                                                            <  1:3714546/148‑1 (MQ=255)
                                                        gcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGACGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         >  2:1210065/1‑149 (MQ=255)
                                                        gcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         <  2:1405908/149‑1 (MQ=255)
                                                        gcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         >  1:410715/1‑149 (MQ=255)
                                                        gcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         <  2:4539947/149‑1 (MQ=255)
                                                        gcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         <  2:4068502/149‑1 (MQ=255)
                                                          tGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTc                                                                                       >  2:1581709/1‑149 (MQ=255)
                                                             aCTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAAc                                                                                    >  2:370999/1‑149 (MQ=255)
                                                               tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGa                                                                                  <  2:3758717/149‑1 (MQ=255)
                                                                 cGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCACTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGc                                                                                <  2:3685977/149‑1 (MQ=255)
                                                                         gCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTa                                                                        <  2:2666357/149‑1 (MQ=255)
                                                                         gCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTa                                                                        <  2:273295/149‑1 (MQ=255)
                                                                           gTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGa                                                                      <  2:1761002/149‑1 (MQ=255)
                                                                               atatTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGaa                                                                  >  1:4718645/1‑149 (MQ=255)
                                                                                tatTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAAc                                                                 >  1:308477/1‑149 (MQ=255)
                                                                                 atTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACt                                                                >  1:834543/1‑149 (MQ=255)
                                                                                  ttCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTg                                                               >  2:1461102/1‑149 (MQ=255)
                                                                                  ttCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTg                                                               >  2:2768292/1‑149 (MQ=255)
                                                                                          cACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCg                                                       <  1:2191043/149‑1 (MQ=255)
                                                                                             gAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAt                                                    >  2:749174/1‑149 (MQ=255)
                                                                                              aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                   >  2:3017259/1‑149 (MQ=255)
                                                                                               aaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  <  2:1124976/149‑1 (MQ=255)
                                                                                               aaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  <  2:3514075/149‑1 (MQ=255)
                                                                                                   aCTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATggtggt                                              >  2:1945008/1‑149 (MQ=255)
                                                                                                       cGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaa                                          <  2:1530442/149‑1 (MQ=255)
                                                                                                        ggCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaac                                         <  2:1661192/149‑1 (MQ=255)
                                                                                                         gCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaaca                                        >  1:2515151/1‑149 (MQ=255)
                                                                                                          cAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaacaa                                       >  2:4198989/1‑149 (MQ=255)
                                                                                                            ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                     >  1:97566/1‑149 (MQ=255)
                                                                                                            ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                     >  1:4365325/1‑149 (MQ=255)
                                                                                                             gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAtt                                    >  2:2076579/1‑149 (MQ=255)
                                                                                                              cGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAttt                                   <  2:1774523/149‑1 (MQ=255)
                                                                                                              cGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAttt                                   <  2:4784070/149‑1 (MQ=255)
                                                                                                                 ggCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                                >  2:3000278/1‑149 (MQ=255)
                                                                                                                    aGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  1:1248883/1‑149 (MQ=255)
                                                                                                                    aGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  2:2785536/1‑149 (MQ=255)
                                                                                                                    aGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  1:4237521/1‑149 (MQ=255)
                                                                                                                      ccATCAATTATATGGTTGAAACGGTATATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGcc                           <  1:3604706/149‑1 (MQ=255)
                                                                                                                       cATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCACCAAATTTCTATCGCCt                          >  1:3358476/1‑149 (MQ=255)
                                                                                                                        aTCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTc                         >  1:2935815/1‑149 (MQ=255)
                                                                                                                        aTCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTc                         >  2:1384476/1‑149 (MQ=255)
                                                                                                                          cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGa                       <  2:1843768/149‑1 (MQ=255)
                                                                                                                           aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGAc                      <  1:2382177/149‑1 (MQ=255)
                                                                                                                           aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGAc                      >  2:2636425/1‑149 (MQ=255)
                                                                                                                            aTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACg                     >  2:2842799/1‑149 (MQ=255)
                                                                                                                            aTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACg                     >  2:3321372/1‑149 (MQ=255)
                                                                                                                               atatGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTg                  >  1:1623105/1‑149 (MQ=255)
                                                                                                                               atatGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCACCTCGACGGTg                  <  2:803187/149‑1 (MQ=255)
                                                                                                                                tatGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGc                 >  1:18450/1‑149 (MQ=255)
                                                                                                                                   ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGa              <  1:860147/149‑1 (MQ=255)
                                                                                                                                   ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGa              <  2:2815077/149‑1 (MQ=255)
                                                                                                                                   ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGa              <  1:1707863/149‑1 (MQ=255)
                                                                                                                                       gAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg          >  1:2164708/1‑149 (MQ=255)
                                                                                                                                        aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGc         >  1:4507177/1‑149 (MQ=255)
                                                                                                                                         aaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCa        >  2:656018/1‑149 (MQ=255)
                                                                                                                                          aCGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCAc       >  2:1011458/1‑149 (MQ=255)
                                                                                                                                            ggTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGc     <  1:645938/149‑1 (MQ=255)
                                                                                                                                             gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt    <  1:3962538/149‑1 (MQ=255)
                                                                                                                                             gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt    >  2:2181494/1‑149 (MQ=255)
                                                                                                                                             gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt    <  1:201246/149‑1 (MQ=255)
                                                                                                                                               agatAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGc  >  1:442715/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337574‑337865

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: