Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,103 C→T 15.2% H112Y (CAT→TAT)  hyaD → hydrogenase 1 maturation protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,1030CT15.2% 440.5 / 58.3 177H112Y (CAT→TAT) hyaDhydrogenase 1 maturation protease
Reads supporting (aligned to +/- strand):  ref base C (82/68);  new base T (10/17);  total (92/85)
Fisher's exact test for biased strand distribution p-value = 9.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 180 total aligned reads displayed.

TATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCT  >  NC_000913/1036004‑1036201
                                                                                                   |                                                                                                  
tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGAc                                                                                                    <  1:1418187/100‑1 (MQ=255)
  tGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACAt                                                                                                  >  1:1013736/1‑100 (MQ=255)
  tGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGGGTTGGCGCGCGCGGATATCCGCGGAc                                                                                                    >  1:2696818/1‑98 (MQ=255)
   gCCGGAGAACGCATTCCGGCTTATCTCCGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACAt                                                                                                  >  1:1146011/1‑99 (MQ=255)
   gCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATc                                                                                                 >  2:3595351/1‑100 (MQ=255)
   gCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATc                                                                                                 <  2:1996461/100‑1 (MQ=255)
          aaCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGc                                                                                          <  2:3269500/100‑1 (MQ=255)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGca                                                                                         >  2:3370141/1‑100 (MQ=255)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGAGGACATCTGCCAGca                                                                                         >  2:4327368/1‑100 (MQ=255)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGTTGTCCGCGGATATCCGCCAGcg                                                                                         <  1:3376332/100‑2 (MQ=38)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  1:1170664/100‑2 (MQ=38)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  2:4294943/100‑2 (MQ=38)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  2:1828099/100‑2 (MQ=38)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  2:2768268/100‑2 (MQ=38)
           aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  2:562160/100‑2 (MQ=38)
            cGCATTCCGGCTTATCTCAGCGCGACGAAAATGAGCCTGCATCATAACAGTGTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGc                                                                                          <  1:536499/98‑1 (MQ=37)
            cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcac                                                                                        >  2:2195191/1‑100 (MQ=255)
            cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                         <  1:1450707/99‑2 (MQ=38)
             gCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcaca                                                                                       >  2:826484/1‑100 (MQ=255)
              cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                                      >  1:3425753/1‑100 (MQ=255)
              cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgcca                                                                                      <  1:2823199/100‑5 (MQ=37)
              cATTCCGGCGTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                                      >  1:3223630/1‑100 (MQ=255)
               aTTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACat                                                                                      >  1:3355602/1‑99 (MQ=255)
                ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgccaac                                                                                    <  1:1555405/100‑7 (MQ=25)
                ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgccaac                                                                                    <  1:379437/100‑7 (MQ=37)
                   cGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAACTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                                 <  1:2443314/100‑1 (MQ=255)
                   cGGCTTATCTCAGCGCGAAGAAAATGAGCCTGAATCAGAACAGTTTCTCCGAAGTGTTTGCGCTGGCGGATATCCGCGGACCTCTGCCAGCACATATTg                                                                                  >  1:2331984/1‑99 (MQ=255)
                     gCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGccc                                                                               <  1:2798479/100‑1 (MQ=255)
                     gCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGccc                                                                               <  2:652754/100‑1 (MQ=255)
                        tATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcg                                                                            <  2:2072716/100‑1 (MQ=255)
                         aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgt                                                                           <  2:1138688/100‑1 (MQ=255)
                           ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                         >  2:3667261/1‑100 (MQ=255)
                           ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                         >  1:4332685/1‑100 (MQ=255)
                           ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCtcgtcg                                                                         >  2:2420353/1‑100 (MQ=255)
                             cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGt                                                                       >  2:745226/1‑100 (MQ=255)
                             cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGt                                                                       >  2:626672/1‑100 (MQ=255)
                              aGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                                      <  1:3022896/100‑1 (MQ=255)
                              aGCGCGAAGAAAATGAGCCTGCATCACAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                                      >  2:2680862/1‑100 (MQ=255)
                               gcgcGAAGCAAATGAGCCTGCATCCGAACAGTTTCTCCGGAGTGTTGGCGCTGGCGGGTATCCGCGGACAGCTGCCAGCACATGTTGCCCccgtctggct                                                                     >  1:1576231/1‑95 (MQ=255)
                                cgcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTg                                                                    <  2:1219066/100‑1 (MQ=255)
                                cgcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTg                                                                    <  1:149789/100‑1 (MQ=255)
                                 gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                   <  2:4269215/100‑1 (MQ=255)
                                 gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                   >  1:189521/1‑100 (MQ=255)
                                 gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                   >  1:951366/1‑100 (MQ=255)
                                  cgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCa                                                                  <  2:1013740/100‑1 (MQ=255)
                                  cgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCa                                                                  <  2:2567547/100‑1 (MQ=255)
                                        aaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGc                                                            >  2:3306502/1‑100 (MQ=255)
                                              gCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCt                                                      <  2:3271561/100‑1 (MQ=255)
                                              gCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCt                                                      >  1:2853986/1‑100 (MQ=255)
                                                cTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                                    <  1:2594342/100‑1 (MQ=255)
                                                 tGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcga                                                   <  2:854823/100‑1 (MQ=255)
                                                 tGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcga                                                   >  2:1939520/1‑100 (MQ=255)
                                                 tGCATCAGAACAGTTTCCCCGACGTGTTGGCGCTGGCGGATATCGGCGGACATCTGCCAGCACACATTCCCCTCGTCGGCCCGCAACCCGCAATGCTc                                                     >  1:2284719/1‑98 (MQ=255)
                                                  gCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                  <  2:3685682/100‑1 (MQ=255)
                                                     tCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCgacgac                                               >  2:1905570/1‑100 (MQ=255)
                                                       aGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTa                                             >  2:3322476/1‑100 (MQ=255)
                                                         aactgcttCTCCGCAGTGTGGGAGCTGGGGGATGTCCGCAGGCAGCTGCCAGCGAATATTGCCCGAGTCGGTGTGCAACCCGCAATGCTc                                                     <  1:348765/84‑1 (MQ=255)
                                                          aCAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcgg                                          <  1:4270883/100‑1 (MQ=255)
                                                             gTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                       <  1:2042705/100‑1 (MQ=255)
                                                               ttCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTa                                     >  1:1495812/1‑100 (MQ=255)
                                                                 ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGc                                   >  1:2942150/1‑100 (MQ=255)
                                                                 ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGc                                   >  1:3411123/1‑100 (MQ=255)
                                                                 ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGc                                   >  1:3568192/1‑100 (MQ=255)
                                                                 ctcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGACTGCAACCCGCAATGCTCGACGACTAcggcgg                                       >  1:3688833/1‑96 (MQ=255)
                                                                   ccGGAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATACTGCCCTCGTCGGTCTGCAACCCGCAATGCTGGACGACTACGGCGGTAGCCt                                 >  1:3968963/1‑100 (MQ=255)
                                                                   ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCt                                 >  2:4071112/1‑100 (MQ=255)
                                                                    cGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                <  2:1005090/100‑1 (MQ=255)
                                                                    cGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCGGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTg                                >  1:882155/1‑100 (MQ=255)
                                                                    cGAAGTGGTGGCGCTGACGAATATCAAAAGACAACTGCCAGCACATATTGACCTCGTAGATCTGCGACCAGCAAAGCCCGACGACTCCGGCGGTAGACt                                 >  1:2313280/1‑99 (MQ=255)
                                                                       gcaatatgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGc                             >  2:2526760/10‑100 (MQ=21)
                                                                        gtgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCg                            >  2:593814/1‑100 (MQ=255)
                                                                         tgtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCCGTCTGCAACCCGCAACGCTCGACGACTACGGCGGTAGCCTGAGCGa                           <  1:3508102/100‑1 (MQ=255)
                                                                          gtTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGaa                          >  1:3683096/1‑100 (MQ=255)
                                                                           ttGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                         >  1:2127926/1‑100 (MQ=255)
                                                                            tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                         <  1:530350/99‑1 (MQ=255)
                                                                            tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                        <  2:973511/100‑1 (MQ=255)
                                                                            tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                        >  1:2575869/1‑100 (MQ=255)
                                                                             tgtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTg                       >  1:585750/4‑100 (MQ=37)
                                                                              gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg                      >  2:831174/3‑100 (MQ=37)
                                                                              gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg                      >  2:4138016/3‑100 (MQ=37)
                                                                               tgcTGGCAGATGTCCGCGGATATCCGCCAGAACGTATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGGGCGAACTGGc                     >  1:2891110/2‑100 (MQ=25)
                                                                               cgcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc                     >  2:1599392/1‑100 (MQ=255)
                                                                               cgcTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc                     >  2:3028402/1‑100 (MQ=255)
                                                                                gcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCa                    >  1:1835026/1‑100 (MQ=38)
                                                                                 cTGGCTGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCa                    <  1:3850158/99‑1 (MQ=255)
                                                                                      ggATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg              >  2:1179209/1‑100 (MQ=255)
                                                                                      ggATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg              <  1:3480168/100‑1 (MQ=255)
                                                                                          atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAc          <  1:2557552/100‑1 (MQ=255)
                                                                                          atCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAc          <  2:2317773/100‑1 (MQ=255)
                                                                                           tCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACt         >  2:4114149/1‑100 (MQ=255)
                                                                                           tCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACt         >  2:717556/1‑100 (MQ=255)
                                                                                            ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCGACCCGCAATGCTCGACGACAACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg        >  1:2390214/1‑100 (MQ=255)
                                                                                            ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg        >  2:1279215/1‑100 (MQ=255)
                                                                                            ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg        <  2:1645875/100‑1 (MQ=255)
                                                                                            ccGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTg        <  1:1536311/100‑1 (MQ=255)
                                                                                             cgcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGc       >  1:3918481/1‑100 (MQ=255)
                                                                                              gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGcc      <  2:2608485/100‑1 (MQ=255)
                                                                                               cgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGccc     <  2:3018060/100‑1 (MQ=255)
                                                                                                ggACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCg    >  2:3871715/1‑100 (MQ=255)
                                                                                                  aCATCTGCAAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTggc      >  1:4164862/1‑94 (MQ=255)
                                                                                                   |                                                                                                  
TATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCT  >  NC_000913/1036004‑1036201

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