Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,720 G→T 100% 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,7200GT99.3% 875.8 / ‑2.0 273E295D (GAG→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (139/132);  total (140/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 273 total aligned reads displayed.

CCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337573‑337865
                                                                                                                                                   |                                                                                                                                                 
cccGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaa                                                                                                                                                 <  1:2555112/150‑1 (MQ=255)
  cGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaaaa                                                                                                                                               >  1:1991934/1‑150 (MQ=255)
    aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACg                                                                                                                                             >  2:4914958/1‑150 (MQ=255)
      cAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGcc                                                                                                                                           <  1:4486535/150‑1 (MQ=255)
         tgttgtTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGGGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCcaca                                                                                                                                        <  1:3459398/150‑1 (MQ=255)
          gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAg                                                                                                                                       <  2:733397/150‑1 (MQ=255)
          gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAg                                                                                                                                       >  2:2684849/1‑150 (MQ=255)
                   gcgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGc                                                                                                                              <  1:387559/150‑1 (MQ=255)
                   gcgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGc                                                                                                                              <  2:1023144/150‑1 (MQ=255)
                    cgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                             >  1:2489920/1‑150 (MQ=255)
                     gATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCaa                                                                                                                            >  1:747900/1‑150 (MQ=255)
                     gATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCaa                                                                                                                            <  1:3694781/150‑1 (MQ=255)
                         gAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTg                                                                                                                        >  2:2386450/1‑150 (MQ=255)
                          aaaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGa                                                                                                                       <  2:5101393/150‑1 (MQ=255)
                           aaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGAc                                                                                                                      >  2:4102503/1‑150 (MQ=255)
                              aTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCAt                                                                                                                   >  2:4628481/1‑150 (MQ=255)
                                  cTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                               >  1:1393457/1‑150 (MQ=255)
                                     gACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCATGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACATCTGAAGGGCAAGCTGACCATCAGTCAc                                                                                                            >  2:623013/1‑150 (MQ=255)
                                         ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  1:3396506/1‑150 (MQ=255)
                                         ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  2:1966746/1‑150 (MQ=255)
                                         ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  2:2735281/1‑150 (MQ=255)
                                         caaTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACTGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACAATCAGTCACGCCt                                                                                                        >  2:1514966/3‑150 (MQ=255)
                                          cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                       <  1:4130709/150‑1 (MQ=255)
                                          cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                       <  1:4892151/150‑1 (MQ=255)
                                           aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCttt                                                                                                      <  2:4004694/150‑1 (MQ=255)
                                           aTGTTCCAGAGTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCttt                                                                                                      <  1:1830884/150‑1 (MQ=255)
                                            tGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTg                                                                                                     >  2:123441/1‑150 (MQ=255)
                                             gTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgc                                                                                                    <  1:2790686/150‑1 (MQ=255)
                                             gTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgc                                                                                                    >  1:4535983/1‑150 (MQ=255)
                                              ttCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgcg                                                                                                   <  1:508372/150‑1 (MQ=255)
                                               tCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgcgc                                                                                                  >  1:4821277/1‑150 (MQ=255)
                                                 cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTg                                                                                                <  2:866846/150‑1 (MQ=255)
                                                 cAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTg                                                                                                >  1:3119715/1‑150 (MQ=255)
                                                    aTTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCa                                                                                             >  1:3964079/1‑150 (MQ=255)
                                                    aTTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCa                                                                                             <  2:3378662/150‑1 (MQ=255)
                                                       gcgcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACg                                                                                          <  1:4300089/150‑1 (MQ=255)
                                                       gcgcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACg                                                                                          >  2:474954/1‑150 (MQ=255)
                                                        cgcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         >  2:1040242/1‑150 (MQ=255)
                                                        cgcTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                         >  1:4605303/1‑150 (MQ=255)
                                                           tGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGGTAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCa                                                                                      >  1:3733062/1‑150 (MQ=255)
                                                           tGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCa                                                                                      >  1:678212/1‑150 (MQ=255)
                                                                 aCGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGc                                                                                >  2:2280719/1‑150 (MQ=255)
                                                                   gACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCaa                                                                              <  2:1416277/150‑1 (MQ=255)
                                                                    aCAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAAc                                                                             <  1:2401285/150‑1 (MQ=255)
                                                                      aaaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                           <  2:286046/150‑1 (MQ=255)
                                                                            gTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAACACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGAt                                                                     <  1:2249763/150‑1 (MQ=255)
                                                                            gTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGAt                                                                     >  1:797308/1‑150 (MQ=255)
                                                                              cGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGa                                                                   >  1:4746123/1‑150 (MQ=255)
                                                                                 tatTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACt                                                                >  2:1563721/1‑150 (MQ=255)
                                                                                  atTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTg                                                               >  1:4944699/1‑150 (MQ=255)
                                                                                       ccTGCACGAACCCACTCCGGCATGCGTTGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGACTGGCTAGCTGAACATCCGTCACGCCTTTGCGCTGGCAAAGCTCAACGAGATACAGGTAGATGAACTGGCGaa                                                          <  2:2403414/150‑1 (MQ=255)
                                                                                        cTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGATAGATAAAACGCCACAGATGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAc                                                         >  1:1877739/1‑150 (MQ=255)
                                                                                         tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAcc                                                        >  1:1606437/1‑150 (MQ=255)
                                                                                         tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAcc                                                        <  2:1984925/150‑1 (MQ=255)
                                                                                         tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAcc                                                        <  1:2303747/150‑1 (MQ=255)
                                                                                           cACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCgg                                                      >  1:2283800/1‑150 (MQ=255)
                                                                                             cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAt                                                    <  2:1991993/150‑1 (MQ=255)
                                                                                             cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAt                                                    <  2:3576444/150‑1 (MQ=255)
                                                                                              gAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                   >  1:4340648/1‑150 (MQ=255)
                                                                                              gAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                   >  1:4248961/1‑150 (MQ=255)
                                                                                               aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGCGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  >  2:2131141/1‑150 (MQ=255)
                                                                                               aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  >  2:3886155/1‑150 (MQ=255)
                                                                                               aaaCCACTCCGACAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  <  2:1444484/150‑1 (MQ=255)
                                                                                                  ccACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggtgg                                               <  1:1029448/150‑1 (MQ=255)
                                                                                                    aCTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGgtggtg                                             >  1:2161469/1‑150 (MQ=255)
                                                                                                     cTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGc                                            <  1:1063298/150‑1 (MQ=255)
                                                                                                       ccGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaa                                          >  2:1698948/1‑150 (MQ=255)
                                                                                                        tggCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaac                                         >  2:4457984/2‑150 (MQ=255)
                                                                                                          gCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaacaa                                       >  2:3474061/1‑150 (MQ=255)
                                                                                                            aGGCGTGGCAGCCCTCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                     >  2:2225216/1‑150 (MQ=255)
                                                                                                            aGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                     <  2:4157407/150‑1 (MQ=255)
                                                                                                            aGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGAGAACCGGATGGTGGTGCAACAAAt                                     >  2:2631884/1‑150 (MQ=255)
                                                                                                             ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAtt                                    <  1:1766658/150‑1 (MQ=255)
                                                                                                             ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAtt                                    <  2:1124339/150‑1 (MQ=255)
                                                                                                                gTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCt                                 <  1:4675502/150‑1 (MQ=255)
                                                                                                                  ggCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                               >  2:4703568/1‑150 (MQ=255)
                                                                                                                   gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTCGATGAACTGGCGAACCGGATGGTGGTGCACCAAATTTCTATc                              >  2:1375485/1‑150 (MQ=255)
                                                                                                                    cAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  1:2705641/1‑150 (MQ=255)
                                                                                                                    cAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  2:2510552/1‑150 (MQ=255)
                                                                                                                    cAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                             >  2:3334643/1‑150 (MQ=255)
                                                                                                                     aGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                            >  2:3912590/1‑150 (MQ=255)
                                                                                                                     aGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                            >  1:203282/1‑150 (MQ=255)
                                                                                                                      gCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGcc                           >  1:1421358/1‑150 (MQ=255)
                                                                                                                       ctataaaTTATTTTGTTGTTAGGGTTGATAATTCGCCAGAGCTGTTGGGGTCGCTGACCCTCAGTCTCGACTTTGAGCTGGCAACGGTCATCGAGCTACTTGTAGTTGTAATGGCGATAGGGATGGTGGTGCAACAAATTTCTTTCGCCt                          <  2:1738063/145‑1 (MQ=255)
                                                                                                                        cATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTc                         <  1:1304850/150‑1 (MQ=255)
                                                                                                                          tCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGa                       <  1:776045/150‑1 (MQ=255)
                                                                                                                           cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGAc                      <  1:2947270/150‑1 (MQ=255)
                                                                                                                            aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACg                     <  1:3635960/150‑1 (MQ=255)
                                                                                                                            aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACg                     >  1:939868/1‑150 (MQ=255)
                                                                                                                              ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGt                   >  2:755879/1‑150 (MQ=255)
                                                                                                                                 tatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGcc                <  2:1609566/150‑1 (MQ=255)
                                                                                                                                    ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAt             >  1:4180237/1‑150 (MQ=255)
                                                                                                                                    ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAt             >  1:4180100/1‑150 (MQ=255)
                                                                                                                                     gTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAtt            <  2:3226522/150‑1 (MQ=255)
                                                                                                                                      ttGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTg           <  1:1677348/150‑1 (MQ=255)
                                                                                                                                       tGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg          >  2:3042838/1‑150 (MQ=255)
                                                                                                                                          aaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCAc       <  1:69493/150‑1 (MQ=255)
                                                                                                                                            cGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGc     <  1:4986603/150‑1 (MQ=255)
                                                                                                                                             ggTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACCAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt    <  1:3166264/150‑1 (MQ=255)
                                                                                                                                               tAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCGGGCAACGCTCAACGAGCAACCGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGc  <  1:2024084/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337573‑337865

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