Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,475 G→A 7.4% L61L (TTG→TTA) ‡ uspD → stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4750GA 7.4% 301.1 / 11.4 94L61L (TTG→TTA) ‡uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base A (4/3);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 101 total aligned reads displayed.

AATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGT  >  NC_000913/4113329‑4113622
                                                                                                                                                  |                                                                                                                                                   
aaTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  2:663553/148‑1 (MQ=255)
          aTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATaa                                                                                                                                          >  4:4750/1‑148 (MQ=255)
               ttACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                                                                     >  4:118179/1‑148 (MQ=255)
                        aaTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCTGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                                                                         <  2:463000/135‑1 (MQ=255)
                        aaTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCTGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                                                                         >  1:463000/1‑135 (MQ=255)
                         aTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCt                                                                                                                           >  3:55883/1‑148 (MQ=255)
                              gCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                                                                      >  2:117836/1‑148 (MQ=255)
                               cccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                                                                     <  4:55883/148‑1 (MQ=255)
                                  tGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                                                                  >  3:364918/1‑148 (MQ=255)
                                   ggAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                                                                 <  3:4750/148‑1 (MQ=255)
                                     aGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                                                                               >  1:665812/1‑148 (MQ=255)
                                          gCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                                                                          <  1:117836/148‑1 (MQ=255)
                                           ccAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                                                                           <  1:839750/114‑1 (MQ=255)
                                           ccAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                                                                           >  2:839750/1‑114 (MQ=255)
                                               aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                                                                     <  4:364918/148‑1 (MQ=255)
                                                cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                                                                     >  4:1013660/1‑131 (MQ=255)
                                                cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                                                                     <  3:1013660/131‑1 (MQ=255)
                                                   aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAAt                                                                                                 <  3:118179/148‑1 (MQ=255)
                                                     tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    >  4:322022/1‑95 (MQ=255)
                                                     tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                    <  3:322022/95‑1 (MQ=255)
                                                     tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                                                                               >  3:125629/1‑148 (MQ=255)
                                                          cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                       <  4:1096561/103‑1 (MQ=255)
                                                          cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                       >  3:1096561/1‑103 (MQ=255)
                                                          cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTgaagaa                                                                                                                                             >  4:829063/1‑97 (MQ=255)
                                                            cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                                                                          >  1:1035396/1‑130 (MQ=255)
                                                            cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                                                                          <  2:1035396/130‑1 (MQ=255)
                                                                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                                                                <  1:240109/136‑1 (MQ=255)
                                                                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                                                                >  2:240109/1‑136 (MQ=255)
                                                                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                                                                    <  2:943156/148‑1 (MQ=255)
                                                                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                                                                    >  1:943156/1‑148 (MQ=255)
                                                                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                                                                 <  1:873037/98‑1 (MQ=255)
                                                                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                                                                 >  2:873037/1‑98 (MQ=255)
                                                                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                                                                  <  1:937842/129‑1 (MQ=255)
                                                                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                                                                  >  2:937842/1‑129 (MQ=255)
                                                                         ttCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                                                                           <  4:125629/148‑1 (MQ=255)
                                                                          tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                                                                <  2:177218/142‑1 (MQ=255)
                                                                          tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                                                                >  1:177218/1‑142 (MQ=255)
                                                                          tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                                                                          >  4:417917/1‑148 (MQ=255)
                                                                              atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                                                                  >  4:786416/1‑72 (MQ=255)
                                                                              atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                                                                  <  3:786416/72‑1 (MQ=255)
                                                                                tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                                                                        <  4:1051892/128‑1 (MQ=255)
                                                                                tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                                                                    >  4:400869/1‑148 (MQ=255)
                                                                                tgatgaTAGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                                                                        >  3:1051892/1‑128 (MQ=255)
                                                                                  atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                       >  3:403580/1‑79 (MQ=255)
                                                                                  atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                       <  4:403580/79‑1 (MQ=255)
                                                                                   tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcnnn                                                                 <  2:665812/148‑4 (MQ=255)
                                                                                      tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                                                                       >  4:1073226/1‑139 (MQ=255)
                                                                                      tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                                                                       <  3:1073226/139‑1 (MQ=255)
                                                                                       ggCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa                                                             <  3:400869/148‑1 (MQ=255)
                                                                                          ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                                                                            >  4:226604/1‑82 (MQ=255)
                                                                                          ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                                                                            <  3:226604/82‑1 (MQ=255)
                                                                                          ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                                                                    <  3:493480/90‑1 (MQ=255)
                                                                                          ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                                                                    >  4:493480/1‑90 (MQ=255)
                                                                                                    tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGTTaa                                                                                                                                 >  3:289898/1‑67 (MQ=255)
                                                                                                    tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGTTaa                                                                                                                                 <  4:289898/67‑1 (MQ=255)
                                                                                                    tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                                                                          >  4:584157/1‑90 (MQ=255)
                                                                                                    tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                                                                          <  3:584157/90‑1 (MQ=255)
                                                                                                      tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgc                                              <  3:417917/148‑1 (MQ=255)
                                                                                                        tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                                                                >  1:974298/1‑112 (MQ=255)
                                                                                                        tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                                                                <  2:974298/112‑1 (MQ=255)
                                                                                                                 cTACTTCCCTGCAACAGATGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                                                                     <  4:716663/114‑1 (MQ=255)
                                                                                                                 cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                                                                     >  3:716663/1‑114 (MQ=255)
                                                                                                                   aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                                                                           <  2:30151/106‑1 (MQ=255)
                                                                                                                   aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                                                                           >  1:30151/1‑106 (MQ=255)
                                                                                                                   aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCa                                 >  4:410563/1‑148 (MQ=255)
                                                                                                                     ttCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaaa                               >  1:463002/1‑148 (MQ=255)
                                                                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAAtt                                      >  2:345227/1‑140 (MQ=255)
                                                                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAAtt                                      <  1:345227/140‑1 (MQ=255)
                                                                                                                        ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAag                            >  4:146705/1‑148 (MQ=255)
                                                                                                                          tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaaa                               >  2:563935/1‑143 (MQ=255)
                                                                                                                          tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaaa                               <  1:563935/143‑1 (MQ=255)
                                                                                                                           gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa                                                     <  3:814004/120‑1 (MQ=255)
                                                                                                                           gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa                                                     >  4:814004/1‑120 (MQ=255)
                                                                                                                            cAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgc                                              <  3:268743/126‑1 (MQ=255)
                                                                                                                            cAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgc                                              >  4:268743/1‑126 (MQ=255)
                                                                                                                              aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAAtt                                      <  4:754080/132‑1 (MQ=255)
                                                                                                                              aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAAtt                                      >  3:754080/1‑132 (MQ=255)
                                                                                                                              aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGt                      <  3:410563/148‑1 (MQ=255)
                                                                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                                                                            >  3:971233/1‑42 (MQ=255)
                                                                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGaa                                         <  1:86136/125‑1 (MQ=255)
                                                                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGaa                                         >  2:86136/1‑125 (MQ=255)
                                                                                                                                      tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaaa                               >  3:1044081/1‑131 (MQ=255)
                                                                                                                                      tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCaaa                               <  4:1044081/131‑1 (MQ=255)
                                                                                                                                        ttcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                                                                   <  3:355748/61‑1 (MQ=255)
                                                                                                                                        ttcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                                                                   >  4:355748/1‑61 (MQ=255)
                                                                                                                                        ttcttcAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAg                       >  1:980859/1‑137 (MQ=255)
                                                                                                                                        ttcttcAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAg                       <  2:980859/137‑1 (MQ=255)
                                                                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa                                                             >  3:619705/1‑95 (MQ=255)
                                                                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa                                                             <  4:619705/95‑1 (MQ=255)
                                                                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa                                                     <  1:1020375/103‑1 (MQ=255)
                                                                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa                                                     >  2:1020375/1‑103 (MQ=255)
                                                                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc        >  4:772033/1‑148 (MQ=255)
                                                                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc        <  4:1029598/147‑1 (MQ=255)
                                                                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc        >  3:1029598/1‑147 (MQ=255)
                                                                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAAAGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCt       >  4:708378/1‑148 (MQ=255)
                                                                                                                                               aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActg                                               >  1:33399/1‑106 (MQ=255)
                                                                                                                                               aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActg                                               <  2:33399/106‑1 (MQ=255)
                                                                                                                                               aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCtgt     <  3:146705/148‑1 (MQ=255)
                                                                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGt  >  4:1041633/1‑148 (MQ=255)
                                                                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGt  <  3:772033/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGT  >  NC_000913/4113329‑4113622

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