Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,211,277 A→G 9.0% V399A (GTG→GCG) ‡ PP_1060 ← glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,211,2770AG 9.0% 176.9 / 5.9 66V399A (GTG→GCG) ‡PP_1060glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base G (3/3);  total (36/30)
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.97e-01

CACCACCAGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAACTCCACGGGC  >  NC_002947/1211132‑1211423
                                                                                                                                                 |                                                                                                                                                  
caccacCAGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTCGGGCCAGTCGGCGCCAATGGCCAGCACGCt                                                                                                                                                <  1:1108314/150‑1 (MQ=255)
 accaccAGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTg                                                                                                                                               <  2:852934/150‑1 (MQ=255)
 accaccAGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTg                                                                                                                                               >  2:1399666/1‑150 (MQ=255)
   caccaGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGTCCAGCGTGCTGGc                                                                                                                                             >  1:293561/1‑150 (MQ=255)
   caccaGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGc                                                                                                                                             >  2:800265/1‑150 (MQ=255)
      caGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCGTGCTGGCCAt                                                                                                                                          <  2:293561/150‑1 (MQ=255)
         cgcgcgcTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTc                                                                                                                                       >  2:63084/1‑150 (MQ=255)
          gcgcgcTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCGTGCTGGccattggc                                                                                                                                      <  1:81662/150‑8 (MQ=255)
            gcgcTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGaa                                                                                                                                    <  1:158313/150‑1 (MQ=255)
             cgcTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAg                                                                                                                                   <  1:657450/150‑1 (MQ=255)
                tGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGgcg                                                                                                                                >  1:1191983/1‑150 (MQ=255)
                    gcAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTgg                                                                                                                            >  1:1010713/1‑149 (MQ=255)
                    gcAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGa                                                                                                                            <  2:75577/150‑1 (MQ=255)
                              ttGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGcc                                                                                                                  <  2:1191983/150‑1 (MQ=255)
                                gCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGAt                                                                                                                                          <  2:663898/124‑1 (MQ=255)
                                gCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGAt                                                                                                                                          >  1:663898/1‑124 (MQ=255)
                                gCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGc                                                                                                                >  1:605109/1‑150 (MQ=255)
                                 cGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCt                                                                                                               <  2:9664/150‑1 (MQ=255)
                                  gTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCGTGCTGGccattg                                                                                                                                        >  2:206531/1‑119 (MQ=255)
                                  gTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCGTGCTGGCCAttg                                                                                                                                        <  1:206531/124‑3 (MQ=255)
                                   tGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTgg                                                                                                             <  2:391956/150‑1 (MQ=255)
                                     gCCACGCCGGGCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGgcg                                                                                                           <  1:429843/150‑1 (MQ=255)
                                              gTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTg                                                                                                  >  1:236179/1‑150 (MQ=255)
                                                        tgttgGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTcgc                                                                                        <  2:605109/150‑1 (MQ=255)
                                                         gttgGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGGCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTcgcg                                                                                       >  2:923849/1‑150 (MQ=255)
                                                          ttgGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGgcg                                                                                                           >  1:784255/1‑129 (MQ=255)
                                                          ttgGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGgcg                                                                                                           <  2:784255/129‑1 (MQ=255)
                                                           tgGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCa                                                                                     <  1:688282/150‑1 (MQ=255)
                                                            ggTGTGGCAGCTTTGCGACTGGATGCAGCCAATGGCAAACATTAAGCGGCGTGCCGAGTTGGCCCAGTCGTCGCCAATGGTACCGACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAg                                                                                    >  1:1169625/1‑150 (MQ=255)
                                                               gtgGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGtt                                                                                 >  2:443784/1‑150 (MQ=255)
                                                                  gCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAg                                                                              >  1:829410/1‑150 (MQ=255)
                                                                  gCAGCTTTGCGACTGGATGCAGCCGATGGAGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGAAAATTCGCCATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCTGGTTCAg                                                                              >  1:1391083/1‑150 (MQ=255)
                                                                       tttGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGACGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGa                                                                         >  1:241097/1‑150 (MQ=255)
                                                                        ttGCGACTGGATTCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGAc                                                                        <  2:459633/150‑1 (MQ=255)
                                                                            gACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAgg                                                                    >  2:8167/1‑150 (MQ=255)
                                                                            gACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAgg                                                                    <  1:91556/150‑1 (MQ=255)
                                                                            gACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAgg                                                                    >  1:1394587/1‑150 (MQ=255)
                                                                                  aTGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGtgttg                                                              <  1:800265/150‑1 (MQ=255)
                                                                                   tGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGtgttgt                                                             <  1:923849/150‑1 (MQ=255)
                                                                                     cAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGc                                                           <  1:63084/150‑1 (MQ=255)
                                                                                      aGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCa                                                          <  1:443784/150‑1 (MQ=255)
                                                                                                     tGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCtt                                           >  2:1093688/1‑150 (MQ=255)
                                                                                                     tGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCtt                                           >  1:345377/1‑150 (MQ=255)
                                                                                                      gAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCtt                                                                                             <  2:1366682/99‑1 (MQ=255)
                                                                                                      gAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCtt                                                                                             >  1:1366682/1‑99 (MQ=255)
                                                                                                      gAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTc                                          >  2:843191/1‑150 (MQ=255)
                                                                                                       aaGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGtgttgt                                                             <  2:619506/130‑1 (MQ=255)
                                                                                                       aaGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGtgttgt                                                             >  1:619506/1‑130 (MQ=255)
                                                                                                       aaGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCa                                         >  2:757501/1‑150 (MQ=255)
                                                                                                                  cGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCAc                              <  2:345377/150‑1 (MQ=255)
                                                                                                                  cGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCAc                              >  1:1203862/1‑150 (MQ=255)
                                                                                                                         gCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGata                       >  1:221875/1‑150 (MQ=255)
                                                                                                                         gCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGata                       >  1:972718/1‑150 (MQ=255)
                                                                                                                         gCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGata                       >  1:1166147/1‑150 (MQ=255)
                                                                                                                          cccAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGgt                                                                  >  1:1305183/1‑106 (MQ=255)
                                                                                                                          cccAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGgt                                                                  <  2:1305183/106‑1 (MQ=255)
                                                                                                                          cccAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGatat                      >  1:940120/1‑150 (MQ=255)
                                                                                                                              gTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGACGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTc                  <  2:241097/150‑1 (MQ=255)
                                                                                                                                cGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCgg                <  1:1399666/150‑1 (MQ=255)
                                                                                                                                 ggCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGt               <  2:1166147/150‑1 (MQ=255)
                                                                                                                                    gcCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGaa            <  2:940120/150‑1 (MQ=255)
                                                                                                                                     ccAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAAc           >  1:787053/1‑150 (MQ=255)
                                                                                                                                        atcgCCAGCGTGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAACTcc        >  2:693239/4‑150 (MQ=255)
                                                                                                                                             cAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAACTCCACggg   >  2:1181693/1‑150 (MQ=255)
                                                                                                                                             cAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGCCACGCCGCTATGGTCGGTGAACTCCACggg   >  2:237413/1‑150 (MQ=255)
                                                                                                                                              aGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAACTCCACGGgc  >  2:999620/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACCACCAGCGCGCGCTGGCGCAGCGGATCTTGCGTGGCCACGCCGGTCGGGCATTTGTTGGTGTGGCAGCTTTGCGACTGGATGCAGCCGATGGCGAACATGAAGCCGCGTGCCGAGTTGGCCCAGTCGGCGCCAATGGCCAGCACGCTGGCGATGTCGAAGGCGCTGACAATCTTGCCGCTGGCGCCCAGCTTGATCTTGTCGCGCAGGTTCAGGCCGACCAGGGTGTTGTGCACGAACAGCAGGCCTTCACGCAGTGGCACGCCGATATGGTCGGTGAACTCCACGGGC  >  NC_002947/1211132‑1211423

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