Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 931,952 C→T 28.9% V33V (GTG→GTA trxB ← thioredoxin reductase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913931,9520CT28.9% 130.4 / 55.0 85V33V (GTG→GTAtrxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (8/17);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATTAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCA  >  NC_000913/931854‑932049
                                                                                                  |                                                                                                 
aTTAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGc                                                                                                   <  2:1408020/99‑1 (MQ=255)
 ttAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtac                                                                                                 >  2:716601/1‑100 (MQ=255)
  taataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtacn                                                                                                <  2:510275/100‑2 (MQ=255)
   aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcag                                                                                               >  2:49250/1‑100 (MQ=255)
       aCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttg                                                                                           >  2:609273/1‑100 (MQ=255)
        cGGACCGGTCAGATCGTTTGGATCGCCAGGCTAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgc                                                                                          <  1:677559/100‑1 (MQ=255)
        cGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgc                                                                                          <  2:1137095/100‑1 (MQ=255)
          gACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcag                                                                                        >  1:463738/1‑100 (MQ=255)
           aCCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcagg                                                                                       <  1:680472/100‑1 (MQ=255)
            ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                      >  1:474193/1‑100 (MQ=255)
             cGGTCAGATCGTTTGGATCTCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcaggtt                                                                                     <  2:2319388/100‑1 (MQ=255)
             cGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcaggtt                                                                                     >  2:679364/1‑100 (MQ=255)
               gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTgg                                                                                   <  1:906174/100‑1 (MQ=255)
               gTCAGATCGTATGGACCCCCAGGCCAGTTTTCCACTTCCGTGGTCGTCGTCCGCTGGCCTCCTTTTTCCATGCCGGTTATCAGTAcaggtagcaggttt                                                                                    <  2:1148107/99‑2 (MQ=255)
                tCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgc                                                                                  <  2:1176318/100‑1 (MQ=255)
                   gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcg                                                                               >  1:1186809/1‑100 (MQ=255)
                   gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcg                                                                               >  1:818410/1‑100 (MQ=255)
                   gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCAATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcg                                                                               >  1:1921174/1‑100 (MQ=255)
                    aTCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcgc                                                                              >  1:2337694/1‑100 (MQ=255)
                     tagTTTGGATCACCAGGCCAGTTTTCCAATTCCGGGGTAGTGATAAGTTGGCCGCCTTTTTCCATGCTGATAATCAGTACAGGTTGCAGGTTAgcgcgcg                                                                             <  2:1035855/98‑1 (MQ=255)
                     tCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcgcg                                                                             <  1:1864436/100‑1 (MQ=255)
                     tCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcgcg                                                                             <  2:1790658/100‑1 (MQ=255)
                            gttCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                      <  1:736230/98‑1 (MQ=255)
                            gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCTCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCATa                                                                      <  2:1351350/100‑1 (MQ=255)
                            gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                      >  2:671648/1‑100 (MQ=255)
                            gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                      >  2:2954/1‑100 (MQ=255)
                            gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGGCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                      <  1:1686630/100‑1 (MQ=255)
                              tCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                      <  1:525119/98‑1 (MQ=255)
                               cGCCAGGCTAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCACGCCGCGTAGAc                                                                   <  2:219586/100‑1 (MQ=255)
                               cGCCAGGCCAGTTTTCCACTTCCGTTGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                   <  1:697464/100‑1 (MQ=255)
                               cGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                   <  2:1540403/100‑1 (MQ=255)
                               cGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                   <  2:1621851/100‑1 (MQ=255)
                                gccaggccagTTTTCCACTTCCGTGGTGGTGTTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGTTTGCAGTTTGGCGCGCGCCGCGTAGAc                                                                   <  2:533008/99‑1 (MQ=255)
                                gccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAca                                                                  >  2:2218771/1‑100 (MQ=255)
                                gccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAca                                                                  >  1:2112812/1‑100 (MQ=255)
                                gccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAca                                                                  >  2:705476/1‑100 (MQ=255)
                                gccaggccagTTTTACACTTCCGTGGTGGTGGTCAGTTGGTCGCCTTTTTCCATGACGGTAATCAGTACAGGTTTTAGGTTGGCGCGCGCCGCGTAGAca                                                                  <  2:1667554/100‑1 (MQ=255)
                                    ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcag                                                              >  1:1554818/1‑100 (MQ=255)
                                     gccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                             <  1:1360538/100‑1 (MQ=255)
                                         gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                         >  1:1188149/1‑100 (MQ=255)
                                          ttttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCgg                                                           <  1:626784/97‑1 (MQ=255)
                                            ttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                      >  1:1732359/1‑100 (MQ=255)
                                            ttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                      >  1:1635392/1‑100 (MQ=255)
                                            ttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                      <  1:1343837/100‑1 (MQ=255)
                                                 cTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgc                                                 >  1:994906/1‑100 (MQ=255)
                                                       tggtggtggTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGcc                                           <  1:1921733/100‑1 (MQ=255)
                                                        ggtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCt                                          <  2:1732338/100‑1 (MQ=255)
                                                        ggtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCt                                          <  2:1635372/100‑1 (MQ=255)
                                                         gtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTg                                         <  1:492913/100‑1 (MQ=255)
                                                          tggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGa                                        <  1:1611821/100‑1 (MQ=255)
                                                           ggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGaa                                       <  1:1660108/100‑1 (MQ=255)
                                                            gtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCTGGGCCTGaa                                       >  2:1576397/1‑99 (MQ=255)
                                                             tggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAcc                                     <  1:554492/100‑1 (MQ=255)
                                                             tggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAcc                                     <  1:309851/100‑1 (MQ=255)
                                                             tggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAcc                                     <  1:679372/100‑1 (MQ=255)
                                                              ggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAccc                                    >  2:919261/1‑100 (MQ=255)
                                                              ggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGATTGCAGATTGGCGAGCGCTGCGTGAACCGCAGCGGTGTATCCCGCCGGGACTg                                         >  2:410271/1‑95 (MQ=255)
                                                                tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                  >  1:1715512/1‑100 (MQ=255)
                                                                tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                  <  2:996743/100‑1 (MQ=255)
                                                                 cAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                                 >  2:494072/1‑100 (MQ=255)
                                                                 cAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                                 <  2:489696/100‑1 (MQ=255)
                                                                      ggCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                            >  1:1392424/1‑100 (MQ=255)
                                                                       gccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGc                           >  1:734668/1‑100 (MQ=255)
                                                                       gccgccTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGc                           >  2:2386099/1‑100 (MQ=255)
                                                                        ccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCa                          >  1:1494678/1‑100 (MQ=255)
                                                                         cgccTTTTTCTTTGTCTTTAATCAGCACAGGATGCAGGTTGGCTCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAg                         >  2:1722921/1‑100 (MQ=255)
                                                                            cTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGttt                      <  2:1788609/100‑1 (MQ=255)
                                                                              ttttCCATGTCGGTAATCATCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGTCTGAACCCAGGATAAGCAGTTTAc                    <  2:671043/100‑1 (MQ=255)
                                                                                ttCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtg                  <  2:889975/100‑1 (MQ=255)
                                                                                 tCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgt                 <  2:1625215/100‑1 (MQ=255)
                                                                                 tCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgt                 <  1:1706728/100‑1 (MQ=255)
                                                                                  ccATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgtg                <  1:51287/100‑1 (MQ=255)
                                                                                   cATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt               >  1:189656/1‑100 (MQ=255)
                                                                                   cATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt               <  1:2386127/100‑1 (MQ=255)
                                                                                    aTGCCGGTAATCAGCACAGGTTGCAGGTGGGCGCGCGCCGCGTAGACAGGAGCGGGGTATTCCGCCGGGCCTGAACCCAGGTTAAGCAGTTTACTGTGtt              <  2:466978/100‑1 (MQ=255)
                                                                                     tGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGttt             <  1:563889/100‑1 (MQ=255)
                                                                                     tGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGttt             >  1:1770060/1‑100 (MQ=255)
                                                                                       ccGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTgg           >  2:1288963/1‑100 (MQ=255)
                                                                                       ccGGTAAGCAGCACAGGTTGCAGGTTGGCGCGCGCCTGGTAGACAGGAGAGGTTTATCCCGCCGGGCCTGTACCCAGGATAAGCAGTTTACTGTGttt             <  2:227258/98‑1 (MQ=255)
                                                                                        cGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGt          <  1:1345039/100‑1 (MQ=255)
                                                                                         ggTAATCAGTACAGGTTGCAGGTTGGGGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTc         <  2:2085734/100‑1 (MQ=255)
                                                                                          gTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCg        <  1:727421/100‑1 (MQ=255)
                                                                                            aaTCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTg      <  1:49251/100‑1 (MQ=255)
                                                                                              tCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc    >  1:654706/1‑100 (MQ=255)
                                                                                              tCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc    >  2:7908/1‑100 (MQ=255)
                                                                                              tCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc    <  2:1715492/100‑1 (MQ=255)
                                                                                              tCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc    >  1:714964/1‑100 (MQ=255)
                                                                                                aGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCa  >  1:2180467/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
ATTAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCA  >  NC_000913/931854‑932049

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: