Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,020,586 G→T 6.0% R305S (CGT→AGT)  yedQ ← putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,020,5860GT 6.0% 607.1 / 12.4 166R305S (CGT→AGT) yedQputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (71/85);  new base T (4/6);  total (75/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 166 total aligned reads displayed.

CGCGCACGCCTTCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTC  >  CP009273/2020442‑2020731
                                                                                                                                                |                                                                                                                                                 
cgcgCACGCCTTCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAgg                                                                                                                                                <  1:235472/148‑1 (MQ=255)
        ccTTCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaaa                                                                                                                                        >  1:977835/1‑148 (MQ=255)
          ttCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAata                                                                                                                                      >  1:2352164/1‑148 (MQ=255)
          ttCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAata                                                                                                                                      <  2:604791/148‑1 (MQ=255)
            cGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATAtt                                                                                                                                    >  2:501342/1‑148 (MQ=255)
                 aaCGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                                                                                                               <  1:2629908/148‑1 (MQ=255)
                      aTGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCATGTTGGa                                                                                                                          <  1:2607327/148‑1 (MQ=255)
                         gACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATgg                                                                                                                       >  2:1381761/1‑148 (MQ=255)
                           cacGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTc                                                                                                                     <  1:2814497/148‑1 (MQ=255)
                             cGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAg                                                                                                                   >  1:212332/1‑148 (MQ=255)
                             cGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAg                                                                                                                   >  1:1924072/1‑148 (MQ=255)
                               cACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGa                                                                                                                 <  1:2871576/148‑1 (MQ=255)
                               cACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGa                                                                                                                 <  2:196577/148‑1 (MQ=255)
                                aCCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                                                                                                                <  2:2857185/148‑1 (MQ=255)
                                aCCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                                                                                                                <  2:722769/148‑1 (MQ=255)
                                    gCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGa                                                                                                            <  2:1631176/148‑1 (MQ=255)
                                    gCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGa                                                                                                            <  2:2622833/148‑1 (MQ=255)
                                     cacaCCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGaa                                                                                                           >  1:595927/1‑148 (MQ=255)
                                      acacCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAag                                                                                                          <  1:2590019/148‑1 (MQ=255)
                                            tCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCa                                                                                                    >  2:3026699/1‑148 (MQ=255)
                                               aaaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaa                                                                                                 <  2:1435159/148‑1 (MQ=255)
                                                aaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaaa                                                                                                >  2:2542093/1‑148 (MQ=255)
                                                aaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaaa                                                                                                >  2:2917619/1‑148 (MQ=255)
                                                 aTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAt                                                                                               >  2:1292180/1‑148 (MQ=255)
                                                  tGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc                                                                                              >  1:363812/1‑148 (MQ=255)
                                                   ggTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtct                                                                                             <  1:3117784/148‑1 (MQ=255)
                                                     tCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCa                                                                                           <  2:2255967/148‑1 (MQ=255)
                                                      ccAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAg                                                                                          >  2:2180430/1‑148 (MQ=255)
                                                          aCGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTa                                                                                      <  2:725472/148‑1 (MQ=255)
                                                          aCGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTa                                                                                      <  2:725492/148‑1 (MQ=255)
                                                             ttCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTg                                                                                   >  2:1310664/1‑148 (MQ=255)
                                                             ttCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTg                                                                                   <  1:1141427/148‑1 (MQ=255)
                                                                 cAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCAt                                                                               >  1:2153204/1‑148 (MQ=255)
                                                                  aGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATa                                                                              >  1:2630661/1‑148 (MQ=255)
                                                                   gCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATag                                                                             >  2:2349224/1‑148 (MQ=255)
                                                                     tAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATagag                                                                           <  2:1292245/148‑1 (MQ=255)
                                                                      aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGc                                                                          <  1:1564086/148‑1 (MQ=255)
                                                                      aaCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGc                                                                          <  1:1086779/148‑1 (MQ=255)
                                                                       aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCt                                                                         >  2:3253302/1‑148 (MQ=255)
                                                                       aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCt                                                                         >  2:1770477/1‑148 (MQ=255)
                                                                         aTAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGa                                                                       <  2:3027982/148‑1 (MQ=255)
                                                                           aGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATa                                                                     <  1:2752088/148‑1 (MQ=255)
                                                                             cGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACt                                                                   <  1:150467/148‑1 (MQ=255)
                                                                                cTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCAc                                                                >  1:3126056/1‑148 (MQ=255)
                                                                                 tGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCAcc                                                               <  2:2563191/148‑1 (MQ=255)
                                                                                   tCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCAt                                                             <  1:1412431/148‑1 (MQ=255)
                                                                                   tCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCAt                                                             <  2:1740023/148‑1 (MQ=255)
                                                                                     cATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                           <  2:780654/148‑1 (MQ=255)
                                                                                     cATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                           <  1:132942/148‑1 (MQ=255)
                                                                                     cATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                           <  1:132884/148‑1 (MQ=255)
                                                                                      aTACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGa                                                          >  2:224551/1‑148 (MQ=255)
                                                                                       tACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAg                                                         >  2:2969564/1‑148 (MQ=255)
                                                                                        aCGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTGACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAgg                                                        <  1:3161734/148‑1 (MQ=255)
                                                                                        aCGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAATATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTCCTGTCATAGAGCTGATACTCACCATCGAgg                                                        >  2:1920506/1‑148 (MQ=255)
                                                                                        aCGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAgg                                                        >  2:1920732/1‑148 (MQ=255)
                                                                                        aCGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAgg                                                        <  1:2292698/148‑1 (MQ=255)
                                                                                          gAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGtt                                                      >  2:814045/1‑148 (MQ=255)
                                                                                           aaTGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGttt                                                     <  2:1214963/148‑1 (MQ=255)
                                                                                           aaTGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGttt                                                     <  1:2299455/148‑1 (MQ=255)
                                                                                            aTGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGtttt                                                    <  1:1873688/148‑1 (MQ=255)
                                                                                            aTGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGtttt                                                    <  1:2727182/148‑1 (MQ=255)
                                                                                              gccgccCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTa                                                  >  2:382168/1‑148 (MQ=255)
                                                                                                cgccCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATc                                                <  2:150779/148‑1 (MQ=255)
                                                                                                  ccccGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGa                                              <  2:691331/148‑1 (MQ=255)
                                                                                                  ccccGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGa                                              >  1:911445/1‑148 (MQ=255)
                                                                                                      gcgTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGc                                          <  2:2996409/148‑1 (MQ=255)
                                                                                                       cgTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCg                                         >  1:147536/1‑148 (MQ=255)
                                                                                                         tgtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGtt                                       >  1:1941398/1‑148 (MQ=255)
                                                                                                          gtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGttt                                      <  2:363812/148‑1 (MQ=255)
                                                                                                          gtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGttt                                      <  1:2607852/148‑1 (MQ=255)
                                                                                                          gtCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGttt                                      >  1:2481038/1‑148 (MQ=255)
                                                                                                           tCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTc                                     <  1:2178650/148‑1 (MQ=255)
                                                                                                            cATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCt                                    >  1:1598682/1‑148 (MQ=255)
                                                                                                              tGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTa                                  <  1:866657/148‑1 (MQ=255)
                                                                                                                 tCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaaa                               <  2:220742/148‑1 (MQ=255)
                                                                                                                 tCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaaa                               >  2:2776540/1‑148 (MQ=255)
                                                                                                                 tCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaaa                               <  1:1041688/148‑1 (MQ=255)
                                                                                                                 tCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaaa                               <  1:1041858/148‑1 (MQ=255)
                                                                                                                 tCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTGTTATCGATGGCGTTTCTTaaaa                               >  2:2776650/1‑148 (MQ=255)
                                                                                                                  ccATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTaaaaa                              <  1:2971572/148‑1 (MQ=255)
                                                                                                                     ttGCCTGCGCCAGCAAGGACAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtg                           <  1:3094445/148‑1 (MQ=255)
                                                                                                                      tGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgc                          <  2:851605/148‑1 (MQ=255)
                                                                                                                        ccTGCGCCAGCAAGGCCAGTTCACTAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgctg                        >  1:3106890/1‑148 (MQ=255)
                                                                                                                        ccTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgctg                        >  2:2477442/1‑148 (MQ=255)
                                                                                                                          tGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCa                      >  1:2466896/1‑148 (MQ=255)
                                                                                                                              ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAACTCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGt                  <  2:94918/148‑1 (MQ=255)
                                                                                                                              ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGt                  >  2:27359/1‑148 (MQ=255)
                                                                                                                              ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGt                  <  2:1794830/148‑1 (MQ=255)
                                                                                                                              ccAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGt                  >  2:1752840/1‑148 (MQ=255)
                                                                                                                              acAGCAAGGCCAGTTCAAGAGGATCAAAAATATTCCCTGTTGTATGGTCAGGATTGGAACAGGTCAAAAATCTCACCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATAGATGGCCTTTCTTAAAAATTCCTGCATAGt                  >  2:27097/2‑148 (MQ=255)
                                                                                                                                 gCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACg               >  2:1138286/1‑148 (MQ=255)
                                                                                                                                  cAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGc              >  1:1907858/1‑148 (MQ=255)
                                                                                                                                  cAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGc              >  1:1823097/1‑148 (MQ=255)
                                                                                                                                  cAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGc              >  1:151718/1‑148 (MQ=255)
                                                                                                                                   aaGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCa             <  1:1125191/148‑1 (MQ=255)
                                                                                                                                   aaGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTNACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCa             >  1:1676695/1‑148 (MQ=255)
                                                                                                                                     ggCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCATAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACg           <  2:2669072/148‑1 (MQ=255)
                                                                                                                                            tCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATAc    >  2:2001311/1‑148 (MQ=255)
                                                                                                                                            tCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATAc    >  1:334756/1‑148 (MQ=255)
                                                                                                                                             cACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACt   <  1:501342/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGCGCACGCCTTCGCTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTC  >  CP009273/2020442‑2020731

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: