Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,105 G→C 7.5% W204C (TGG→TGC) ‡ sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1050GC 7.5% 367.0 / 18.4 106W204C (TGG→TGC) ‡sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (51/47);  new base C (6/2);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 109 total aligned reads displayed.

GAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGC  >  NC_000913/2652958‑2653253
                                                                                                                                                   |                                                                                                                                                    
gAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg                                                                                                                                                      >  3:960526/1‑148 (MQ=255)
            tGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTGGCCCCGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgc                                                                                                                                         >  2:914808/1‑149 (MQ=255)
              tattatTGGCAAGACATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcgc                                                                                                                                       >  2:211833/1‑149 (MQ=255)
               attattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTgcgc                                                                                                                                       >  4:312313/1‑148 (MQ=255)
               attattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGcgcg                                                                                                                                      >  2:905202/1‑149 (MQ=255)
                 tattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGa                                                                                                                                     >  3:1042364/1‑148 (MQ=255)
                   ttGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAg                                                                                                                                   <  4:1042364/148‑1 (MQ=255)
                         aGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcat                                                                                                                                             >  3:117604/1‑130 (MQ=255)
                         aGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcat                                                                                                                                             <  4:117604/132‑3 (MQ=255)
                             aTGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaa                                                                                                                         >  4:147056/1‑148 (MQ=255)
                              tGAAAATACGGCGCAAATTATTGATGCCTGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaa                                                                                                                        <  3:747384/148‑1 (MQ=255)
                              tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                                       >  1:727037/1‑149 (MQ=255)
                              tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                                       >  2:84389/1‑149 (MQ=255)
                              tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaa                                                                                                                        <  4:513364/148‑1 (MQ=255)
                              tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaa                                                                                                                        >  3:320888/1‑148 (MQ=255)
                                  aaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGa                                                                                                                    <  3:680245/148‑1 (MQ=255)
                                  aaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                   <  1:608000/149‑1 (MQ=255)
                                  aaTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                   >  1:860933/1‑149 (MQ=255)
                                   aTACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAcc                                                                                                                  <  2:231578/149‑1 (MQ=255)
                                     aCGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCg                                                                                                                 >  3:813969/1‑148 (MQ=255)
                                      cGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCATGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGa                                                                                                                <  4:598518/148‑1 (MQ=255)
                                      cGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGAt                                                                                                               <  1:136706/149‑1 (MQ=255)
                                           cAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACt                                                                                                          <  2:425744/149‑1 (MQ=255)
                                             aaTTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTg                                                                                                         >  3:148877/1‑148 (MQ=255)
                                              attattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGa                                                                                                       >  1:305805/1‑149 (MQ=255)
                                               ttattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                                                                      <  1:965268/149‑1 (MQ=255)
                                                tattGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATg                                                                                                     <  2:467489/149‑1 (MQ=255)
                                                  ttGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGc                                                                                                    >  3:124479/1‑148 (MQ=255)
                                                  ttGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCg                                                                                                   >  1:684729/1‑149 (MQ=255)
                                                   tGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGa                                                                                                  <  1:427419/149‑1 (MQ=255)
                                                     aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGata                                                                                                >  2:855617/1‑149 (MQ=255)
                                                       gcccgcccGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGatat                                                                                               <  3:527757/148‑1 (MQ=255)
                                                           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                          >  2:680397/1‑149 (MQ=255)
                                                           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                          <  1:211833/149‑1 (MQ=255)
                                                           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                          <  2:731634/149‑1 (MQ=255)
                                                           gcccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                          <  2:1044998/149‑1 (MQ=255)
                                                            cccgGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAtttttt                                                                                          >  4:107071/1‑148 (MQ=255)
                                                                    cACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGaa                                                                                                                                    >  2:379795/1‑98 (MQ=255)
                                                                    cACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGaa                                                                                                                                    <  1:379795/98‑1 (MQ=255)
                                                                    cACGTTTTAACGAAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc                                                                                 <  2:256868/149‑1 (MQ=255)
                                                                      cGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggcg                                                                                >  3:575490/1‑148 (MQ=255)
                                                                         tttAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCa                                                                             >  4:803737/1‑148 (MQ=255)
                                                                           tAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCt                                                                          >  1:779263/1‑149 (MQ=255)
                                                                            aaCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTa                                                                         >  1:941609/1‑149 (MQ=255)
                                                                             aCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTa                                                                         <  4:124479/148‑1 (MQ=255)
                                                                                cAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                                                                      <  4:148877/148‑1 (MQ=255)
                                                                                cAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGa                                                                      >  3:205444/1‑148 (MQ=255)
                                                                                  gAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                   >  1:195115/1‑101 (MQ=255)
                                                                                  gAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGAc                                                                                                                   <  2:195115/101‑1 (MQ=255)
                                                                                   aaGTTGATGAACCTCGCCCACGTTTACGTCGCGGACATATTCCCCGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAGTTTTTGGTCGCGGCGTCAGCTACGACaaa                                                                  <  2:832944/149‑1 (MQ=255)
                                                                                    aGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAAc                                                                 <  2:977351/149‑1 (MQ=255)
                                                                                        gatgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAAt                                                             >  1:213559/1‑149 (MQ=255)
                                                                                          tgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                           <  1:401903/117‑1 (MQ=255)
                                                                                          tgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                           >  2:401903/1‑117 (MQ=255)
                                                                                            aaCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTAt                                                          >  4:765238/1‑148 (MQ=255)
                                                                                            aaCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTAt                                                          <  4:882938/148‑1 (MQ=255)
                                                                                                tCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCa                                                     >  2:345210/1‑149 (MQ=255)
                                                                                                    ccAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTg                                                 >  2:888187/1‑149 (MQ=255)
                                                                                                      aGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGc                                                <  4:960526/148‑1 (MQ=255)
                                                                                                      aGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCg                                               >  2:396301/1‑149 (MQ=255)
                                                                                                         tttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGct                                            >  2:314543/1‑149 (MQ=255)
                                                                                                              gTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTgg                                        <  4:924165/148‑1 (MQ=255)
                                                                                                                cgcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgtg                                      >  4:21493/1‑148 (MQ=255)
                                                                                                                  cgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTa                                    <  4:813969/148‑1 (MQ=255)
                                                                                                                    gACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAAc                                  >  4:335453/1‑148 (MQ=255)
                                                                                                                       atatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgt                                       <  2:633716/140‑1 (MQ=255)
                                                                                                                       atatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGgt                                       >  1:633716/1‑140 (MQ=255)
                                                                                                                       atatTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCa                              >  1:136440/1‑149 (MQ=255)
                                                                                                                          ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGa                                                                                                       <  3:916368/73‑1 (MQ=255)
                                                                                                                          ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGa                                                                                                       >  4:916368/1‑73 (MQ=255)
                                                                                                                          ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                      <  4:893322/122‑1 (MQ=255)
                                                                                                                          ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTc                                                      >  3:893322/1‑122 (MQ=255)
                                                                                                                            cccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGt                         <  1:914808/149‑1 (MQ=255)
                                                                                                                              cGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTgg                       <  2:779263/149‑1 (MQ=255)
                                                                                                                                gAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtt                     >  2:727182/1‑149 (MQ=255)
                                                                                                                                 aGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGttt                    >  1:762449/1‑149 (MQ=255)
                                                                                                                                  gCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtttt                   <  2:755116/149‑1 (MQ=255)
                                                                                                                                    aCTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTg                  <  4:467544/148‑1 (MQ=255)
                                                                                                                                     cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt                >  1:185301/1‑149 (MQ=255)
                                                                                                                                     cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt                <  2:727037/149‑1 (MQ=255)
                                                                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt                <  4:492071/148‑1 (MQ=255)
                                                                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt                <  3:803737/148‑1 (MQ=255)
                                                                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa               >  1:811466/1‑149 (MQ=255)
                                                                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa               >  2:845469/1‑149 (MQ=255)
                                                                                                                                      tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa               <  1:396301/149‑1 (MQ=255)
                                                                                                                                      caccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGtt                >  4:433177/6‑148 (MQ=255)
                                                                                                                                       gAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTa               >  4:774543/1‑148 (MQ=255)
                                                                                                                                       accagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTg                                         <  4:249670/118‑1 (MQ=255)
                                                                                                                                       accagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTg                                         >  3:249670/5‑122 (MQ=255)
                                                                                                                                        ccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGc             >  2:852288/4‑149 (MQ=255)
                                                                                                                                         aTGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGc             <  4:221755/148‑1 (MQ=255)
                                                                                                                                         aTGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGc             >  3:568936/1‑148 (MQ=255)
                                                                                                                                          agTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCa            >  4:1021875/2‑148 (MQ=255)
                                                                                                                                          agTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCa            >  4:9257/2‑148 (MQ=255)
                                                                                                                                            ttCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGTGTCAGCTCCGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTc         >  2:189699/1‑149 (MQ=255)
                                                                                                                                                gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                 <  1:772992/137‑1 (MQ=255)
                                                                                                                                                gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGt                 >  2:772992/1‑137 (MQ=255)
                                                                                                                                                gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGa     >  1:608252/1‑149 (MQ=255)
                                                                                                                                                  ggACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACg   >  1:532941/1‑149 (MQ=255)
                                                                                                                                                   gACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGc  <  1:345210/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGC  >  NC_000913/2652958‑2653253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: