Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,520,565 T→A 10.4% L21Q (CTG→CAG)  yncG → glutathione S‑transferase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,520,5650TA10.4% 780.7 / 66.5 260L21Q (CTG→CAG) yncGglutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base T (119/114);  new base A (10/17);  total (129/131)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 260 total aligned reads displayed.

CAACTGACGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAGCGGCG  >  CP009273/1520419‑1520710
                                                                                                                                                  |                                                                                                                                                 
cAACTGACGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTg                                                                                                                                                  >  2:1870549/1‑148 (MQ=255)
 aaCTGACGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAgg                                                                                                                                                 <  1:610594/148‑1 (MQ=255)
  aCTGACGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGc                                                                                                                                                <  1:2309367/148‑1 (MQ=255)
      aCGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGat                                                                                                                                            >  1:2234717/1‑148 (MQ=255)
        gAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGatat                                                                                                                                          <  1:3310174/148‑1 (MQ=255)
         aaaaaaTGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATAtt                                                                                                                                         >  2:2154515/1‑148 (MQ=255)
             aaTGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCtt                                                                                                                                     <  1:2826975/148‑1 (MQ=255)
             aaTGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCtt                                                                                                                                     <  1:1146643/148‑1 (MQ=255)
             aaTGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCtt                                                                                                                                     >  1:1851625/1‑148 (MQ=255)
               tGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTAt                                                                                                                                   <  2:2103197/148‑1 (MQ=255)
                gACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATc                                                                                                                                  >  2:63932/1‑148 (MQ=255)
                 aCAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCa                                                                                                                                 <  1:1630202/148‑1 (MQ=255)
                       ccTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGt                                                                                                                           <  2:1203279/148‑1 (MQ=255)
                         ttttGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGGGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTg                                                                                                                         <  2:2359272/148‑1 (MQ=255)
                          tttGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGa                                                                                                                        >  1:3033018/1‑148 (MQ=255)
                             gCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGt                                                                                                                     <  1:1970955/148‑1 (MQ=255)
                              cTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTc                                                                                                                    <  2:2198473/148‑1 (MQ=255)
                                 gTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGt                                                                                                                 <  1:1758003/148‑1 (MQ=255)
                                    aaCCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGt                                                                                                              <  1:492974/148‑1 (MQ=255)
                                     aCCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGtt                                                                                                             >  2:2996957/1‑148 (MQ=255)
                                      ccTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGttt                                                                                                            <  1:2483864/148‑1 (MQ=255)
                                      ccTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGttt                                                                                                            <  2:2713898/148‑1 (MQ=255)
                                       cTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGtttt                                                                                                           >  2:1370643/1‑148 (MQ=255)
                                       cTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGtttt                                                                                                           >  2:1314442/1‑148 (MQ=255)
                                          tgtACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGAc                                                                                                        <  2:3013419/148‑1 (MQ=255)
                                          tgtACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGAc                                                                                                        <  1:2420520/148‑1 (MQ=255)
                                            tACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCa                                                                                                      <  2:3222428/148‑1 (MQ=255)
                                             aCTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCAc                                                                                                     >  1:1818285/1‑148 (MQ=255)
                                             aCTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCAc                                                                                                     >  1:142437/1‑148 (MQ=255)
                                              cTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACg                                                                                                    >  2:3099646/1‑148 (MQ=255)
                                                gTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAg                                                                                                  <  2:1966130/148‑1 (MQ=255)
                                                 tCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAgg                                                                                                 >  1:29537/1‑148 (MQ=255)
                                                      cacTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCt                                                                                            >  2:2726187/1‑148 (MQ=255)
                                                        cTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTc                                                                                          <  2:2778294/148‑1 (MQ=255)
                                                         ttAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGATCTTCa                                                                                         >  1:1220538/1‑148 (MQ=255)
                                                          tAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCAc                                                                                        >  2:3204479/1‑148 (MQ=255)
                                                               tttAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAg                                                                                   <  2:1229478/148‑1 (MQ=255)
                                                                ttAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGt                                                                                  <  2:1011859/148‑1 (MQ=255)
                                                                ttAAAAGATTGTGAGGGACATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGt                                                                                  <  1:1002139/148‑1 (MQ=255)
                                                                  aaaaGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTa                                                                                >  1:760141/1‑148 (MQ=255)
                                                                   aaaGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTAt                                                                               >  2:2015056/1‑148 (MQ=255)
                                                                   aaaGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTAt                                                                               <  2:117006/148‑1 (MQ=255)
                                                                     aGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTg                                                                             >  1:687869/1‑148 (MQ=255)
                                                                      gATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGa                                                                            >  2:509641/1‑148 (MQ=255)
                                                                       aTTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGaa                                                                           >  1:1593564/1‑148 (MQ=255)
                                                                         tgtgAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGaaaa                                                                         <  1:2395417/148‑1 (MQ=255)
                                                                         tgtgAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGaaaa                                                                         <  1:479078/148‑1 (MQ=255)
                                                                         tgtgAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGaaaa                                                                         >  2:105820/1‑148 (MQ=255)
                                                                           tgAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAAcc                                                                       <  2:2678856/148‑1 (MQ=255)
                                                                            gAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAAccc                                                                      <  2:584353/148‑1 (MQ=255)
                                                                               gggCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGa                                                                   <  1:2203280/148‑1 (MQ=255)
                                                                               gggCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGa                                                                   <  1:641822/148‑1 (MQ=255)
                                                                                    tatGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCg                                                              >  1:1185854/1‑148 (MQ=255)
                                                                                    tatGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCg                                                              <  2:2830431/148‑1 (MQ=255)
                                                                                      tGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGtt                                                            <  2:476651/148‑1 (MQ=255)
                                                                                       gATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTtg                                                           >  1:2627191/1‑148 (MQ=255)
                                                                                       gATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTtg                                                           >  1:3255872/1‑148 (MQ=255)
                                                                                        aTTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTtgt                                                          <  1:1208841/148‑1 (MQ=255)
                                                                                         ttAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTtgtg                                                         >  2:1670229/1‑148 (MQ=255)
                                                                                         ttAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTtgtg                                                         <  2:1929835/148‑1 (MQ=255)
                                                                                          tAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTgtgt                                                        >  1:2554279/1‑148 (MQ=255)
                                                                                          tAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTgtgt                                                        <  1:294156/148‑1 (MQ=255)
                                                                                           aaGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTc                                                       >  2:2630266/1‑148 (MQ=255)
                                                                                             ggTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAg                                                     <  1:2056206/148‑1 (MQ=255)
                                                                                             ggTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAg                                                     >  2:783571/1‑148 (MQ=255)
                                                                                                gtATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGtt                                                  <  1:736324/148‑1 (MQ=255)
                                                                                                 tATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTc                                                 <  2:2311870/148‑1 (MQ=255)
                                                                                                  aTGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTcc                                                >  2:1412673/1‑148 (MQ=255)
                                                                                                   tGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCg                                               >  1:313558/1‑148 (MQ=255)
                                                                                                    ggCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGa                                              >  2:857641/1‑148 (MQ=255)
                                                                                                        tACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACt                                          >  1:1543585/1‑148 (MQ=255)
                                                                                                        tACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACt                                          >  2:3046705/1‑148 (MQ=255)
                                                                                                         aCCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTg                                         <  1:1066167/148‑1 (MQ=255)
                                                                                                          cccGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTgg                                        <  2:2475620/148‑1 (MQ=255)
                                                                                                           ccGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGc                                       >  1:3220286/1‑148 (MQ=255)
                                                                                                             ggATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCAt                                     <  2:435610/148‑1 (MQ=255)
                                                                                                              gATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCAtt                                    >  2:1477003/1‑148 (MQ=255)
                                                                                                               aTGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTg                                   >  1:1674173/1‑148 (MQ=255)
                                                                                                                  gggCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGaa                                >  2:1043031/1‑148 (MQ=255)
                                                                                                                   ggCTCTACAATCAGTGAGCTGATGCTAACACTGGCCGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGaaa                               <  2:335771/148‑1 (MQ=255)
                                                                                                                       ctACAATCAGTGAGCTGATGCTAACACTGTCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGa                           <  2:2271695/148‑1 (MQ=255)
                                                                                                                       ctACAATCAGTGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGa                           <  2:3289174/148‑1 (MQ=255)
                                                                                                                           aaTCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACgag                       <  1:820228/148‑1 (MQ=255)
                                                                                                                             tCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGAt                     >  1:2752716/1‑148 (MQ=255)
                                                                                                                              cAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGAtt                    <  1:2652316/148‑1 (MQ=255)
                                                                                                                              cAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGAtt                    >  1:3091048/1‑148 (MQ=255)
                                                                                                                              cAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGAtt                    >  1:399707/1‑148 (MQ=255)
                                                                                                                               aGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTa                   <  2:955744/148‑1 (MQ=255)
                                                                                                                               aGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTa                   <  2:1617671/148‑1 (MQ=255)
                                                                                                                                 tGAGCTGATGCTAACACAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGTGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATg                 >  1:555061/1‑148 (MQ=255)
                                                                                                                                        aTGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGgaga          >  2:2972463/1‑148 (MQ=255)
                                                                                                                                          gCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACa        >  1:2498962/1‑148 (MQ=255)
                                                                                                                                           cTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAg       <  2:808151/148‑1 (MQ=255)
                                                                                                                                            tAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACCgc      >  2:2609656/1‑148 (MQ=255)
                                                                                                                                              acacTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcgg    >  2:2065420/1‑148 (MQ=255)
                                                                                                                                              acacTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcgg    >  1:2771846/1‑148 (MQ=255)
                                                                                                                                              acacAGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcgg    >  1:403476/1‑148 (MQ=255)
                                                                                                                                               cacTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCAATGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcggc   >  1:2063882/1‑148 (MQ=255)
                                                                                                                                                acTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcggcg  >  2:116237/1‑148 (MQ=255)
                                                                                                                                                acTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAgcggcg  >  2:1790604/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CAACTGACGAAAAAATGACAAATCCTTTTGCTGGTTAACCTGTGTACTGTCCTACACTTAATCTTTAAAAGATTGTGAGGGGCATATGATTAAGGTGTATGGCGTACCCGGATGGGGCTCTACAATCAGTGAGCTGATGCTAACACTGGCTGATATTCCTTATCAGTTCGTTGATGTCAGTGGTTTTGACCACGAGGGAGCTTCACGCGAGTTATTGAAAACCCTGAACCCGTTGTGTCAGGTTCCGACACTGGCATTGGAAAATGACGAGATTATGACGGAGACAGCGGCG  >  CP009273/1520419‑1520710

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: