Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,797,363 A→C 8.9% L198R (CTC→CGC)  gltB ← L‑glutamate synthase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,797,3630AC 8.9% 201.0 / 11.8 79L198R (CTC→CGC) gltBL‑glutamate synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base C (2/5);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TGATGGTGTTGATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTCAGGCCTTC  >  NC_002947/5797214‑5797511
                                                                                                                                                     |                                                                                                                                                    
tGATGGTGTTGATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCga                                                                                                                                                      <  1:735183/150‑1 (MQ=255)
      tgtTGATCTCGCTGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATcgc                                                                                                                                                >  1:413807/1‑150 (MQ=255)
       gtTGATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATcgcg                                                                                                                                               >  1:136283/1‑150 (MQ=255)
        ttGATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATcgcgc                                                                                                                                              >  2:482048/1‑150 (MQ=255)
         tGATCTCGCTGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATcgcgcg                                                                                                                                             <  2:413807/150‑1 (MQ=255)
          gATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGTCGCGAGATCGCGCgg                                                                                                                                            >  2:989633/1‑150 (MQ=255)
             ctcGCCGTTGTGGGTGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgt                                                                                                                                         <  2:310735/150‑3 (MQ=255)
             ctcGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgt                                                                                                                                         <  1:539317/150‑3 (MQ=255)
             ctcGCCGCTGTGGGCGAGGAATCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgt                                                                                                                                         <  2:615962/150‑3 (MQ=255)
              tcGGCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgtt                                                                                                                                        <  2:152671/150‑4 (MQ=255)
              tcGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatc                                                                                                                                        >  2:1179266/1‑150 (MQ=255)
                gCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatcat                                                                                                                                      <  2:380348/150‑1 (MQ=255)
                gCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatcat                                                                                                                                      >  2:550212/1‑150 (MQ=255)
                gCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatcat                                                                                                                                      >  1:328583/1‑150 (MQ=255)
                     tgtgGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGc                                                                                                                                 <  2:1326395/150‑1 (MQ=255)
                      gtgGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTTGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGcc                                                                                                                                <  1:374274/150‑1 (MQ=255)
                      gtgGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGcc                                                                                                                                >  2:1275362/1‑150 (MQ=255)
                              ggAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAga                                                                                                                        <  1:947376/150‑1 (MQ=255)
                              ggAAGCGGAACGGCTGCGCCACCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCAGCATCAGGCCTTTGTAga                                                                                                                        <  2:731183/150‑1 (MQ=255)
                                   cGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGatgatg                                                                                                                   <  2:1423291/150‑1 (MQ=255)
                                     gAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgt                                                                                                                                         >  1:191008/1‑124 (MQ=255)
                                     gAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGCGATCTCGCGGcgt                                                                                                                                         <  2:191008/126‑3 (MQ=255)
                                      aaCGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTc                                                                                                                >  2:948577/1‑150 (MQ=255)
                                      aaCGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTc                                                                                                                <  2:1080750/150‑1 (MQ=255)
                                          gCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTtgt                                                                                                            <  1:1075650/150‑1 (MQ=255)
                                          gCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTtgt                                                                                                            <  1:806509/150‑1 (MQ=255)
                                          gCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTtgt                                                                                                            >  2:1092355/1‑150 (MQ=255)
                                          gCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTtgt                                                                                                            <  2:328583/150‑1 (MQ=255)
                                          gCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTtgt                                                                                                            <  1:1092355/150‑1 (MQ=255)
                                            tGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGc                                                                                                          <  1:1275834/150‑1 (MQ=255)
                                                  aGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAg                                                                                                    >  1:1173673/1‑150 (MQ=255)
                                                  aGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAg                                                                                                    <  1:989633/150‑1 (MQ=255)
                                                  aGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAg                                                                                                    <  2:428246/150‑1 (MQ=255)
                                                    cGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatcat                                                                                                                                      >  2:606605/1‑114 (MQ=255)
                                                    cGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGcatcat                                                                                                                                      <  1:606605/114‑1 (MQ=255)
                                                    cGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCt                                                                                                  >  2:304450/1‑150 (MQ=255)
                                                         aTTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGa                                                                                             >  1:1305283/1‑150 (MQ=255)
                                                               gtagtgtgtTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGTTCGCGCGGCATCATCATGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGt                                                                                       <  1:555303/145‑1 (MQ=255)
                                                               gCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGt                                                                                       >  1:1219014/1‑150 (MQ=255)
                                                               gCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGt                                                                                       <  1:550212/150‑1 (MQ=255)
                                                                 aGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTgg                                                                                     <  1:333723/150‑1 (MQ=255)
                                                                   ggTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTc                                                                                   <  1:482048/150‑1 (MQ=255)
                                                                      gtTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGc                                                                                >  2:429119/1‑150 (MQ=255)
                                                                       ttGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGcg                                                                               <  1:1108864/150‑1 (MQ=255)
                                                                             gagaAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTgg                                                                         <  2:717627/150‑1 (MQ=255)
                                                                                  gcgcTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATc                                                                                                                                                  >  1:126596/1‑72 (MQ=255)
                                                                                  gcgcTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATc                                                                                                                                                  <  2:126596/72‑1 (MQ=255)
                                                                                  gcgcTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGtt                                                                    <  2:108465/150‑1 (MQ=255)
                                                                                   cgcTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGCAGGTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTg                                                                   <  1:275096/150‑1 (MQ=255)
                                                                                   cgcTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTg                                                                   <  1:304450/150‑1 (MQ=255)
                                                                                      tggtggAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAg                                                                                              >  1:232024/1‑120 (MQ=255)
                                                                                      tggtggAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAg                                                                                              <  2:232024/120‑1 (MQ=255)
                                                                                           gAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACcga                                                           >  1:1080639/1‑150 (MQ=255)
                                                                                               acGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGt                                                                                    >  2:827316/1‑121 (MQ=255)
                                                                                               acGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGt                                                                                    <  1:827316/121‑1 (MQ=255)
                                                                                                 gCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGacgaacg                                                     >  2:1334600/1‑150 (MQ=255)
                                                                                                      tCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAacgacgac                                                <  1:1179266/150‑1 (MQ=255)
                                                                                                      tCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAacgacgac                                                <  1:429119/150‑1 (MQ=255)
                                                                                                        gcgGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACgg                                              <  2:136283/150‑1 (MQ=255)
                                                                                                              tGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTg                                        <  2:61444/150‑1 (MQ=255)
                                                                                                                 aGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCGGGGCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGACGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAAc                                     >  1:591225/1‑150 (MQ=255)
                                                                                                                   gcgcTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAg                                   >  1:1367022/1‑150 (MQ=255)
                                                                                                                   gcgcTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAg                                   >  1:914787/1‑150 (MQ=255)
                                                                                                                     gcTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTg                                                                                                 <  2:25092/86‑1 (MQ=255)
                                                                                                                     gcTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTg                                                                                                 >  1:25092/1‑86 (MQ=255)
                                                                                                                     gcTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCt                                 <  1:948577/150‑1 (MQ=255)
                                                                                                                         gTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGGCACCGACGAACGACGACGGGCACTGAACAGCTTGAt                             >  2:643664/1‑150 (MQ=255)
                                                                                                                          tCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTg                                                                                                                            <  1:499355/54‑1 (MQ=255)
                                                                                                                          tCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTg                                                                                                                            >  2:499355/1‑54 (MQ=255)
                                                                                                                              ccAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTg                                                                          >  1:1050898/1‑100 (MQ=255)
                                                                                                                              ccAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTg                                                                          <  2:1050898/100‑1 (MQ=255)
                                                                                                                              ccAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACgg                                              <  2:252928/128‑1 (MQ=255)
                                                                                                                              ccAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACgg                                              >  1:252928/1‑128 (MQ=255)
                                                                                                                                   tCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTcc                   >  1:33117/1‑150 (MQ=255)
                                                                                                                                      ggATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTgg                <  2:914787/150‑1 (MQ=255)
                                                                                                                                            aaCGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTc          >  1:1257314/1‑150 (MQ=255)
                                                                                                                                              tgccgcGCGATCTCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTCAg        >  2:642439/2‑150 (MQ=255)
                                                                                                                                               gccgcGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTCAgg       <  2:1080639/150‑1 (MQ=255)
                                                                                                                                                   cgAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGAt                                                                                            >  2:606783/1‑61 (MQ=255)
                                                                                                                                                   cgAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGAt                                                                                            <  1:606783/61‑1 (MQ=255)
                                                                                                                                                    gagaTCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTCAGGCCTTc  <  2:1219014/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATGGTGTTGATCTCGCCGTTGTGGGCGAGGAAGCGGAACGGCTGCGCCAGCGGCCATTTCGGCAGGGTGTTGGTGGAGAAGCGCTGGTGGAACACGCAGATCGCGGTTTGCAGGCGCTCGTCACCCAGGTCTGGATAAAACGCCGCGAGATCGCGCGGCATCATCAGGCCTTTGTAGATGATGGTCTTGTGCGAGAAGCTGCAGATGTAGTGGTCGGCGTCGTGGGCGTTGGCCACCGACGAACGACGACGGGCACTGAACAGCTTGATGGCGAATTCCTGGTCGCTCAGGCCTTC  >  NC_002947/5797214‑5797511

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: