Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,952,569 C→T 36.6% A29V (GCG→GTG)  ilvE → branched‑chain amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,952,5690CT36.6% 88.8 / 84.5 90A29V (GCG→GTG) ilvEbranched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (21/36);  new base T (15/18);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AAAGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952470‑3952668
                                                                                                   |                                                                                                   
aaaGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCAcgc                                                                                                     >  1:1741222/1‑100 (MQ=255)
aaaGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCAcgc                                                                                                     >  1:1378343/1‑100 (MQ=255)
  aGGAATCTAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGCAGACGCGAAGGTGCATGTGATGTCGCACgcgc                                                                                                   >  1:1352526/1‑100 (MQ=255)
   ggAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGCTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACgcgc                                                                                                   >  2:59560/1‑99 (MQ=255)
    gAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTg                                                                                                 <  1:8965/100‑1 (MQ=255)
    gAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTg                                                                                                 >  2:949953/1‑100 (MQ=255)
     aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTg                                                                                                 <  2:1184923/99‑1 (MQ=255)
     aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGc                                                                                                <  1:1055825/100‑1 (MQ=255)
     aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTACAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGtgctgc                                                                                                >  2:2184261/1‑100 (MQ=255)
      atatAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGc                                                                                                <  1:196854/99‑1 (MQ=255)
      atatAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCa                                                                                               <  2:2245180/100‑1 (MQ=255)
      atatAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCa                                                                                               <  2:2199431/100‑1 (MQ=255)
        ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                             <  2:1136510/100‑1 (MQ=255)
        ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACt                                                                                             <  2:1557027/100‑1 (MQ=255)
        ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTATAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                             >  2:2021504/1‑100 (MQ=255)
         taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTa                                                                                            >  1:619303/1‑100 (MQ=255)
          aaaaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTAt                                                                                           >  2:2124611/1‑100 (MQ=255)
           aaaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATg                                                                                          >  1:928630/1‑100 (MQ=255)
           aaaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATg                                                                                          <  1:1373993/100‑1 (MQ=255)
            aaaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATgg                                                                                         <  2:1747744/100‑1 (MQ=255)
             aaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGc                                                                                        >  1:1371961/1‑100 (MQ=255)
               tGACCACGAAGAAAGCTGATTACATTTGGTTGAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATgg                                                                                         <  1:974649/97‑1 (MQ=255)
                  ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTc                                                                                   >  1:989930/1‑100 (MQ=255)
                   cACTAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCg                                                                                  <  1:170999/100‑1 (MQ=255)
                     cGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCgg                                                                                 <  2:1793906/99‑1 (MQ=255)
                      gaagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGtt                                                                               <  2:668775/100‑1 (MQ=255)
                       aagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGCTGTCGCACGCGCTGCACTATGGCACTTCGGttt                                                                              >  2:875585/1‑100 (MQ=255)
                       aagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttt                                                                              <  1:1145584/100‑1 (MQ=255)
                       aagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttt                                                                              <  1:401392/100‑1 (MQ=255)
                       aagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGttt                                                                              <  2:734561/100‑1 (MQ=255)
                       aagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGttt                                                                              <  1:2008826/100‑1 (MQ=255)
                         gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGttttt                                                                            <  1:2141250/100‑1 (MQ=255)
                          aaaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTg                                                                           <  2:1983905/100‑1 (MQ=255)
                          aaaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTg                                                                           <  1:825584/100‑1 (MQ=255)
                           aaGCTGATTACATTTGGTTCAATGGGGATATGGTTCGCTGGGAAGACTCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGa                                                                          <  2:550503/100‑1 (MQ=255)
                           aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGa                                                                          <  1:1789280/100‑1 (MQ=255)
                             gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                        <  1:895166/100‑1 (MQ=255)
                              cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  1:1942516/100‑1 (MQ=255)
                              cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:2233049/100‑1 (MQ=255)
                              cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       >  1:1434435/1‑100 (MQ=255)
                              cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  1:747230/100‑1 (MQ=255)
                                   tACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTGGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATc                                                                   >  2:2059543/1‑99 (MQ=255)
                                      aTTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGtt                                                               >  1:1044401/1‑100 (MQ=255)
                                      aTTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGtt                                                               <  2:1247433/100‑1 (MQ=255)
                                        tttgTTTAATGGTGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGc                                                             <  1:580177/97‑1 (MQ=255)
                                        ttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGc                                                             <  2:1788163/100‑1 (MQ=255)
                                           gTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTAc                                                          >  1:2003663/1‑100 (MQ=255)
                                            ttCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACg                                                         <  2:542642/100‑1 (MQ=255)
                                             tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGa                                                        <  1:1192218/100‑1 (MQ=255)
                                             tCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGa                                                        >  2:331823/1‑100 (MQ=255)
                                               aaTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACt                                                      >  1:363356/1‑100 (MQ=255)
                                                aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCAGTTTTTGAAGGCATCCGTTGCTACGACTc                                                     >  1:1573376/1‑100 (MQ=255)
                                                  ggggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGc                                                   <  2:1184125/100‑1 (MQ=255)
                                                     gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcaca                                                <  1:1777394/100‑1 (MQ=255)
                                                       gaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACTCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcaca                                                <  2:1551949/98‑1 (MQ=255)
                                                         tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            <  1:134315/100‑1 (MQ=255)
                                                         tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            <  2:1680607/100‑1 (MQ=255)
                                                         tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            >  1:187205/1‑100 (MQ=255)
                                                          ggTGCGGTGGGAAGACGCGAAGGTGGATGTGATGTCGCACGTGCTGAACTATGGCACTTCGGGTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            <  2:811794/99‑1 (MQ=255)
                                                            ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                         <  1:952709/100‑1 (MQ=255)
                                                            ttCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGAcc                                         >  2:1486734/1‑100 (MQ=255)
                                                             tCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCg                                        <  1:434845/100‑1 (MQ=255)
                                                                 tGGGAAGACGCGAAGTTGCATGTGATGTCGCATTTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTAGGACTCGCACAAAGGACCGGTTg                                    <  1:1104901/100‑1 (MQ=255)
                                                                 tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    >  2:1719957/1‑100 (MQ=255)
                                                                 tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:849375/100‑1 (MQ=255)
                                                                 tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    >  2:1688288/1‑100 (MQ=255)
                                                                  gggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                                   >  1:1282700/1‑100 (MQ=255)
                                                                   ggAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGgtagta                                  >  1:641978/1‑100 (MQ=255)
                                                                    gAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                                 >  1:2026609/1‑100 (MQ=255)
                                                                    gAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                                 >  2:1639348/1‑100 (MQ=255)
                                                                       gACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTcc                              <  1:1122217/100‑1 (MQ=255)
                                                                        aCGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                             >  1:971467/1‑100 (MQ=255)
                                                                        aCGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                             <  2:371350/100‑1 (MQ=255)
                                                                         cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  1:318135/1‑100 (MQ=255)
                                                                         cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            <  2:1996951/100‑1 (MQ=255)
                                                                         cgcgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  2:852256/1‑100 (MQ=255)
                                                                           cgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTTATTGCATTCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                          >  2:1392606/1‑100 (MQ=255)
                                                                           cgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                          >  1:1311008/1‑100 (MQ=255)
                                                                               ggTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGt                      >  2:319310/1‑100 (MQ=255)
                                                                                gTGCATGTGATGTCGCATGCGCTGCACTATGGCACTTTGGATTTTGAAGGGATCCGTTGCTACGACTCGCACAAAGGACCGGTGGTATTGCGGCCTCGTg                     >  2:2005609/1‑100 (MQ=255)
                                                                                gTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTg                     <  1:1762317/100‑1 (MQ=255)
                                                                                gTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTg                     >  1:1250243/1‑100 (MQ=255)
                                                                                  gCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTTAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   <  2:1434425/100‑1 (MQ=255)
                                                                                  gCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   >  1:1420657/1‑100 (MQ=255)
                                                                                  gCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACCACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   <  2:1378337/100‑1 (MQ=255)
                                                                                   cATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                  <  1:160745/100‑1 (MQ=255)
                                                                                     tgtgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCat                <  2:971465/100‑1 (MQ=255)
                                                                                     tgtgATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCat                <  2:1498443/100‑1 (MQ=255)
                                                                                      gtgATGTCGCACGCGCTGCAATATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCata               >  2:1367725/1‑100 (MQ=255)
                                                                                        gATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg             <  1:1695388/100‑1 (MQ=255)
                                                                                             cGCACGCGCTGCACTATGGCGCTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        <  1:1110872/100‑1 (MQ=255)
                                                                                               cACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc      >  1:1809703/1‑100 (MQ=255)
                                                                                                aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc      <  2:1592775/99‑1 (MQ=255)
                                                                                                aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCt     >  2:1902991/1‑100 (MQ=255)
                                                                                                aCGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGCATTCCGCCATCGTGAGCATATGCAGCGTCt     >  2:945874/1‑100 (MQ=255)
                                                                                                  gcgcTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTg    <  2:641978/99‑1 (MQ=255)
                                                                                                   cgcTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa  <  2:298530/100‑1 (MQ=255)
                                                                                                   cgcTGCACGATGGCACTTCGTTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa  <  1:1755804/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAAGGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952470‑3952668

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: