Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,672 G→C 8.0% intergenic (+20/+43) hslO → / ← yhgE heat shock protein Hsp33/DUF4153 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,6720GC 8.0% 292.5 / 14.6 87intergenic (+20/+43)hslO/yhgEheat shock protein Hsp33/DUF4153 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base C (4/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAG  >  NC_000913/3530525‑3530813
                                                                                                                                                   |                                                                                                                                             
gCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCagagag                                                                                                                                               <  3:303358/148‑1 (MQ=255)
      tGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTg                                                                                                                                        >  5:111901/1‑149 (MQ=255)
           gAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGACGGCTCTGCCGGa                                                                                                                                   >  7:258119/1‑149 (MQ=255)
           gAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGACGGCTCTGCCGGa                                                                                                                                   >  7:322789/1‑149 (MQ=255)
            aaGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCGGa                                                                                                                                   >  4:305846/1‑148 (MQ=255)
            aaGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCGGa                                                                                                                                   <  3:305846/148‑1 (MQ=255)
             aGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGAtt                                                                                                                                 <  5:380355/149‑1 (MQ=255)
              gATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGAttt                                                                                                                                <  3:303968/149‑1 (MQ=255)
               aTGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGAtttt                                                                                                                               <  4:140323/149‑1 (MQ=255)
                tGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTa                                                                                                                              >  1:393267/1‑149 (MQ=255)
                  gagAAATTGACATGCATTGTGATTACTGCGGTAACCACTTTCTGTTCATTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                             <  6:163119/146‑1 (MQ=255)
                  gCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                             <  7:82439/148‑1 (MQ=255)
                  gCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                             >  1:397448/1‑148 (MQ=255)
                   cGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGa                                                                                                                           <  8:203448/149‑1 (MQ=255)
                     aaaTTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGAc                                                                                                                          >  8:159360/1‑148 (MQ=255)
                       aTTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTaa                                                                                                                       <  6:111901/149‑1 (MQ=255)
                          gACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                     <  2:136302/148‑1 (MQ=255)
                          gACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                     <  2:393267/148‑1 (MQ=255)
                            cATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGACGGCTCTGCCGGATTTTACGACTAATacac                                                                                                                  <  8:258119/149‑1 (MQ=255)
                            cATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGACGGCTCTGCCGGATTTTACGACTAATacac                                                                                                                  <  8:322789/149‑1 (MQ=255)
                             aTGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTa                                                                                                                              >  7:318170/1‑136 (MQ=255)
                             aTGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTa                                                                                                                              <  8:318170/136‑1 (MQ=255)
                             aTGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACAcc                                                                                                                 >  1:63291/1‑149 (MQ=255)
                              tGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCt                                                                                                                >  4:184631/1‑149 (MQ=255)
                               gCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCt                                                                                                                >  6:20139/1‑148 (MQ=255)
                                 aTTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAAAGAGTATGTCCGCCAGAGAGCGCTCTGCCGGATATAACGACTAATACACCTCAt                                                                                                             >  4:36493/1‑149 (MQ=255)
                                  ttGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATg                                                                                                            >  8:53622/1‑149 (MQ=255)
                                  ttGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATg                                                                                                            >  4:197410/1‑149 (MQ=255)
                                     tgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgcgc                                                                                                         >  2:69088/1‑149 (MQ=255)
                                     tgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgcgc                                                                                                         >  8:348154/1‑149 (MQ=255)
                                           cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                   <  3:184631/149‑1 (MQ=255)
                                            tGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCa                                                                                                  >  2:369918/1‑149 (MQ=255)
                                              cGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACt                                                                                                <  7:159360/149‑1 (MQ=255)
                                                   aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTc                                                                                           <  1:369918/149‑1 (MQ=255)
                                                   aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTc                                                                                           <  1:256672/149‑1 (MQ=255)
                                                     cACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                            >  2:407706/1‑146 (MQ=255)
                                                     cACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                            >  2:144181/1‑146 (MQ=255)
                                                     cACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                            >  5:336655/1‑146 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGAc                                                                                                                          >  7:402281/1‑115 (MQ=255)
                                                      aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGAc                                                                                                                          <  8:402281/115‑1 (MQ=255)
                                                          tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTca                                                                                    >  6:161707/1‑149 (MQ=255)
                                                          tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTca                                                                                    >  2:385104/1‑149 (MQ=255)
                                                                 aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTAATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAg                                                                             >  2:116975/1‑149 (MQ=255)
                                                                      gATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCttt                                                                        <  1:116975/149‑1 (MQ=255)
                                                                       aTGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                            <  4:169390/128‑1 (MQ=255)
                                                                       aTGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCAtt                                                                                            >  3:169390/1‑128 (MQ=255)
                                                                        tGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCa                                                                      >  1:44265/1‑149 (MQ=255)
                                                                        tGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCAGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCa                                                                      <  1:60884/149‑1 (MQ=255)
                                                                         ggATATTGCTGAAATCCGCTACAACGCGTCTCCGGCAGTTCGGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCGGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCa                                                                      <  3:36493/148‑1 (MQ=255)
                                                                          gttatTGCTGATCTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTATTGAGTTTGGCCGGCAGAGAGGGCGCTGCCGGATTTTACGACTAATACCCCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAc                                                                     <  3:51546/146‑1 (MQ=255)
                                                                          gATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCAc                                                                     <  2:44265/148‑1 (MQ=255)
                                                                           atatTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCa                                                                   <  1:407706/149‑1 (MQ=255)
                                                                              ttGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATc                                                                 <  1:144181/148‑1 (MQ=255)
                                                                              ttGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATc                                                                 <  2:63291/148‑1 (MQ=255)
                                                                                 cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCAc                                                              <  5:20139/148‑1 (MQ=255)
                                                                                  tGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTa                                                            >  2:202366/1‑149 (MQ=255)
                                                                                   gAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTa                                                            >  3:291435/1‑148 (MQ=255)
                                                                                    aaaTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACa                                                          <  6:336655/149‑1 (MQ=255)
                                                                                    aaaTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCGTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACa                                                          <  7:348154/149‑1 (MQ=255)
                                                                                            aacaacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTc                                                                       >  1:159585/1‑128 (MQ=255)
                                                                                            aacaacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTc                                                                       <  2:159585/128‑1 (MQ=255)
                                                                                              caacGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGcc                                                >  5:133362/1‑149 (MQ=255)
                                                                                              caacGCGTCTCCGGCAGATCCGCAAGTTCATTAAGGAGTATGTCCGGCAGAGCGGGCTCTGCCGGATTTTACGACTAATACCCGTGATGCGCTTTGTCCCCCATTCAGATGCCCATTCAGCCTTTCACCATCCACTACAAGATCgcgc                                                 >  6:353538/1‑148 (MQ=255)
                                                                                                  gcgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGc                                            <  8:399916/149‑1 (MQ=255)
                                                                                                  gcgTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGc                                            <  4:291435/149‑1 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTAACATTCCGCCTTTCACCATCCACTACAAGATCGCGCCACGCtt                                          >  7:158955/1‑147 (MQ=255)
                                                                                                      ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTa                                        <  6:133362/149‑1 (MQ=255)
                                                                                                            cAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTc                                  <  5:161707/149‑1 (MQ=255)
                                                                                                                   gCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCATATTCACATTCAGCCTTTGACCATCCCCTACAAGATCGCGCCACGATTTAGGCTTCGTccc                             >  6:55748/1‑147 (MQ=255)
                                                                                                                           aTTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTc                                                                                           >  1:260096/1‑77 (MQ=255)
                                                                                                                           aTTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTc                                                                                           <  2:260096/77‑1 (MQ=255)
                                                                                                                            ttAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCAt                  <  7:53622/149‑1 (MQ=255)
                                                                                                                             tAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATc                 >  8:35764/1‑149 (MQ=255)
                                                                                                                               aTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATcc                <  1:202366/148‑1 (MQ=255)
                                                                                                                               aTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt               <  3:299013/149‑1 (MQ=255)
                                                                                                                               aTGAGTATGTCCGGCAGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTTAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt               <  7:247632/149‑1 (MQ=255)
                                                                                                                                tGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACaa                                                         >  6:59901/1‑106 (MQ=255)
                                                                                                                                tGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACaa                                                         <  5:59901/106‑1 (MQ=255)
                                                                                                                                 gAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCCGTTTGCCATCCTTc             >  5:251965/1‑149 (MQ=255)
                                                                                                                                   gTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCa                                                               >  5:118048/1‑97 (MQ=255)
                                                                                                                                   gTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCa                                                               <  6:118048/97‑1 (MQ=255)
                                                                                                                                    tATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCa          <  1:385104/149‑1 (MQ=255)
                                                                                                                                       gTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt       >  6:286429/1‑149 (MQ=255)
                                                                                                                                        tCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg      >  8:239506/1‑149 (MQ=255)
                                                                                                                                          cGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgt     >  4:364977/1‑148 (MQ=255)
                                                                                                                                           ggCAGAGAGGGCTGTGCCGGATTTTAGGACTAATACACCTGATTCGCTTTGTCACTCATTCAGTTTCACATTCAGCCTTTCACCATCCACTAGAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCttt     <  7:97276/147‑3 (MQ=255)
                                                                                                                                           ggCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGATTAGACCACAGTTTTCCATCCAACGCACCtgt     >  8:97276/1‑147 (MQ=255)
                                                                                                                                             cAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGTTTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAg  >  1:403389/1‑148 (MQ=255)
                                                                                                                                                  agGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAg         >  8:257104/1‑136 (MQ=255)
                                                                                                                                                  agGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAg         <  7:257104/136‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAG  >  NC_000913/3530525‑3530813

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