Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,013 G→A 7.4% G12D (GGC→GAC)  waaC → ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,0130GA 7.4% 202.0 / 4.9 67G12D (GGC→GAC) waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base G (35/27);  new base A (3/2);  total (38/29)
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.79e-01

AAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAC  >  NC_000913/3795865‑3796156
                                                                                                                                                    |                                                                                                                                               
aaaGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATggg                                                                                                                                                 <  2:1415553/149‑1 (MQ=255)
      cGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGt                                                                                                                                           >  2:835013/1‑149 (MQ=255)
          atggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCGCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtctc                                                                                                                                       <  2:729365/147‑1 (MQ=255)
          aaGGTGACGCTGCGTTGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtctc                                                                                                                                       <  2:1448279/149‑1 (MQ=255)
          aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtctc                                                                                                                                       >  2:107908/1‑149 (MQ=255)
          aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAAATCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtctc                                                                                                                                       <  2:1242616/149‑1 (MQ=255)
           cggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGCCCATCGACGATGttttaa                                                                                                                                      <  1:1275035/148‑5 (MQ=255)
             gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCAt                                                                                                                                    >  1:1060521/1‑149 (MQ=255)
             gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCAt                                                                                                                                    >  2:415021/1‑149 (MQ=255)
                aCGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACg                                                                                                                                 <  1:1529818/149‑1 (MQ=255)
                  gCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGTAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGtt                                                                                                                               <  2:1570095/149‑1 (MQ=255)
                  gCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGtt                                                                                                                               <  1:340173/149‑1 (MQ=255)
                             ttATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGcactca                                                                                                                    >  1:1341229/1‑149 (MQ=255)
                              tATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGcactcac                                                                                                                   <  1:514122/149‑1 (MQ=255)
                              tATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGcactcac                                                                                                                   <  1:1367288/149‑1 (MQ=255)
                                  aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGAt                                                                                                               >  1:712630/1‑149 (MQ=255)
                                    cAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGc                                                                                                             >  2:175767/1‑149 (MQ=255)
                                     agagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGcc                                                                                                            <  1:1035422/149‑1 (MQ=255)
                                         cTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagc                                                                                                        >  2:1425718/1‑149 (MQ=255)
                                           tAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTTCCCGCACTCACTGATGCCcagcag                                                                                                      <  1:107908/149‑1 (MQ=255)
                                           tAATCGACATTACTCCCCAGCGCGAACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCATGACGGATGCGGGTTTTGATCGTTAAAACATCCTCGCTGGGCGATGTTCTCCATACGTTGCACGCACTCACTGATGCCcagcag                                                                                                      >  1:1140721/1‑149 (MQ=255)
                                            aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTAACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACAACGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAgg                                                                                                     <  1:1000287/149‑1 (MQ=255)
                                             aTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                    >  1:905933/1‑149 (MQ=255)
                                              tCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCa                                                                                                   >  2:1532097/1‑149 (MQ=255)
                                               cGACATTACTCCCCAGCGCGTACTGGAAGAACTCAGCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCaa                                                                                                  >  1:705387/1‑149 (MQ=255)
                                                gACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAAt                                                                                                 >  2:746306/1‑149 (MQ=255)
                                                      aCTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAgg                                                                                           <  2:1060521/149‑1 (MQ=255)
                                                        tCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGa                                                                                         <  1:1023794/149‑1 (MQ=255)
                                                            cAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGacact                                                                                     >  2:257723/1‑145 (MQ=255)
                                                              gcgcGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGt                                                                                   <  2:864414/149‑1 (MQ=255)
                                                                    aCTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACt                                                                             >  1:489341/1‑149 (MQ=255)
                                                                        gaagcaCTCAACGCGCTATTGTGACAAGAGGAAGCCTGACGGATGCGGGTTTGGATCGTTACAACATCGTCGATGGGCGATGTACTCCATACGGTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTATGACTgggg                                                                         >  1:978034/1‑148 (MQ=255)
                                                                        gaagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                         >  2:561493/1‑149 (MQ=255)
                                                                        gaagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                         >  2:1045435/1‑149 (MQ=255)
                                                                        gaagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                         >  1:590369/1‑149 (MQ=255)
                                                                              cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGga                                                                   >  2:373591/1‑149 (MQ=255)
                                                                               tCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                  >  2:1526140/1‑149 (MQ=255)
                                                                               tCAACGCACTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                  >  1:8622/1‑149 (MQ=255)
                                                                                  aCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGGGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaagaa                                                               <  2:1333709/149‑1 (MQ=255)
                                                                                      gcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                           <  1:1190338/149‑1 (MQ=255)
                                                                                        tATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTACATACGTTGCCCGCACTCACTCATGCCCAGCCGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTc                                                         >  2:257688/1‑149 (MQ=255)
                                                                                               aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                  <  1:737532/149‑1 (MQ=255)
                                                                                               aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGAAGGCAATCCCAGGGATTAAGATTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                  >  1:113272/1‑149 (MQ=255)
                                                                                                cAAGAGGACGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCTGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                 <  2:155129/149‑1 (MQ=255)
                                                                                                cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGTTTAAGTTTGACTGGGTGGTGGAAGAAGCGATCGCCCAGAt                                                 >  2:668452/1‑149 (MQ=255)
                                                                                                    agGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcct                                             <  1:394754/149‑1 (MQ=255)
                                                                                                      gAAGCCTGACGGATGTGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttc                                           >  1:496218/1‑149 (MQ=255)
                                                                                                      gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAACTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttc                                           <  1:119458/149‑1 (MQ=255)
                                                                                                                       ggTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTg                          <  2:1333638/149‑1 (MQ=255)
                                                                                                                        gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGa                         <  2:8622/149‑1 (MQ=255)
                                                                                                                            tGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGa                     <  1:835013/149‑1 (MQ=255)
                                                                                                                              aTCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                   >  1:898450/1‑149 (MQ=255)
                                                                                                                              aTCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGt                   >  1:1280712/1‑149 (MQ=255)
                                                                                                                               tCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGtt                  >  1:1500597/1‑149 (MQ=255)
                                                                                                                                 gTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAt                <  2:392780/149‑1 (MQ=255)
                                                                                                                                  ttAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt               >  2:1142020/1‑149 (MQ=255)
                                                                                                                                  ttAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt               >  2:610948/1‑149 (MQ=255)
                                                                                                                                  ttAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt               >  2:638302/1‑149 (MQ=255)
                                                                                                                                  ttAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt               >  2:720801/1‑149 (MQ=255)
                                                                                                                                     gaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCTCTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtggt                                                                      <  1:260926/90‑1 (MQ=255)
                                                                                                                                     gaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtggt                                                                      >  2:260926/2‑91 (MQ=255)
                                                                                                                                     gaaCATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACGCGGTGGTCGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCt            >  1:140265/2‑149 (MQ=255)
                                                                                                                                     aaaCATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGACGTAttgcc            >  2:287170/1‑146 (MQ=255)
                                                                                                                                     aaaCATCGGCGATGGGCGAGGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCt            <  1:1063633/149‑1 (MQ=255)
                                                                                                                                        cATCGCCCATCGACGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtgtg         >  2:872272/1‑149 (MQ=255)
                                                                                                                                          tcgtcgATGGGCGATGTTCTCCATACGGTGCCCGCACTCACTGATGCCCAGCAGGCATTCCCTGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCAGGCTGCCGTTGAGCGAGTTTTTCCTGTGGc       <  1:257688/149‑1 (MQ=255)
                                                                                                                                               gATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc  <  2:705387/149‑1 (MQ=255)
                                                                                                                                               gATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc  >  1:1172816/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAC  >  NC_000913/3795865‑3796156

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