Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,006 T→C 7.5% S10P (TCG→CCG) ‡ 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,0060TC 7.5% 190.6 / 5.6 67S10P (TCG→CCG) ‡waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base C (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 = 4.88e-01

TCACAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGC  >  NC_000913/3795861‑3796151
                                                                                                                                                 |                                                                                                                                                 
tCACAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGa                                                                                                                                                >  1:1428529/1‑149 (MQ=255)
  acaAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATg                                                                                                                                              <  1:850909/149‑1 (MQ=255)
    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:1045435/1‑149 (MQ=255)
                                                                            gaagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggt                                                                    >  2:561493/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)
                                                                                                                                                 |                                                                                                                                                 
TCACAAAGTGCGTAAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGC  >  NC_000913/3795861‑3796151

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: ‑