Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,245,537 G→T 100% N286K (AAC→AAA plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,245,5370GT99.5% 686.6 / ‑4.2 215N286K (AAC→AAAplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (96/117);  total (96/119)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 215 total aligned reads displayed.

GTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245393‑4245681
                                                                                                                                                |                                                                                                                                                
      cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGCACGCGCGCAGCTTTATGGACg                                                                                                                                         >  2:2020914/1‑148 (MQ=255)
      cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACg                                                                                                                                         >  2:3744035/1‑148 (MQ=255)
      cccGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACg                                                                                                                                         >  2:277053/1‑148 (MQ=255)
        cggcggCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGtt                                                                                                                                       >  2:1839851/1‑148 (MQ=255)
         ggcggcGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTg                                                                                                                                      <  2:2108935/148‑1 (MQ=255)
           cggcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGAt                                                                                                                                    >  1:1474300/1‑148 (MQ=255)
             gcgATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGc                                                                                                                                  <  1:2239288/148‑1 (MQ=255)
                atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               >  1:1714948/1‑148 (MQ=255)
                atatGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCt                                                                                                                               >  2:3672404/1‑148 (MQ=255)
                     cggcggCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTaa                                                                                                                          <  2:767974/148‑1 (MQ=255)
                          gcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCg                                                                                                                     >  1:1493096/1‑148 (MQ=255)
                          gcACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCg                                                                                                                     <  1:982228/148‑1 (MQ=255)
                             ccAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGtt                                                                                                                  <  1:317068/148‑1 (MQ=255)
                              cAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTc                                                                                                                 <  2:796275/148‑1 (MQ=255)
                                 ccccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAg                                                                                                              >  2:1063463/1‑148 (MQ=255)
                                 ccccTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGCTCCAg                                                                                                              >  2:1480342/1‑148 (MQ=255)
                                     tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                          <  1:3790556/148‑1 (MQ=255)
                                     tGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGa                                                                                                          >  1:3582505/1‑148 (MQ=255)
                                      ggTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGaa                                                                                                         >  1:3730990/1‑148 (MQ=255)
                                      ggTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGaa                                                                                                         >  1:3123707/1‑148 (MQ=255)
                                       gTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                        <  2:2816421/148‑1 (MQ=255)
                                       gTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAg                                                                                                        <  2:2816409/148‑1 (MQ=255)
                                         gATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGcc                                                                                                      <  2:2958351/148‑1 (MQ=255)
                                           tACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCa                                                                                                    >  1:3818162/1‑148 (MQ=255)
                                           tACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCa                                                                                                    >  1:3087168/1‑148 (MQ=255)
                                            aCAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAg                                                                                                   >  1:2756933/1‑148 (MQ=255)
                                             cAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGa                                                                                                  <  1:2327124/148‑1 (MQ=255)
                                              aGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGaa                                                                                                 <  2:2186681/148‑1 (MQ=255)
                                                cACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATa                                                                                               <  1:857394/148‑1 (MQ=255)
                                                    tAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgt                                                                                           <  1:3352366/148‑1 (MQ=255)
                                                      aGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtca                                                                                         >  1:1561979/1‑148 (MQ=255)
                                                       gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcag                                                                                        <  1:2184021/148‑1 (MQ=255)
                                                       gAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcag                                                                                        >  2:380110/1‑148 (MQ=255)
                                                        aaaGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagt                                                                                       >  1:3021630/1‑148 (MQ=255)
                                                          agcagcAGGTAGTCGATGTGACTGCGGTGTCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATTCGgtcagtca                                                                                     <  2:3553430/148‑1 (MQ=255)
                                                          agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     >  2:2664940/1‑148 (MQ=255)
                                                          agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     >  2:1772237/1‑148 (MQ=255)
                                                          agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     <  2:2258325/148‑1 (MQ=255)
                                                          agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     >  2:3663397/1‑148 (MQ=255)
                                                          agcagcAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGgtcagtca                                                                                     <  2:3308300/148‑1 (MQ=255)
                                                           gcagcaGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                    >  2:1922369/1‑148 (MQ=255)
                                                           gcagcaGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGtcagtcag                                                                                    >  2:870941/1‑148 (MQ=255)
                                                            cagcagGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTTAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAgg                                                                                   <  2:2062603/148‑1 (MQ=255)
                                                                 ggTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAat                                                                              <  1:902692/148‑1 (MQ=255)
                                                                  gTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAatc                                                                             <  1:548489/148‑1 (MQ=255)
                                                                      tCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCAtct                                                                         >  2:2974251/1‑148 (MQ=255)
                                                                          tgtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTa                                                                     >  2:3574963/1‑148 (MQ=255)
                                                                          tgtgACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTa                                                                     <  1:1839851/148‑1 (MQ=255)
                                                                              aCTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAagag                                                                 >  1:1020615/1‑148 (MQ=255)
                                                                                tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               <  2:471441/148‑1 (MQ=255)
                                                                                tGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGaa                                                               <  1:2871701/148‑1 (MQ=255)
                                                                                      ggCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgc                                                         <  2:257053/148‑1 (MQ=255)
                                                                                        cAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTcgccg                                                       >  1:1452373/1‑148 (MQ=255)
                                                                                           ggCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCaa                                                    >  2:3011835/1‑148 (MQ=255)
                                                                                            gCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAt                                                   <  1:1312570/148‑1 (MQ=255)
                                                                                             cacaTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtc                                                  >  1:1552167/1‑148 (MQ=255)
                                                                                               caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                >  1:2812821/1‑148 (MQ=255)
                                                                                               caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAAtctc                                                >  1:1367466/1‑148 (MQ=255)
                                                                                               caTATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCACCGCAAtctc                                                >  2:409025/1‑148 (MQ=255)
                                                                                                atatACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATctct                                               >  2:2642818/1‑148 (MQ=255)
                                                                                                 tataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCtt                                              <  1:1685316/148‑1 (MQ=255)
                                                                                                  ataCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTc                                             >  1:317698/1‑148 (MQ=255)
                                                                                                   taCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTcc                                            >  1:3268535/1‑148 (MQ=255)
                                                                                                   taCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTAGACGCAATCTCttac                                            >  1:507898/1‑146 (MQ=255)
                                                                                                    aCCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCa                                           >  2:751255/1‑148 (MQ=255)
                                                                                                       aGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCa                                        <  2:1714294/148‑1 (MQ=255)
                                                                                                        gCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAg                                       <  2:2614573/148‑1 (MQ=255)
                                                                                                          tCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGACGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTg                                     <  1:2595940/148‑1 (MQ=255)
                                                                                                             tGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCac                                  >  2:1067848/1‑147 (MQ=255)
                                                                                                                ccGTCGTGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcg                               >  2:1940034/1‑148 (MQ=255)
                                                                                                                   tcgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGt                            >  1:3563158/1‑148 (MQ=255)
                                                                                                                    cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           <  1:2420279/148‑1 (MQ=255)
                                                                                                                    cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           <  1:2078040/148‑1 (MQ=255)
                                                                                                                    cgtGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGtt                           <  2:3248653/148‑1 (MQ=255)
                                                                                                                      tGGGCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTct                         >  1:482432/1‑148 (MQ=255)
                                                                                                                       gggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctg                        <  1:2913031/148‑1 (MQ=255)
                                                                                                                        ggCCAGCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTctgc                       <  2:608246/148‑1 (MQ=255)
                                                                                                                             gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  >  2:2690014/1‑148 (MQ=255)
                                                                                                                             gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  <  2:2592713/148‑1 (MQ=255)
                                                                                                                             gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  <  1:265594/148‑1 (MQ=255)
                                                                                                                             gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                  >  1:2046916/1‑148 (MQ=255)
                                                                                                                             gCTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCGTCCATCAGTGCAATCGCGTTCTGCTgcgc                  <  2:674260/148‑1 (MQ=255)
                                                                                                                              cTGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCt                 >  2:3471705/1‑148 (MQ=255)
                                                                                                                               tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAGGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  2:524109/148‑1 (MQ=255)
                                                                                                                               tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  2:1354562/148‑1 (MQ=255)
                                                                                                                               tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  2:3240632/148‑1 (MQ=255)
                                                                                                                               tGGCGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCtt                <  2:2893813/148‑1 (MQ=255)
                                                                                                                                  cGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt             <  2:2111083/148‑1 (MQ=255)
                                                                                                                                  cGAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt             <  2:182213/148‑1 (MQ=255)
                                                                                                                                   gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  2:1529755/148‑1 (MQ=255)
                                                                                                                                   gAACGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTc            <  1:2505788/148‑1 (MQ=255)
                                                                                                                                    aaCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCa           <  2:2444342/148‑1 (MQ=255)
                                                                                                                                    aaCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCa           <  1:3663397/148‑1 (MQ=255)
                                                                                                                                     aCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCAt          >  2:1589846/1‑148 (MQ=255)
                                                                                                                                     aCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCAt          <  1:1559391/148‑1 (MQ=255)
                                                                                                                                     aCGCGCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCAt          <  1:130895/148‑1 (MQ=255)
                                                                                                                                      cgcgCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg         >  1:538770/1‑148 (MQ=255)
                                                                                                                                      cgcgCTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATg         <  1:1063804/148‑1 (MQ=255)
                                                                                                                                       gcgcTCAGCTTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATgg        <  2:1282450/148‑1 (MQ=255)
                                                                                                                                           tCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgaga    <  2:3393144/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GTTGATCCCGGCGGCGATATGCGGCGGCACCAGCCCCTGGTGATACAGCACGTAAGAAAGCAGCAGGTAGTCCATGTGACTGCGGTGGCAAGGCACATATACCAGCTCATGACCGTCGTGGGCCAGCTGGCGAACGCGCTCAGCGTTATGGACGTTGATGCCCTGGTAAAGTCGGTTCCAGGTGAAGCCCAGAATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATT  >  CP009273/4245393‑4245681

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: