Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 568,464 G→T 23.6% E30* (GAA→TAA)  ninE → DLP12 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273568,4640GT23.6% 610.3 / 160.1 250E30* (GAA→TAA) ninEDLP12 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (81/110);  new base T (28/31);  total (109/141)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 252 total aligned reads displayed.

TGCCACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACC  >  CP009273/568315‑568612
                                                                                                                                                     |                                                                                                                                                    
  ccACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGaa                                                                                                                                                    >  1:1454658/1‑150 (MQ=255)
   cACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTaaa                                                                                                                                                   >  2:1055009/1‑150 (MQ=255)
   cACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGaaa                                                                                                                                                   <  2:2539295/150‑1 (MQ=255)
    aCCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAAt                                                                                                                                                  <  2:1892935/150‑1 (MQ=255)
       tGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATAcc                                                                                                                                               >  2:473632/1‑150 (MQ=255)
        gACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCa                                                                                                                                              <  2:3339843/150‑1 (MQ=255)
           tATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAaca                                                                                                                                           <  2:98676/150‑1 (MQ=255)
             ttCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACTAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACt                                                                                                                                         <  1:3053839/150‑1 (MQ=255)
             ttCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACt                                                                                                                                         <  1:3737323/150‑1 (MQ=255)
              tCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTg                                                                                                                                        >  1:1761611/1‑150 (MQ=255)
                cacaGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCt                                                                                                                                      <  1:1847543/150‑1 (MQ=255)
                 acaGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTc                                                                                                                                     <  2:3590104/150‑1 (MQ=255)
                 acaGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTc                                                                                                                                     <  2:1397083/150‑1 (MQ=255)
                 acaGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTc                                                                                                                                     <  1:1676385/150‑1 (MQ=255)
                      gCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGata                                                                                                                                <  1:911701/150‑1 (MQ=255)
                      gCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGata                                                                                                                                <  1:3127106/150‑1 (MQ=255)
                       cTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGatat                                                                                                                               <  2:3010765/150‑1 (MQ=255)
                              ttCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCa                                                                                                                        >  2:65118/1‑150 (MQ=255)
                              ttCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCa                                                                                                                        <  2:3815462/150‑1 (MQ=255)
                              ttCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCa                                                                                                                        <  1:2864918/150‑1 (MQ=255)
                                 cAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCt                                                                                                                     <  1:3191975/150‑1 (MQ=255)
                                  aGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCtt                                                                                                                    <  2:2971/150‑1 (MQ=255)
                                  aGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCtt                                                                                                                    <  1:3304171/150‑1 (MQ=255)
                                  aGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCtt                                                                                                                    >  1:2707031/1‑150 (MQ=255)
                                    ctgctgCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTgg                                                                                                                  >  1:3891373/1‑150 (MQ=255)
                                      gctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGtt                                                                                                                <  1:154840/150‑1 (MQ=255)
                                       ctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACCCTGCTCGGATATACCGCCAGCGGGGTTg                                                                                                               >  1:3014496/1‑150 (MQ=255)
                                       ctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTg                                                                                                               <  2:2874472/150‑1 (MQ=255)
                                        tgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGa                                                                                                              >  2:3658444/1‑150 (MQ=255)
                                        tgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGa                                                                                                              <  1:3538437/150‑1 (MQ=255)
                                        tgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGa                                                                                                              <  1:2004716/150‑1 (MQ=255)
                                         gcgcGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGAt                                                                                                             <  1:806188/150‑1 (MQ=255)
                                         gcgcGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGAt                                                                                                             <  2:1315485/150‑1 (MQ=255)
                                            cgACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaa                                                                                                          >  1:1479327/1‑150 (MQ=255)
                                            cgACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaa                                                                                                          <  2:3299910/150‑1 (MQ=255)
                                              aCAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaaaa                                                                                                        >  2:3228797/1‑150 (MQ=255)
                                               cAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaaaaa                                                                                                       <  2:2377213/150‑1 (MQ=255)
                                                aaGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGAGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAAt                                                                                                      <  2:1063549/150‑1 (MQ=255)
                                                 aGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATg                                                                                                     <  1:3625259/150‑1 (MQ=255)
                                                 aGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATg                                                                                                     >  1:407435/1‑150 (MQ=255)
                                                 aGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATg                                                                                                     >  1:2582334/1‑150 (MQ=255)
                                                  gCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATgg                                                                                                    >  1:3654074/1‑150 (MQ=255)
                                                  gCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATgg                                                                                                    >  1:3026427/1‑150 (MQ=255)
                                                     ggaggCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTg                                                                                                 <  1:1342154/150‑1 (MQ=255)
                                                       aggCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCg                                                                                               >  2:4291799/1‑150 (MQ=255)
                                                        ggCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGa                                                                                              <  2:1188472/150‑1 (MQ=255)
                                                          cGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACt                                                                                            >  1:393017/1‑150 (MQ=255)
                                                           gTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTg                                                                                           <  1:326451/150‑1 (MQ=255)
                                                           gTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTg                                                                                           >  1:203340/1‑150 (MQ=255)
                                                               tGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcag                                                                                       <  2:2076355/150‑1 (MQ=255)
                                                               tGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcag                                                                                       >  1:1187130/1‑150 (MQ=255)
                                                                    acacCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagg                                                                                  >  2:3146414/1‑150 (MQ=255)
                                                                     cacCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagga                                                                                 >  2:2059677/1‑150 (MQ=255)
                                                                     cacCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagga                                                                                 <  2:2211817/150‑1 (MQ=255)
                                                                          cTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGt                                                                            >  2:629872/1‑150 (MQ=255)
                                                                          cTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGt                                                                            >  1:4091402/1‑150 (MQ=255)
                                                                           ttATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTc                                                                           <  2:210435/150‑1 (MQ=255)
                                                                              ttCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATg                                                                        <  2:678095/150‑1 (MQ=255)
                                                                               tCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATgg                                                                       >  1:4011865/1‑150 (MQ=255)
                                                                                cGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGACATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTTCGACTGGCAGCAAGGAGGAGTCATGGc                                                                      <  2:3642264/150‑1 (MQ=255)
                                                                                  tgtCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTg                                                                    <  2:1813754/150‑1 (MQ=255)
                                                                                  tgtCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTg                                                                    >  2:2629336/1‑150 (MQ=255)
                                                                                    tCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGAt                                                                  >  1:639090/1‑150 (MQ=255)
                                                                                     cATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGAtt                                                                 >  2:142912/1‑150 (MQ=255)
                                                                                        gAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTga                                                              >  1:771916/1‑150 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                             <  2:1299483/150‑1 (MQ=255)
                                                                                           cGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGaa                                                           <  1:913591/150‑1 (MQ=255)
                                                                                            ggACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTGGAGaaa                                                          >  1:1089890/1‑150 (MQ=255)
                                                                                              acacATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAag                                                        >  1:4076522/1‑150 (MQ=255)
                                                                                               cacaTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagc                                                       >  2:713314/1‑150 (MQ=255)
                                                                                                 caTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagcag                                                     <  2:1838678/150‑1 (MQ=255)
                                                                                                  aTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagcagc                                                    >  1:2905638/1‑150 (MQ=255)
                                                                                                  aTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagcagc                                                    >  1:3483067/1‑150 (MQ=255)
                                                                                                  aTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagcagc                                                    <  1:2360886/150‑1 (MQ=255)
                                                                                                    cTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAgcgc                                                  >  1:4133111/1‑150 (MQ=255)
                                                                                                     tACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGcgcg                                                 >  2:2186400/1‑150 (MQ=255)
                                                                                                      aCAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGt                                                >  2:1533078/1‑150 (MQ=255)
                                                                                                       cAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTg                                               >  2:3720791/1‑150 (MQ=255)
                                                                                                        agagTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTgg                                              >  2:4018732/1‑150 (MQ=255)
                                                                                                        agagTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTgg                                              <  1:699987/150‑1 (MQ=255)
                                                                                                           gTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCg                                           <  2:1557619/150‑1 (MQ=255)
                                                                                                               cAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGaa                                       >  2:505962/1‑150 (MQ=255)
                                                                                                               cAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGaa                                       >  1:1313418/1‑150 (MQ=255)
                                                                                                                aaaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAAt                                      <  2:4132613/150‑1 (MQ=255)
                                                                                                                aaaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAAt                                      <  1:474814/150‑1 (MQ=255)
                                                                                                                aaaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAAt                                      >  1:88889/1‑150 (MQ=255)
                                                                                                                 aaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATg                                     >  1:361084/1‑150 (MQ=255)
                                                                                                                            gCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGaa                          >  1:1830470/1‑150 (MQ=255)
                                                                                                                              gTAAACCTGAGCTGAAGCCATCAGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCTACTGGAAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGAAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                        >  1:4042781/1‑150 (MQ=255)
                                                                                                                               tAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATcc                       <  2:1191707/150‑1 (MQ=255)
                                                                                                                                   ccTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTgg                   <  2:1761611/150‑1 (MQ=255)
                                                                                                                                    cTGAGCTGAAGCCATCCTAAAGCCCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGc                  <  1:1038537/150‑1 (MQ=255)
                                                                                                                                    cTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGc                  <  1:1906567/150‑1 (MQ=255)
                                                                                                                                     tGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCg                 >  2:3534073/1‑150 (MQ=255)
                                                                                                                                     tGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCg                 >  1:1199452/1‑150 (MQ=255)
                                                                                                                                      gAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGt                >  2:2640182/1‑150 (MQ=255)
                                                                                                                                        gCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTAt              >  1:4355689/1‑150 (MQ=255)
                                                                                                                                            aaGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTaa          >  1:2207584/1‑150 (MQ=255)
                                                                                                                                              gCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATg        <  2:2812919/150‑1 (MQ=255)
                                                                                                                                               ccATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATgg       <  2:3995960/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCCACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACC  >  CP009273/568315‑568612

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: