Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,952,569 C→T 74.7% 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,5690CT74.7% 271.8 / 101.6 166A29V (GCG→GTG) ilvEbranched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (27/15);  new base T (61/62);  total (89/77)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 170 total aligned reads displayed.

GGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952473‑3952668
                                                                                                |                                                                                                   
 gAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTg                                                                                                 >  2:3733023/1‑100 (MQ=255)
  aaTATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGtgctgc                                                                                                <  2:225784/100‑1 (MQ=255)
    tataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCAc                                                                                              <  2:3986922/100‑1 (MQ=255)
    tataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCAc                                                                                              >  1:3445730/1‑100 (MQ=255)
     ataAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACt                                                                                             >  2:3336502/1‑100 (MQ=255)
      taAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTa                                                                                            >  2:1053819/1‑100 (MQ=255)
          aaTGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGc                                                                                        <  2:1595091/100‑1 (MQ=255)
          aaTGACCAAGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGCAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGc                                                                                        <  1:873070/100‑1 (MQ=255)
            tGACCACGAAGAAAGCTGATTACAGTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCAc                                                                                      <  1:203694/100‑1 (MQ=255)
               ccACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTc                                                                                   >  2:1165997/1‑100 (MQ=255)
                cACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCg                                                                                  >  1:2057424/1‑100 (MQ=255)
                cACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCg                                                                                  >  1:1756181/1‑100 (MQ=255)
                 aCGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCgg                                                                                 <  1:3309985/100‑1 (MQ=255)
                 aCGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCgg                                                                                 >  1:2519677/1‑100 (MQ=255)
                   gaagaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGtt                                                                               >  2:510055/1‑100 (MQ=255)
                     agaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGtttt                                                                             <  1:750269/100‑1 (MQ=255)
                      gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttttt                                                                            <  2:3546446/100‑1 (MQ=255)
                      gaaAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGttttt                                                                            <  2:1790562/100‑1 (MQ=255)
                       aaaGCTGATTACATTTGGTTCAATGGGGAGATGGCTCGCTGGGAAGACGCGAAGTTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTg                                                                           <  1:4352187/100‑1 (MQ=255)
                        aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGa                                                                          <  1:1353508/100‑1 (MQ=255)
                        aaGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGa                                                                          >  2:4048357/1‑100 (MQ=255)
                         aGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGaa                                                                         <  2:1795230/100‑1 (MQ=255)
                          gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                        <  2:2857524/100‑1 (MQ=255)
                          gCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAg                                                                        <  1:3471135/100‑1 (MQ=255)
                           cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:3224792/100‑1 (MQ=255)
                           cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       >  2:462676/1‑100 (MQ=255)
                           cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  2:3391382/100‑1 (MQ=255)
                           cTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAgg                                                                       <  1:2525682/100‑1 (MQ=255)
                              aTTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCAt                                                                    >  1:2599452/1‑100 (MQ=255)
                              aTTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCAt                                                                    >  2:788951/1‑100 (MQ=255)
                               ttACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATc                                                                   <  2:3724973/100‑1 (MQ=255)
                                 aCATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCg                                                                 >  1:2942472/1‑100 (MQ=255)
                                     ttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTGGCAAGTGCTGCACTGTGGCAGTTCGGTTTTTGAAGGAATCCGGTg                                                              >  1:4033979/1‑99 (MQ=255)
                                     ttGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTgn                                                             <  2:1582646/100‑2 (MQ=255)
                                       ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTa                                                           >  1:1654088/1‑100 (MQ=255)
                                       ggTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTa                                                           >  1:846527/1‑100 (MQ=255)
                                        gTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTAc                                                          <  1:3217777/100‑1 (MQ=255)
                                         ttCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACg                                                         >  1:1203320/1‑100 (MQ=255)
                                         ttCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACg                                                         <  1:1962577/100‑1 (MQ=255)
                                          tCAATGGAGAGATGGTTCGCTTGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCTGTTTTTGAAGGCATCCGTTGCTACg                                                         <  1:890282/99‑1 (MQ=255)
                                             aTGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTc                                                     >  1:3348554/1‑100 (MQ=255)
                                              tGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCg                                                    >  2:3384444/1‑100 (MQ=255)
                                                gggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATTGCACTTCGGGTTTTGAAGGCATCCGTTGCTACGACTCGca                                                  <  1:1387706/100‑1 (MQ=255)
                                                gggAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTCATGTCGAACATGCTGCACTATGGCACGTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGca                                                  <  1:2973993/100‑1 (MQ=255)
                                                  gagaTGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGcaca                                                <  2:627461/100‑1 (MQ=255)
                                                      tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            <  2:2120783/100‑1 (MQ=255)
                                                      tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            >  2:1054672/1‑100 (MQ=255)
                                                      tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            >  2:890611/1‑100 (MQ=255)
                                                      tGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAgg                                            >  2:2311564/1‑100 (MQ=255)
                                                       ggTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGa                                           >  2:3864053/1‑100 (MQ=255)
                                                       ggTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGAATCCGTTGCTACGACTCGCACAAAGGa                                           >  1:2374848/1‑100 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       >  1:4178883/1‑100 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       <  2:60078/100‑1 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       >  1:1708735/1‑100 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       <  2:3758491/100‑1 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       >  1:1818997/1‑100 (MQ=255)
                                                           cGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCgg                                       >  2:2522074/1‑100 (MQ=255)
                                                            gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                      >  2:1483635/1‑100 (MQ=255)
                                                            gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                      >  1:4292010/1‑100 (MQ=255)
                                                            gCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGt                                      >  1:3569644/1‑100 (MQ=255)
                                                             cTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTGTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGtt                                     >  1:480712/1‑100 (MQ=255)
                                                              tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:2821056/100‑1 (MQ=255)
                                                              tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:2333431/100‑1 (MQ=255)
                                                              tGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTg                                    <  1:3923617/100‑1 (MQ=255)
                                                              tGGGAAGACGCCAAGGTGCATGTGATGTCGCACGTGCTGCACTATCACACCTCGGTTTTTGAAGGCAGCCGTTGCTACGATTCGCACAAAGGACCGGtgg                                    <  1:1910217/100‑3 (MQ=255)
                                                               gggAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGt                                   >  2:2369916/1‑100 (MQ=255)
                                                                ggAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                  >  2:3297708/1‑100 (MQ=255)
                                                                ggAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTa                                  <  1:3883714/100‑1 (MQ=255)
                                                                 gAAGACGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                                 >  2:1878732/1‑100 (MQ=255)
                                                                 gAAGACGCGAAGGCGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAt                                 >  1:730500/1‑100 (MQ=255)
                                                                     aCGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                             <  2:4178914/100‑1 (MQ=255)
                                                                     aCGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                             <  1:3297689/100‑1 (MQ=255)
                                                                     aCGCGAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCg                             <  1:3805979/100‑1 (MQ=255)
                                                                      cgcgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGc                            >  2:2409134/1‑100 (MQ=255)
                                                                        cgAAGGTGCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                          >  1:3844473/1‑100 (MQ=255)
                                                                        cgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACGAAGGACCGGTTGTATTCCGcc                           >  1:3273828/1‑99 (MQ=255)
                                                                        cgAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCa                          >  2:1409381/1‑100 (MQ=255)
                                                                            ngTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCg                       >  1:1773900/2‑99 (MQ=255)
                                                                               gCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   <  2:2168389/100‑1 (MQ=255)
                                                                               gCATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAg                   >  2:1427171/1‑100 (MQ=255)
                                                                                cATGTGATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                  <  1:4048327/100‑1 (MQ=255)
                                                                                cATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGc                  >  2:1712931/1‑100 (MQ=255)
                                                                                  tgtgATGTCACACGCGCTGTAGTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACCCGCACAAAGGACCGGTTGTATTCCGCCGTCGTGAg                   <  1:1527194/97‑1 (MQ=255)
                                                                                     gATGTCGCACGTTCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGAATATg             <  1:895697/100‑1 (MQ=255)
                                                                                     gATGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg             <  1:134780/100‑1 (MQ=255)
                                                                                     gATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGTCCGGTTGTATTCCGCCATCGTGAGCATATg             <  1:610255/100‑1 (MQ=255)
                                                                                       tatCGCACGTGCTCCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCa           >  2:2996724/3‑100 (MQ=255)
                                                                                       tGTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGCCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATg             <  1:547931/98‑1 (MQ=255)
                                                                                        gTCGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAg          >  2:523084/1‑100 (MQ=255)
                                                                                          cGCACGTGCTGCACTATGGCCATTCGGTTTTTGACGGCATCCGTTGCTACGACTCGCCCAAAGGACCGGGTGTATTCCGCCATCGTGAGCATATGc            >  1:4151230/1‑96 (MQ=255)
                                                                                          cGCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGCTGTATTCCGCCATCGTCAGCATATGCa           >  1:1248492/1‑97 (MQ=255)
                                                                                          cGCACGTGCTGCACTATGGCACGTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        <  1:2678820/100‑1 (MQ=255)
                                                                                          cGCACGTGCTGCACTATCGAACTTCGGTTTTTGAATTCATCCGTTGCTACGCCTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        <  1:1045062/100‑1 (MQ=255)
                                                                                          cGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCg        >  2:1116516/1‑100 (MQ=255)
                                                                                           gCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCGTATGCAGCGt       >  1:2886203/1‑100 (MQ=255)
                                                                                           gCACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGt       <  1:1414883/100‑1 (MQ=255)
                                                                                            cACGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc      <  2:2886220/100‑1 (MQ=255)
                                                                                            cACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTc      >  1:4373898/1‑100 (MQ=255)
                                                                                             aCGTGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTAGTCCGCCATCGTGAGCATATGCAGCGTCt     <  1:1030450/100‑1 (MQ=255)
                                                                                                cgcTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCa  <  2:3278758/100‑1 (MQ=255)
                                                                                                |                                                                                                   
GGAATATAAAAATGACCACGAAGAAAGCTGATTACATTTGGTTCAATGGGGAGATGGTTCGCTGGGAAGACGCGAAGGTGCATGTGATGTCGCACGCGCTGCACTATGGCACTTCGGTTTTTGAAGGCATCCGTTGCTACGACTCGCACAAAGGACCGGTTGTATTCCGCCATCGTGAGCATATGCAGCGTCTGCA  >  NC_000913/3952473‑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: