Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,685,381 G→A 7.8% intergenic (+20/‑39) PP_5593 → / → PP_5594 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,685,3810GA 7.8% 182.6 / 4.9 64intergenic (+20/‑39)PP_5593/PP_5594hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base A (2/3);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TAACAAGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685232‑3685517
                                                                                                                                                     |                                                                                                                                        
tAACAAGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCCGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCg                                                                                                                                          >  2:655032/1‑150 (MQ=255)
    aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      >  2:1514491/1‑150 (MQ=255)
    aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  1:470096/150‑1 (MQ=255)
    aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCAAGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  2:1532498/150‑1 (MQ=255)
        aCCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGAcc                                                                                                                                  <  1:835517/150‑1 (MQ=255)
         cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  2:985200/145‑1 (MQ=255)
         cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      >  1:985200/1‑145 (MQ=255)
         cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 >  2:1794331/1‑150 (MQ=255)
         cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 <  1:1836110/150‑1 (MQ=255)
          ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  1:659947/150‑1 (MQ=255)
          ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                >  1:1457553/1‑150 (MQ=255)
          ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  2:1688412/150‑1 (MQ=255)
          ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGAAGTTACGGAGCGAccat                                                                                                                                >  2:1564557/1‑148 (MQ=255)
           cGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1737561/150‑1 (MQ=255)
                aaTTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTg                                                                                                                          <  2:44124/150‑1 (MQ=255)
                        cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATc                                                                                                                  >  2:725702/1‑150 (MQ=255)
                        cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       >  1:344882/1‑143 (MQ=255)
                        cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       <  2:344882/145‑3 (MQ=255)
                            gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 <  2:444025/147‑1 (MQ=255)
                            gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 >  1:444025/1‑147 (MQ=255)
                             aCCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACTGAGGCGCAGTGCCGGGGCGACATGGCACTGCAATTATCCAAgg                                                                                                             <  2:275201/150‑1 (MQ=255)
                               cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           >  1:1520962/1‑150 (MQ=255)
                               cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           <  2:1925457/150‑1 (MQ=255)
                                 tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  2:1633268/1‑128 (MQ=255)
                                 tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1633268/128‑1 (MQ=255)
                                            ccGACGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCAACAAGGACGGCCGTTCGTAACTGACCGAGTCTTAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATgcgc                                                                                              <  1:725702/150‑1 (MQ=255)
                                                     gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       >  2:886213/1‑116 (MQ=255)
                                                     gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       <  1:886213/116‑1 (MQ=255)
                                                      tCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCAc                                                                                    <  1:1564557/150‑1 (MQ=255)
                                                       cGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACt                                                                                   <  2:104050/150‑1 (MQ=255)
                                                            aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 <  2:247053/147‑1 (MQ=255)
                                                            aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 >  1:247053/1‑147 (MQ=255)
                                                             ggcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcat                                                                             <  1:431262/150‑1 (MQ=255)
                                                              gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:357792/150‑1 (MQ=255)
                                                              gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1457553/150‑1 (MQ=255)
                                                                  ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  2:254098/1‑146 (MQ=255)
                                                                  ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:254098/146‑1 (MQ=255)
                                                                  ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        <  2:709126/102‑2 (MQ=38)
                                                                  ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        >  1:709126/1‑101 (MQ=38)
                                                                               ggcAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTgcg                                                           >  2:1698902/1‑148 (MQ=255)
                                                                                   aCATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTg                                                       <  1:1794010/150‑1 (MQ=255)
                                                                                      tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            <  2:1095197/110‑1 (MQ=255)
                                                                                      tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            >  1:1095197/1‑110 (MQ=255)
                                                                                           cACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcgg                                               >  2:1584896/1‑150 (MQ=255)
                                                                                                 cgcgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATGATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGttt                                         >  1:1105209/1‑150 (MQ=255)
                                                                                                   cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                <  2:1465636/141‑1 (MQ=255)
                                                                                                   cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                >  1:1465636/1‑141 (MQ=255)
                                                                                                       tCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACcg                                   >  2:862278/1‑150 (MQ=255)
                                                                                                         gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            >  2:224630/1‑123 (MQ=255)
                                                                                                         gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            <  1:224630/123‑1 (MQ=255)
                                                                                                          ggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGa                                >  2:404611/1‑150 (MQ=255)
                                                                                                             cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1345835/103‑1 (MQ=255)
                                                                                                             cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  1:1345835/1‑103 (MQ=255)
                                                                                                              acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  2:531596/51‑1 (MQ=255)
                                                                                                              acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  1:531596/1‑51 (MQ=255)
                                                                                                              acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCCTCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1520962/150‑1 (MQ=255)
                                                                                                              acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1426930/150‑1 (MQ=255)
                                                                                                                aaGGGCGGCCGTTCGTAACTGACCGAGGAGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGa                          >  2:1677846/1‑150 (MQ=255)
                                                                                                                 aGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGac                         <  1:1794331/150‑1 (MQ=255)
                                                                                                                           ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       >  2:1273805/1‑126 (MQ=255)
                                                                                                                           ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       <  1:1273805/126‑1 (MQ=255)
                                                                                                                            tCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTaa              >  1:1931327/1‑150 (MQ=255)
                                                                                                                                  cTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTc        <  1:712385/150‑1 (MQ=255)
                                                                                                                                      ccGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAc    <  1:862278/150‑1 (MQ=255)
                                                                                                                                        gAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAccc  <  2:1931327/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
TAACAAGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685232‑3685517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑