Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,099,959 C→T 5.5% A35T (GCC→ACC)  PP_3607 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,099,9590CT 5.5% 205.5 / 2.4 73A35T (GCC→ACC) PP_3607membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/40);  new base T (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

GTCAGCATGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGCCCA  >  NC_002947/4099811‑4100106
                                                                                                                                                    |                                                                                                                                                   
gTCAGCATGACCATCTGCGCACGCAGGCCAGAGCCACGGCGCTCGTGGATAAACACCCGCAAGGCTGAGGGGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCAGTGCCGACCAGGAGCAACACTACCTGCTTGTAGCCGGCa                                                                                                                                                    >  1:826168/1‑150 (MQ=255)
    gCATGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttga                                                                                                                                                >  2:617115/1‑150 (MQ=255)
     cATGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgag                                                                                                                                               <  1:1429594/150‑1 (MQ=255)
      attaCCATCTGCGCACGCAGGCCCGCGACACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgagg                                                                                                                                              <  2:399654/147‑1 (MQ=255)
      aTGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgagg                                                                                                                                              <  1:1781087/150‑1 (MQ=255)
      aTGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgagg                                                                                                                                              <  2:623445/150‑1 (MQ=255)
       tGACCATCTGGGCACGCAGGCCCGCGCTACGGCGCTGGTTGATAAACACCCGCCACGCTGAGGTGAAAACGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCTACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggt                                                                                                                                             >  2:645218/1‑150 (MQ=255)
       tGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggg                                                                                                                                             >  2:316971/1‑149 (MQ=255)
          ccATCTGCGCACGCAGGCCCGCGGCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttga                                                                                                                                          <  2:142268/150‑1 (MQ=255)
          ccATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttga                                                                                                                                          <  1:1133838/150‑1 (MQ=255)
          ccATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttga                                                                                                                                          <  2:1229959/150‑1 (MQ=255)
          ccATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttga                                                                                                                                          <  1:1693058/150‑1 (MQ=255)
          ccATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCGGCTTGAAGCCGGCAttgaggttga                                                                                                                                          <  2:1286403/150‑1 (MQ=255)
           cATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGGCCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttgag                                                                                                                                         <  2:164429/150‑1 (MQ=255)
            aTCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCGGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttgaggttgagg                                                                                                                                        >  1:844478/1‑150 (MQ=255)
             tCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGa                                                                                                                                       >  2:749039/1‑150 (MQ=255)
               tGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAAc                                                                                                                                     <  1:1667286/150‑1 (MQ=255)
                               cgcCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTCAGGTTGAGGACCCAGGCACTGCTCAGCa                                                                                                                     <  2:1059489/150‑1 (MQ=255)
                                gcCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCAt                                                                                                                    >  2:337642/1‑150 (MQ=255)
                                    cGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCATCAtcaag                                                                                                                <  2:1756631/150‑1 (MQ=255)
                                    cGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCAtcaag                                                                                                                <  1:1804875/150‑1 (MQ=255)
                                    cGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCAtcaag                                                                                                                <  1:389150/150‑1 (MQ=255)
                                      gcgcTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCAtcaaggt                                                                                                              >  2:1874152/1‑150 (MQ=255)
                                        gcTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCAGCCCAACGCTAAGACGACGTGCTTGATGCCGGCATACAGACAAAGCTACAAGGCGCTGGGTATCAtcaaggtca                                                                                                            >  1:513016/1‑150 (MQ=255)
                                            gTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGc                                                                                                        >  2:1025854/1‑150 (MQ=255)
                                               gATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAt                                                                                                     >  2:1733603/1‑150 (MQ=255)
                                                aTAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGa                                                                                                      <  1:1512680/148‑1 (MQ=255)
                                                aTAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGa                                                                                                      >  2:1512680/1‑148 (MQ=255)
                                                  aaaCACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGTGACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCaa                                                                                                  >  1:1349017/1‑150 (MQ=255)
                                                     cccccGCCAGGCCTAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCATCAACACTCCGTGCTGGAAGCCGGCATTGAGGTTGAGGAACCAGCCACTTGCCAGCATCAAGGTCAAGGCGATGAAtcc                                                                                               <  2:758113/148‑3 (MQ=255)
                                                      acCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGgc                                                                                              <  1:316971/150‑1 (MQ=255)
                                                      acCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGgc                                                                                              <  1:1765388/150‑1 (MQ=255)
                                                        ccGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGc                                                                                                                                >  2:1248785/1‑114 (MQ=255)
                                                        ccGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCAGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGc                                                                                                                                <  1:1248785/114‑1 (MQ=255)
                                                             aGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTCGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGa                                                                                       <  1:150373/150‑1 (MQ=255)
                                                                     gTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGt                                                                               >  2:566353/1‑150 (MQ=255)
                                                                            cGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAg                                                                        >  1:900533/1‑150 (MQ=255)
                                                                                ggcggcGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGt                                                                    <  1:1874152/150‑1 (MQ=255)
                                                                                 gcggcgTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTg                                                                   <  1:1025854/150‑1 (MQ=255)
                                                                                  cggcgTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTgg                                                                  <  1:566353/150‑1 (MQ=255)
                                                                                              tcatCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGATGCAGGCACTGGCCAGCATCAAGGTCAAGGCGANCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGaaa                                                      <  1:422296/147‑1 (MQ=255)
                                                                                              gCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGaaa                                                      <  1:532113/150‑1 (MQ=255)
                                                                                                gTCGGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAgg                                                    <  2:844478/150‑1 (MQ=255)
                                                                                                gTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAGCAAGGGCGCTCCGAGACGACGAGGTTCAGGCGTGGCAGTTGCGGAAAgg                                                    >  2:1365541/1‑150 (MQ=255)
                                                                                                  cAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTGCCTGCTTGAAGCCGGCATTGAGGTGAAGGAACCAGGCACTGGCCAGCATCCAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGcc                                                  <  1:720147/150‑1 (MQ=255)
                                                                                                     gCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttga                                                                                                                                                >  2:1361728/1‑53 (MQ=255)
                                                                                                     gCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCAttga                                                                                                                                                <  1:1361728/53‑1 (MQ=255)
                                                                                                             cGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCa                                             >  2:20766/1‑144 (MQ=255)
                                                                                                             cGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCa                                             <  1:20766/144‑1 (MQ=255)
                                                                                                               gTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTc                                     <  1:617115/150‑1 (MQ=255)
                                                                                                                tGCCGACCAGCACCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCACGGCGATCAATGGCGCGCCGAAGCGCCGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCAGCAAAATACTCg                                    <  2:1057531/150‑1 (MQ=255)
                                                                                                                         gcagcaACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGt                           >  1:901187/1‑150 (MQ=255)
                                                                                                                          cagcaACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGtt                          <  2:1415325/150‑1 (MQ=255)
                                                                                                                            gcaACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATa                             <  2:506509/145‑1 (MQ=255)
                                                                                                                            gcaACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATa                             >  1:506509/1‑145 (MQ=255)
                                                                                                                                cacTACCTGCTTGAAGCCGGTATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAt                                                                                                     >  1:1645077/1‑69 (MQ=255)
                                                                                                                                cacTACCTGCTTGAAGCCGGTATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAt                                                                                                     >  1:1644952/1‑69 (MQ=255)
                                                                                                                                cacTACCTGCTTGAAGCCGGTATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAt                                                                                                     <  2:1645077/69‑1 (MQ=255)
                                                                                                                                cacTACCTGCTTGAAGCCGGTATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGAt                                                                                                     <  2:1644952/69‑1 (MQ=255)
                                                                                                                                cacTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAacgg                    >  2:1702795/1‑150 (MQ=255)
                                                                                                                                   tACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCa                                                                         <  1:1890358/94‑1 (MQ=255)
                                                                                                                                   tACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCa                                                                         >  2:1890358/1‑94 (MQ=255)
                                                                                                                                     ccTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACggg               <  1:1733603/150‑1 (MQ=255)
                                                                                                                                       tGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGa             >  1:1315465/1‑150 (MQ=255)
                                                                                                                                         cTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGt           <  1:352173/150‑1 (MQ=255)
                                                                                                                                            gAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGtata        >  2:1784713/1‑150 (MQ=255)
                                                                                                                                             aaGCCGGCATTGAGGTTGAGGAACCCGGCACTGGCCAGCATAAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTCCATGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGtataa       <  1:645218/150‑2 (MQ=255)
                                                                                                                                                ccGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGcc    >  1:787976/1‑150 (MQ=255)
                                                                                                                                                 cGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGccc   >  2:1820610/1‑150 (MQ=255)
                                                                                                                                                 cGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGccc   >  2:1820218/1‑150 (MQ=255)
                                                                                                                                                 cGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGccc   >  2:1596377/1‑150 (MQ=255)
                                                                                                                                                  ggCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGCCca  <  1:749039/150‑1 (MQ=255)
                                                                                                                                                  ggCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGCCca  <  2:787976/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCAGCATGACCATCTGCGCACGCAGGCCCGCGCCACGGCGCTCGTTGATAAACACCCGCCAGGCTGAGGTGAAACCGAAGGCGGCGTGGTACAGCGTCAGGCCAAGGGCGGTGCCGACCAGCAGCAACACTACCTGCTTGAAGCCGGCATTGAGGTTGAGGAACCAGGCACTGGCCAGCATCAAGGTCAAGGCGATCAATGGCGCGCCGAGGCGACGTGGTTCAGGCGTGGCAGTTGCGGAAAGGCCCATCAAAATACTCGGACAATAGTTAAACGGACGGGGAGTATACGCCCA  >  NC_002947/4099811‑4100106

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