Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,709,897 G→A 5.3% L252L (CTG→CTA) ‡ rsxD → SoxR [2Fe‑2S] reducing system protein RsxD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,709,8970GA 5.3% 239.6 / 2.2 75L252L (CTG→CTA) ‡rsxDSoxR [2Fe‑2S] reducing system protein RsxD
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base A (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

ATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGGTCGTCTTA  >  NC_000913/1709754‑1710045
                                                                                                                                               |                                                                                                                                                    
aTTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTa                                                                                                                                                <  2:377670/150‑1 (MQ=255)
  tATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATg                                                                                                                                              <  2:1429192/150‑1 (MQ=255)
   aTGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCCAGCGTTacta                                                                                                                                             <  1:1370863/150‑4 (MQ=255)
    tGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgcg                                                                                                                                            <  2:843879/150‑1 (MQ=255)
     gCAATATCCGATCTACTGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgcgc                                                                                                                                           <  1:892161/150‑1 (MQ=255)
     gCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATgcgc                                                                                                                                           <  2:1263458/150‑1 (MQ=255)
      cAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCa                                                                                                                                          >  1:1097993/1‑150 (MQ=255)
         tatCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATg                                                                                                                                       >  1:996952/1‑150 (MQ=255)
          atCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGt                                                                                                                                      <  2:542123/150‑1 (MQ=255)
            ccGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTg                                                                                                                                    <  1:1012252/150‑1 (MQ=255)
             cGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTgg                                                                                                                                   >  2:1154395/1‑150 (MQ=255)
                   aCAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgt                                                                                                                             <  1:155516/150‑1 (MQ=255)
                   aCAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgt                                                                                                                             <  1:264964/150‑1 (MQ=255)
                    aaGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgtt                                                                                                                            <  2:418296/149‑1 (MQ=255)
                     aGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGgttg                                                                                                                           <  2:80181/150‑1 (MQ=255)
                            ttCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCAc                                                                                                                    >  2:620955/1‑150 (MQ=255)
                             tCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCAcc                                                                                                                   <  1:1543861/150‑1 (MQ=255)
                                    ggcgCTGGTTGGCATTGGGTAAATCTCGCCTGGCTCGCGGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTATTAAAGATGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAacc                                                                                                            >  2:491503/1‑149 (MQ=255)
                                       gcTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCCAGCGTTacta                                                                                                                                             <  2:578209/114‑4 (MQ=255)
                                       gcTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCCAGCGTTacta                                                                                                                                             >  1:578209/1‑111 (MQ=255)
                                         gggTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTTGCATATTCCCCTCATCTTCTTAGTAACGCTGGGTTTATGCGCAATGTTTGGCTGGTTGTTGTCATCAGAAACAGTggg                                                                                                       <  2:399278/149‑2 (MQ=255)
                                          ggTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGGTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGgca                                                                                                      <  2:1250357/150‑1 (MQ=255)
                                           gTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTTCTATGTCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGgcag                                                                                                     <  2:1389214/150‑1 (MQ=255)
                                              ggCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCAc                                                                                                  >  1:843089/1‑150 (MQ=255)
                                                     ggTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGAACCGCAAAt                                                                                           >  1:3349/1‑150 (MQ=255)
                                                        aaaTCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTggg                                                                                                                                  >  1:1671917/1‑108 (MQ=255)
                                                        aaaTCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTggg                                                                                                                                  <  2:1671917/108‑1 (MQ=255)
                                                               gCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCCTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATctgctg                                                                                 <  2:602602/150‑1 (MQ=255)
                                                                ccTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGt                                                                                >  1:568348/1‑150 (MQ=255)
                                                                      tggctggcGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAg                                                                          >  1:687978/1‑150 (MQ=255)
                                                                       ggctggcGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGc                                                                         <  2:137434/150‑1 (MQ=255)
                                                                         ctggcGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGa                                                                       >  2:1147727/1‑150 (MQ=255)
                                                                           ggcggcGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGAcc                                                                     <  1:161978/150‑1 (MQ=255)
                                                                            gcggcgTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCa                                                                    <  2:807093/150‑1 (MQ=255)
                                                                             cggcgTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCAt                                                                   <  2:1308064/150‑1 (MQ=255)
                                                                               gcgTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGc                                                                 >  1:1337350/1‑150 (MQ=255)
                                                                                    tGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGgc                                                            <  2:128525/150‑1 (MQ=255)
                                                                                      gTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGgcgc                                                          >  2:1330059/1‑150 (MQ=255)
                                                                                          cTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACTGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTc                                                      <  2:366040/150‑1 (MQ=255)
                                                                                           tATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCt                                                     <  1:136337/150‑1 (MQ=255)
                                                                                            aTGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGgcgc                                                          <  2:1632946/144‑1 (MQ=255)
                                                                                            aTGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGgcgc                                                          >  1:1632946/1‑144 (MQ=255)
                                                                                               gCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTAt                                                 >  1:930951/1‑150 (MQ=255)
                                                                                                  gAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTAtttt                                              <  2:177231/150‑1 (MQ=255)
                                                                                                  gAAAGCGATTCGCTGGCACATTCCCCTCTGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGGCTGGAGCGACCATGCTCGGCGCATTCTTTAtttt                                              <  2:948591/150‑1 (MQ=255)
                                                                                                      gCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTtgact                                          >  2:7472/1‑150 (MQ=255)
                                                                                                       cGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTtgactg                                         >  1:994616/1‑150 (MQ=255)
                                                                                                        gATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTtgactga                                        >  1:388912/1‑150 (MQ=255)
                                                                                                          ttCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGAcc                                      >  2:446826/1‑150 (MQ=255)
                                                                                                             gCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCgg                                   >  2:683245/1‑150 (MQ=255)
                                                                                                              cTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGt                                  >  1:37145/1‑150 (MQ=255)
                                                                                                                  cATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTAcc                              >  2:745461/1‑150 (MQ=255)
                                                                                                                     atTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCt                           <  1:986858/150‑1 (MQ=255)
                                                                                                                        ccccTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCt                        <  1:1381/150‑1 (MQ=255)
                                                                                                                         cccTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       >  1:481054/1‑150 (MQ=255)
                                                                                                                         cccTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       >  2:385116/1‑150 (MQ=255)
                                                                                                                           cTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACg                     >  2:844665/1‑150 (MQ=255)
                                                                                                                            tCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGa                    >  1:1403409/1‑150 (MQ=255)
                                                                                                                             cAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGAc                   >  2:565736/1‑150 (MQ=255)
                                                                                                                              aggttcttAGTATCGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCAGCAGAAACACTGGCAGCACCGCTAATTCATCTGCTGTCTGGAGCGACCTTGCTCGTCGCATTCTTTATTTTGACTGACCCGGTTACCTCTTCTACGAcc                  <  1:896606/147‑1 (MQ=255)
                                                                                                                              aGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGAcc                  >  2:510271/1‑150 (MQ=255)
                                                                                                                                   gcgcatAACGCCAGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCg             >  1:1176101/6‑150 (MQ=255)
                                                                                                                                   cttAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCg             <  2:1097993/150‑1 (MQ=255)
                                                                                                                                   cttAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCg             >  1:1291681/1‑150 (MQ=255)
                                                                                                                                   cttAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACACGGTTACAGCTTCTACGACCAATCg             >  2:415530/1‑150 (MQ=255)
                                                                                                                                    ttAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGt            <  1:1147727/150‑1 (MQ=255)
                                                                                                                                      aGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTgg          <  2:930951/150‑1 (MQ=255)
                                                                                                                                      aGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCAACTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTgg          >  1:1466594/1‑150 (MQ=255)
                                                                                                                                       gTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGgt         <  2:996952/150‑1 (MQ=255)
                                                                                                                                        tAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       >  2:1089222/1‑135 (MQ=255)
                                                                                                                                        tAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       <  1:1089222/135‑1 (MQ=255)
                                                                                                                                           cGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGgtcgtc     >  1:544022/1‑150 (MQ=255)
                                                                                                                                            gCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGGTCGTCt    <  2:1337350/150‑1 (MQ=255)
                                                                                                                                            gCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCAACTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGGTCGTCt    <  2:1466594/150‑1 (MQ=255)
                                                                                                                                              tGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGGTCGTCTTa  <  1:432274/150‑1 (MQ=255)
                                                                                                                                               ggCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       <  2:332431/128‑1 (MQ=255)
                                                                                                                                               ggCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTa                       >  1:332431/1‑128 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ATTATGCAATATCCGATCTACAGCGGTATTCTGGCGGGCGCTGGTTGGCAATGGGTAAATCTCGCCTGGCTGGCTGGCGGCGTATGGTTGCTATGGCAGAAAGCGATTCGCTGGCATATTCCCCTCAGCTTCTTAGTAACGCTGGCGTTATGCGCAATGTTGGGCTGGTTGTTCTCACCAGAAACACTGGCAGCACCGCAAATTCATCTGCTGTCTGGAGCGACCATGCTCGGCGCATTCTTTATTTTGACTGACCCGGTTACCGCTTCTACGACCAATCGTGGTCGTCTTA  >  NC_000913/1709754‑1710045

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