Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,024,313 G→C 6.4% S279R (AGC→AGG opdP ← glycine‑glutamate dipeptide porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,024,3130GC 6.4% 246.4 / 5.3 77S279R (AGC→AGGopdPglycine‑glutamate dipeptide porin
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base C (2/3);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTCCAC  >  NC_002947/1024164‑1024462
                                                                                                                                                     |                                                                                                                                                     
tCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAgg                                                                                                                                                       >  1:462007/1‑150 (MQ=255)
tCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAgg                                                                                                                                                       <  1:41931/150‑1 (MQ=255)
   ttGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTg                                                                                                                                                    <  2:730111/150‑1 (MQ=255)
   ttGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGCCTg                                                                                                                                                    >  2:1405556/1‑150 (MQ=255)
        gTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCAc                                                                                                                                               >  1:1074320/1‑150 (MQ=255)
         ttGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGCGTGTTGGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACt                                                                                                                                              <  1:1022931/150‑1 (MQ=255)
          tGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGCCTGAGCAccg                                                                                                                                             <  1:1323299/150‑3 (MQ=255)
          tGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGCCTGAGCAccg                                                                                                                                             <  1:186767/150‑3 (MQ=255)
           gACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTgc                                                                                                                                            <  2:1323773/150‑1 (MQ=255)
           gACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTgc                                                                                                                                            >  1:284002/1‑150 (MQ=255)
           gACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGCCTGAGCAccgg                                                                                                                                            <  1:1405556/150‑4 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                         >  2:231550/1‑150 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                         >  2:1041096/1‑150 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                         <  2:102987/150‑1 (MQ=255)
                gctgcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGt                                                                                                                                       <  1:372682/150‑1 (MQ=255)
                      aggaggCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGTGCTGTCGCCGa                                                                                                                                 >  1:9813/1‑150 (MQ=255)
                            cgcTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACAt                                                                                                                           <  1:1424634/150‑1 (MQ=255)
                             gcTGAGGGTGTGGGCTTGATGGGTCCGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATg                                                                                                                          <  1:362035/150‑1 (MQ=255)
                             gcTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATg                                                                                                                          >  2:362035/1‑150 (MQ=255)
                              cTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATgg                                                                                                                         <  2:1158841/150‑1 (MQ=255)
                                 aGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTg                                                                                                                      <  2:427907/150‑1 (MQ=255)
                                  gggTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAgg                                                                                                                                                       <  1:526701/116‑1 (MQ=255)
                                  gggTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAgg                                                                                                                                                       >  2:526701/1‑116 (MQ=255)
                                          gCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAag                                                                                                             <  1:231550/150‑1 (MQ=255)
                                            ttGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGgc                                                                                                                    >  2:6653/1‑141 (MQ=255)
                                            ttGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGgc                                                                                                                    <  1:6653/141‑1 (MQ=255)
                                            ttGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                           <  2:678995/150‑1 (MQ=255)
                                            ttGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                           <  1:593401/150‑1 (MQ=255)
                                               aTGCGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAgtagta                                                                                                        <  2:462007/150‑1 (MQ=255)
                                                      aGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTc                                                                                                 >  2:634994/1‑150 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGa                                                                                             >  2:639649/1‑150 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGa                                                                                             >  2:379261/1‑150 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGa                                                                                             >  2:382877/1‑150 (MQ=255)
                                                           cactcaCAGGCTGTAATTTTCGTTTTCGATCTTTCCCATTTTCTGGCTACCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTTAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGCTTCCAGaa                                                                                            <  2:518836/145‑1 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                            >  1:786193/1‑150 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCAGAAGTAGTACTGGTTCCAGaa                                                                                            >  2:865276/1‑150 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGACGGCATAGACGTGCTCTTAGGAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                            <  1:878765/150‑1 (MQ=255)
                                                           cAGTGACAGGCGGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGTGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                            <  1:181673/150‑1 (MQ=255)
                                                               gACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTcc                                                                                        <  2:786193/150‑1 (MQ=255)
                                                                 cAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTc                                                                                      <  1:183408/150‑1 (MQ=255)
                                                                  aGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTc                                                                                      <  1:595616/149‑1 (MQ=255)
                                                                  aGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTc                                                                                      >  2:595616/1‑149 (MQ=255)
                                                                      tggaaGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGGAGCGACGGTATTGCCTGAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGGGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCt                                                                                 >  2:760645/4‑150 (MQ=255)
                                                                        tAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTtg                                                                               <  1:848412/150‑1 (MQ=255)
                                                                                    tCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTgg                                                                                                                     >  1:734787/1‑100 (MQ=255)
                                                                                    tCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTgg                                                                                                                     <  2:734787/100‑1 (MQ=255)
                                                                                     cGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGc                                                                  <  1:1276712/150‑1 (MQ=255)
                                                                                      gATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCt                                                                 >  1:1389879/1‑150 (MQ=255)
                                                                                       aTCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                <  1:1041096/150‑1 (MQ=255)
                                                                                        tCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTgg                                                               >  1:557626/1‑150 (MQ=255)
                                                                                              ccATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                <  1:951274/143‑1 (MQ=255)
                                                                                              ccATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                <  1:951168/143‑1 (MQ=255)
                                                                                              ccATGTTCTTGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                >  2:951274/1‑143 (MQ=255)
                                                                                              ccATGTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                >  2:951168/1‑143 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagtag                                                                                                          <  1:404892/100‑1 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagtag                                                                                                          >  2:404892/1‑100 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                           >  2:1097059/1‑99 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                           <  1:1097059/99‑1 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAg                                                        <  1:379261/150‑1 (MQ=255)
                                                                                                       ggCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGa                                                <  2:1074320/150‑1 (MQ=255)
                                                                                                        gCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGTGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGaa                                               <  2:9813/150‑1 (MQ=255)
                                                                                                          tGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACg                                             >  2:1189752/1‑150 (MQ=255)
                                                                                                              ggCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTg                                         >  2:650826/1‑150 (MQ=255)
                                                                                                                 gTCGACGGTACTGTAGTAGTTCAGGCCGGTGCTCGGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCTCCCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTAAACGGCTggaa                                      <  2:709277/150‑3 (MQ=255)
                                                                                                                    gACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAg                                                                                                                                                  <  2:238392/39‑1 (MQ=255)
                                                                                                                    gACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAg                                                                                                                                                  >  1:238392/1‑39 (MQ=255)
                                                                                                                       ggTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGAt                                <  2:1389879/150‑1 (MQ=255)
                                                                                                                          cTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACACGGTGAGACCCGGGCACTGTAACGTCGGGGGGTTGATGcc                             >  2:373028/1‑150 (MQ=255)
                                                                                                                                         gCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcg              >  1:1319802/1‑150 (MQ=255)
                                                                                                                                          gcagTGCTCAGCCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTCAACGAGCTGAATTAGTTGCTTGCACGGCAGGTAGTTGATGCCGGCGGTGCTGGCAcgg             >  2:397643/4‑150 (MQ=255)
                                                                                                                                           cGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTTGTAGTTGATGCCGGCGGTGCTGGCAcggt            <  1:634994/150‑1 (MQ=255)
                                                                                                                                           cGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcggt            <  2:557626/150‑1 (MQ=255)
                                                                                                                                              tGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcg              <  1:216851/145‑1 (MQ=255)
                                                                                                                                              tGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcg              >  2:216851/1‑145 (MQ=255)
                                                                                                                                                cTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGctct                                                   >  1:186141/1‑106 (MQ=255)
                                                                                                                                                cTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGctct                                                   <  2:186141/106‑1 (MQ=255)
                                                                                                                                                cTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGtt       <  1:382877/150‑1 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                <  1:901223/92‑1 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                                >  2:901223/1‑92 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGttt      >  1:1210363/1‑150 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGttt      >  2:266693/1‑150 (MQ=255)
                                                                                                                                                     gCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTCCAc  >  1:620652/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTCCAC  >  NC_002947/1024164‑1024462

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