Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,959,004 G→A 100% I131I (ATC→ATT)  lgt ← phosphatidylglycerol‑prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,959,0040GA99.5% 745.9 / ‑4.4 211I131I (ATC→ATT) lgtphosphatidylglycerol‑prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (117/93);  total (117/94)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-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 211 total aligned reads displayed.

GTTGGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCG  >  CP009273/2958857‑2959149
                                                                                                                                                   |                                                                                                                                                 
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 >  2:744683/1‑150 (MQ=255)
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 <  1:594618/150‑1 (MQ=255)
gttggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAAt                                                                                                                                                 >  2:3037365/1‑150 (MQ=255)
  tggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaa                                                                                                                                               >  1:4057059/1‑150 (MQ=255)
  tggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaa                                                                                                                                               <  2:367427/150‑1 (MQ=255)
   ggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATaaa                                                                                                                                              >  2:820060/1‑150 (MQ=255)
          agcagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCagag                                                                                                                                       >  1:2100741/1‑150 (MQ=255)
            cagcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGAc                                                                                                                                     >  2:3436674/1‑150 (MQ=255)
             agcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGAcc                                                                                                                                    >  2:4019786/1‑150 (MQ=255)
             agcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGAcc                                                                                                                                    >  2:3948895/1‑150 (MQ=255)
              gcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCt                                                                                                                                   >  2:483108/1‑150 (MQ=255)
              gcaAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCt                                                                                                                                   <  2:2510668/150‑1 (MQ=255)
               caAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAAAGAGTGGTGCAATAAAATCAGAGACCTg                                                                                                                                  >  2:3290932/1‑150 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 >  1:555255/1‑150 (MQ=255)
                aaaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTgg                                                                                                                                 >  1:3231225/1‑150 (MQ=255)
                 aaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGga                                                                                                                                <  2:2300704/150‑1 (MQ=255)
                 aaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGga                                                                                                                                >  1:3965050/1‑150 (MQ=255)
                   atatCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGgaag                                                                                                                              <  2:3999076/150‑1 (MQ=255)
                         tctGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGaa                                                                                                                        >  1:1213153/1‑150 (MQ=255)
                         tctGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGaa                                                                                                                        <  1:3132607/150‑1 (MQ=255)
                          ctGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAAc                                                                                                                       >  1:1162928/1‑150 (MQ=255)
                            gTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAAGGt                                                                                                                     <  2:2392407/150‑1 (MQ=255)
                            gTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGt                                                                                                                     >  2:377776/1‑150 (MQ=255)
                             tACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGtt                                                                                                                    <  1:99836/150‑1 (MQ=255)
                              aCGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGttt                                                                                                                   >  1:3544959/1‑150 (MQ=255)
                              aCGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGttt                                                                                                                   >  1:3544982/1‑150 (MQ=255)
                                gggAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAg                                                                                                                 >  1:1942172/1‑150 (MQ=255)
                                  gAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTa                                                                                                               <  1:1352446/150‑1 (MQ=255)
                                   aGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAc                                                                                                              <  1:4037208/150‑1 (MQ=255)
                                       aGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAcggcg                                                                                                          >  2:3132326/1‑150 (MQ=255)
                                       aGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAcggcg                                                                                                          >  2:2725132/1‑150 (MQ=255)
                                        gggAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAcggcgg                                                                                                         >  2:2807951/1‑150 (MQ=255)
                                        gggAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTAcggcgg                                                                                                         <  2:2005654/150‑1 (MQ=255)
                                          gAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACggcgggc                                                                                                       >  2:1155758/1‑150 (MQ=255)
                                          gAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACggcgggc                                                                                                       >  2:1155762/1‑150 (MQ=255)
                                            aCAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGa                                                                                                     >  2:3416477/1‑150 (MQ=255)
                                            aCAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGa                                                                                                     >  2:2159788/1‑150 (MQ=255)
                                             cAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGaa                                                                                                    >  2:2447822/1‑150 (MQ=255)
                                                   ggCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAat                                                                                              <  1:915365/150‑1 (MQ=255)
                                                     cAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAatca                                                                                            <  1:302794/150‑1 (MQ=255)
                                                       aaCGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAatcatc                                                                                          >  1:1111942/1‑150 (MQ=255)
                                                       aaCGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAatcatc                                                                                          >  1:3751941/1‑150 (MQ=255)
                                                         cGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATcac                                                                                        <  1:2952148/150‑1 (MQ=255)
                                                          gggAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATcacc                                                                                       >  1:2454886/1‑150 (MQ=255)
                                                            gAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATcaccac                                                                                     >  1:1492283/1‑150 (MQ=255)
                                                             aaGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACg                                                                                    >  2:1209953/1‑150 (MQ=255)
                                                              aGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGa                                                                                   >  1:3402066/1‑150 (MQ=255)
                                                              aGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGa                                                                                   >  1:3402352/1‑150 (MQ=255)
                                                               gTTCGGGTCATAGCGGCCCCACATTTCACCGTTAATAATGTTGCCCAGAAGCCCGGCACCAAGAGCAAACGGAATGTGTTGTGCAATAAAATCAGAGACCTGGAAGAATGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGAt                                                                                  <  1:2015529/150‑1 (MQ=255)
                                                                 tCGGGTCACCGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATaa                                                                                <  1:2468668/150‑1 (MQ=255)
                                                                   gggTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACg                                                                              <  2:2600767/150‑1 (MQ=255)
                                                                    ggTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGc                                                                             >  2:2095200/1‑150 (MQ=255)
                                                                    ggTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGc                                                                             >  2:1320480/1‑150 (MQ=255)
                                                                      tCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCa                                                                           >  1:3474730/1‑150 (MQ=255)
                                                                       cAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCaa                                                                          <  1:3167419/150‑1 (MQ=255)
                                                                        aaCGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAAt                                                                         >  1:1452485/1‑150 (MQ=255)
                                                                        aaCGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAAt                                                                         <  2:4157472/150‑1 (MQ=255)
                                                                         aCGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATc                                                                        <  1:463066/150‑1 (MQ=255)
                                                                            cgGCCCCAAAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAgg                                                                     >  1:3569420/1‑150 (MQ=255)
                                                                               ccccACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGgccg                                                                  <  1:735135/150‑1 (MQ=255)
                                                                                        tCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGaaa                                                         >  2:72714/1‑150 (MQ=255)
                                                                                         cACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAg                                                        <  1:208349/150‑1 (MQ=255)
                                                                                          aCCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGa                                                       >  2:3508853/1‑150 (MQ=255)
                                                                                            cGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACa                                                     >  1:2848007/1‑150 (MQ=255)
                                                                                             gctaataaCGTTGCCCAGACGCCCGGCACCAAGACCAACCGGAACGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTGAGTACGGCGGGCGAAGATAATCATCACCACGATACCGCCAATCAGGCCGCCGTGGAAAGACAt                                                    <  2:4128340/148‑1 (MQ=255)
                                                                                              ttAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATg                                                   <  1:2936853/150‑1 (MQ=255)
                                                                                                aataaAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgcc                                                 <  2:3345158/150‑1 (MQ=255)
                                                                                                aataaAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgcc                                                 >  2:1476299/1‑150 (MQ=255)
                                                                                                aataaAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgcc                                                 >  1:1074750/1‑150 (MQ=255)
                                                                                                   aaaGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATgccgcc                                              >  1:3458555/1‑150 (MQ=255)
                                                                                                    aaGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGccgccg                                             <  1:3881536/150‑1 (MQ=255)
                                                                                                       ttGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTcc                                          <  2:3832896/150‑1 (MQ=255)
                                                                                                       gtGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCCCGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTcc                                          <  2:1854218/149‑1 (MQ=255)
                                                                                                          cccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAg                                       >  2:2944431/1‑150 (MQ=255)
                                                                                                           ccAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGa                                      >  2:287515/1‑150 (MQ=255)
                                                                                                            cAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGac                                     >  1:946206/1‑150 (MQ=255)
                                                                                                             aGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGaca                                    >  1:2277863/1‑150 (MQ=255)
                                                                                                              gACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGacac                                   <  1:831001/150‑1 (MQ=255)
                                                                                                              gACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGacac                                   <  1:831038/150‑1 (MQ=255)
                                                                                                               aCGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACg                                  <  2:804873/150‑1 (MQ=255)
                                                                                                                cGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACgg                                 >  1:2676707/1‑150 (MQ=255)
                                                                                                                cGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGAGGAAAGACATGCCGCCGTCCCAGACACgg                                 >  2:289622/1‑150 (MQ=255)
                                                                                                                    cGGCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACa                             <  2:1747823/150‑1 (MQ=255)
                                                                                                                     ggCACCAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAg                            >  1:1696380/1‑150 (MQ=255)
                                                                                                                       cACCAAGACCAAACGGACTGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGAt                          >  2:3298763/1‑150 (MQ=255)
                                                                                                                          cAAGACCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACa                       >  1:1192324/1‑150 (MQ=255)
                                                                                                                              aCCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCgg                   >  2:3172626/1‑150 (MQ=255)
                                                                                                                              aCCAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCgg                   <  2:3156531/150‑1 (MQ=255)
                                                                                                                               ccAAACGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGa                  <  1:129516/150‑1 (MQ=255)
                                                                                                                                 aaaCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATc                <  1:771043/150‑1 (MQ=255)
                                                                                                                                 aaaCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATc                <  1:2567974/150‑1 (MQ=255)
                                                                                                                                 aaaCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATc                >  2:778079/1‑150 (MQ=255)
                                                                                                                                  aaCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCg               <  2:689401/150‑1 (MQ=255)
                                                                                                                                   aCGGAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCgg              <  2:1660391/150‑1 (MQ=255)
                                                                                                                                      gAATGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCa           >  2:1955154/1‑150 (MQ=255)
                                                                                                                                        aTGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATa         <  2:167352/150‑1 (MQ=255)
                                                                                                                                         tGAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATaa        <  1:810665/150‑1 (MQ=255)
                                                                                                                                          gAGTGGTGCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATaaa       <  2:3374054/150‑1 (MQ=255)
                                                                                                                                             tggtgCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTg    >  1:3863675/1‑150 (MQ=255)
                                                                                                                                               gtgCAATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCg  >  1:4078995/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTGGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGACCAAACGGAATGAGTGGTGCGATAAAATCAGAGACCTGGAAGAAGGAACGTTTAGTACGGCGGGCGAAGATAATCATCACCACGATAACGCCAATCAGGCCGCCGTGGAAAGACATGCCGCCGTCCCAGACACGGAACAGATACAGCGGATCGGCCATAAACTGCG  >  CP009273/2958857‑2959149

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: ‑