Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,247 A→T 100% V57E (GTG→GAG)  yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,2470AT98.7% 1040.9 / 2.2 297V57E (GTG→GAG) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (144/149);  total (145/152)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
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 297 total aligned reads displayed.

CGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTAC  >  CP009273/3994101‑3994394
                                                                                                                                                  |                                                                                                                                                   
cGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTccc                                                                                                                                                  >  1:1166219/1‑150 (MQ=255)
cGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTccc                                                                                                                                                  >  2:4101757/1‑150 (MQ=255)
cGATGTACAAACCGAGCGGCGCATAGCTTTCATCAGAGATGCCGCCGGTCATCGCCCTCTGCACCTGACCCTGACTATGAAGATCCGCCAGCAGAGCATGTAACGTTCCCCCCGCTTCGCCAGTGGCTTTTTGCCAGTACATTTTCAccc                                                                                                                                                  <  2:138146/150‑1 (MQ=255)
   tGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCAGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGct                                                                                                                                               <  2:1434218/150‑1 (MQ=255)
     tACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGctct                                                                                                                                             >  1:1355542/1‑150 (MQ=255)
     tACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGctct                                                                                                                                             >  1:1499100/1‑150 (MQ=255)
       cAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGt                                                                                                                                           <  1:390164/150‑1 (MQ=255)
        aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                          >  1:2189280/1‑150 (MQ=255)
        aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGtt                                                                                                                                          <  1:379716/150‑1 (MQ=255)
            cGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAg                                                                                                                                      <  2:1595206/150‑1 (MQ=255)
              aGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGa                                                                                                                                    >  1:1273955/1‑150 (MQ=255)
                cggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAtt                                                                                                                                  <  2:3424454/150‑1 (MQ=255)
                cggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAtt                                                                                                                                  <  1:188343/150‑1 (MQ=255)
                   cgcATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGATTGAGGATTAAc                                                                                                                               <  2:1968085/150‑1 (MQ=255)
                       tAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTGTGAGGATTAACGGta                                                                                                                           >  2:922071/1‑150 (MQ=255)
                        aGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtct                                                                                                                          >  2:1708906/1‑148 (MQ=255)
                         gCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGtata                                                                                                                         <  1:1769/150‑1 (MQ=255)
                          cTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCACCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATAc                                                                                                                        >  1:686920/1‑150 (MQ=255)
                             tcatcaTAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGcc                                                                                                                     <  1:946648/150‑1 (MQ=255)
                              catcatAGATGCCTCCGTTCATCGCCATCTGAACCTGACCCTTACTATTAATATCCGCCAGCAGAGCATGTAAAGTTCCCCACGCTTCGCCATTGGCTTTTTGCCATTACATTTTCTCCAGCTCTGTTTGAGGATTAACGGTATACGCCt                                                                                                                    <  2:665778/150‑1 (MQ=255)
                              catcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCt                                                                                                                    <  1:4369678/150‑1 (MQ=255)
                              catcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCt                                                                                                                    <  1:3919789/150‑1 (MQ=255)
                               atcatAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCGTTTTGCCAGTACCTTTTCTCCCGCTCTGTTGGAGGATTAACGGTATACTCCTg                                                                                                                   >  2:712118/1‑150 (MQ=255)
                                 catAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                                 >  1:3619015/1‑150 (MQ=255)
                                  atAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTAc                                                                                                                >  1:1630378/1‑150 (MQ=255)
                                   tAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACg                                                                                                               >  2:2486179/1‑150 (MQ=255)
                                      aTGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTc                                                                                                            >  2:1881264/1‑150 (MQ=255)
                                       tGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCa                                                                                                           <  1:1087332/150‑1 (MQ=255)
                                           gccgTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGt                                                                                                       >  2:2293890/1‑150 (MQ=255)
                                            ccgTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTc                                                                                                      <  2:2396590/150‑1 (MQ=255)
                                                 cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATc                                                                                                 >  1:2747103/1‑150 (MQ=255)
                                                 cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATc                                                                                                 <  2:116505/150‑1 (MQ=255)
                                                 cATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATc                                                                                                 <  1:1972969/150‑1 (MQ=255)
                                                   tCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTg                                                                                               <  1:967715/150‑1 (MQ=255)
                                                    cGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTga                                                                                              >  1:3681572/1‑150 (MQ=255)
                                                       cATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGagag                                                                                           >  1:3001303/1‑150 (MQ=255)
                                                          cTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgc                                                                                        >  2:379349/1‑150 (MQ=255)
                                                          cTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgc                                                                                        >  2:710904/1‑150 (MQ=255)
                                                           tGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcg                                                                                       <  2:806646/150‑1 (MQ=255)
                                                            gCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATATGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcgc                                                                                      <  1:3740337/150‑1 (MQ=255)
                                                            gCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcgc                                                                                      >  1:1887592/1‑150 (MQ=255)
                                                             cACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCa                                                                                     >  1:120791/1‑150 (MQ=255)
                                                                ctgaccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatc                                                                                  <  2:882745/150‑1 (MQ=255)
                                                                   accctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcatc                                                                               >  2:2600291/1‑150 (MQ=255)
                                                                    ccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAAtcatca                                                                              >  1:3541402/1‑150 (MQ=255)
                                                                    ccctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAAtcatca                                                                              >  1:903294/1‑150 (MQ=255)
                                                                         acTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGc                                                                         <  2:2166438/150‑1 (MQ=255)
                                                                         acTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGc                                                                         >  1:2458840/1‑150 (MQ=255)
                                                                           tATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCaa                                                                       <  2:4633347/150‑1 (MQ=255)
                                                                            aTTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAAc                                                                      <  2:2758671/150‑1 (MQ=255)
                                                                             ttAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACa                                                                     >  2:4189585/1‑150 (MQ=255)
                                                                             ttAAGATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACa                                                                     <  2:4353281/150‑1 (MQ=255)
                                                                              tAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTTGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAg                                                                    <  2:2857971/150‑1 (MQ=255)
                                                                              tAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAg                                                                    <  1:3558543/150‑1 (MQ=255)
                                                                              tAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAg                                                                    <  2:1692443/150‑1 (MQ=255)
                                                                                 tatCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCaa                                                                 >  1:4451542/1‑150 (MQ=255)
                                                                                   tCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAAc                                                               <  1:1575682/150‑1 (MQ=255)
                                                                                    ccGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACa                                                              <  1:2253259/150‑1 (MQ=255)
                                                                                     cGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACaa                                                             >  2:3289351/1‑150 (MQ=255)
                                                                                      gCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAg                                                            <  1:1011564/150‑1 (MQ=255)
                                                                                      gCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAg                                                            <  1:3293774/150‑1 (MQ=255)
                                                                                           cagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAg                                                       <  2:3896338/150‑1 (MQ=255)
                                                                                            agagCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                      >  1:2340467/1‑150 (MQ=255)
                                                                                               gCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAg                                                   >  2:4257438/1‑150 (MQ=255)
                                                                                                     aaCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAgg                                             <  2:3910387/150‑1 (MQ=255)
                                                                                                     aaCGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAgg                                             <  2:2844197/150‑1 (MQ=255)
                                                                                                        gTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                          >  1:2586278/1‑150 (MQ=255)
                                                                                                            cccACGCTTCGCCATTGGATTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaa                                      <  1:2028583/150‑1 (MQ=255)
                                                                                                             ccACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaa                                     >  1:3397603/1‑150 (MQ=255)
                                                                                                              cACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGaaaa                                    <  2:4533277/150‑1 (MQ=255)
                                                                                                                 gCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATc                                 >  1:3054956/1‑150 (MQ=255)
                                                                                                                  cTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATcc                                <  1:3966236/150‑1 (MQ=255)
                                                                                                                   ttCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                               <  1:4260464/150‑1 (MQ=255)
                                                                                                                   ttCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                               <  2:3643148/150‑1 (MQ=255)
                                                                                                                     cGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGtt                             <  1:2667954/150‑1 (MQ=255)
                                                                                                                     cGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAATCCGttt                            >  2:2718965/1‑150 (MQ=255)
                                                                                                                      gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                            <  2:3215724/150‑1 (MQ=255)
                                                                                                                        cATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTga                          >  1:3997478/1‑150 (MQ=255)
                                                                                                                            ggCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAtt                      <  2:1977927/150‑1 (MQ=255)
                                                                                                                              cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                    <  2:3255341/150‑1 (MQ=255)
                                                                                                                               tttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCaa                   >  1:1037183/1‑150 (MQ=255)
                                                                                                                               tttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTGTGAGATTCaa                   >  1:1680053/1‑150 (MQ=255)
                                                                                                                                ttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                  <  1:1875012/150‑1 (MQ=255)
                                                                                                                                  ttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGt                <  2:3552358/150‑1 (MQ=255)
                                                                                                                                  ttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGt                <  1:319228/150‑1 (MQ=255)
                                                                                                                                   tGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg               >  2:2928398/1‑150 (MQ=255)
                                                                                                                                   tGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg               <  2:2087449/150‑1 (MQ=255)
                                                                                                                                    gCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGa              <  2:2973704/150‑1 (MQ=255)
                                                                                                                                    gCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGa              >  1:4629735/1‑150 (MQ=255)
                                                                                                                                    gCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGa              >  2:1077313/1‑150 (MQ=255)
                                                                                                                                     ccAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGAt             <  2:3069801/150‑1 (MQ=255)
                                                                                                                                      cAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATc            <  2:1269784/150‑1 (MQ=255)
                                                                                                                                       aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa           >  2:608372/1‑150 (MQ=255)
                                                                                                                                       aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa           >  2:2323230/1‑150 (MQ=255)
                                                                                                                                       aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa           >  2:2503376/1‑150 (MQ=255)
                                                                                                                                          aCATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTc        >  1:551080/1‑150 (MQ=255)
                                                                                                                                           cATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc       <  2:2194370/150‑1 (MQ=255)
                                                                                                                                            aTTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt      <  1:2849752/150‑1 (MQ=255)
                                                                                                                                             ttttCTCCCGCTCTGTTTGAGGAGTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCtt     >  1:4044515/1‑150 (MQ=255)
                                                                                                                                               ttCTCCCGCTCTGTTTGAGGATTAACGGTATAGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTAc  >  2:4440101/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTAC  >  CP009273/3994101‑3994394

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: