Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,400 T→C R98R (CGT→CGC)  ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4000TC95.5% 461.2 / 7.7 156R98R (CGT→CGC) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (80/69);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 174 total aligned reads displayed.

GCTTCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGG  >  NC_000913/3950312‑3950489
                                                                                        |                                                                                         
gCTTCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCGCCAGCGCCTATATccgcccgc                                                                                       <  1:4489617/93‑1 (MQ=255)
 cTTCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCt                                                                                      <  1:1153912/93‑1 (MQ=255)
  ttCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTg                                                                                     >  2:4629162/1‑93 (MQ=255)
  ttCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTg                                                                                     <  2:3181045/93‑1 (MQ=255)
    cccGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGAt                                                                                   <  1:3747778/93‑1 (MQ=255)
    cccGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGAt                                                                                   <  1:3112737/93‑1 (MQ=255)
     ccGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATc                                                                                  <  1:3984932/93‑1 (MQ=255)
      cGGTTTCGCAGAGTATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCt                                                                                 <  2:2396535/93‑1 (MQ=255)
      cGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCt                                                                                 <  2:2396534/93‑1 (MQ=255)
       ggTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCtt                                                                                >  2:3434219/1‑93 (MQ=255)
       ggTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCtt                                                                                >  2:2233381/1‑93 (MQ=255)
       ggTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCtt                                                                                >  2:2573653/1‑93 (MQ=255)
        gTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtc                                                                               <  1:4128945/93‑1 (MQ=255)
         tttCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtcg                                                                              <  1:742969/93‑1 (MQ=255)
         tttCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtcg                                                                              >  2:2521012/1‑93 (MQ=255)
           tCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtcgtc                                                                            <  1:1707780/93‑1 (MQ=255)
           tCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtcgtc                                                                            <  2:1056781/93‑1 (MQ=255)
            cGCAGAGGATTGATGAGCTGATGGAAGCTTGTCTGGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTtcgtcg                                                                           <  1:1119602/93‑1 (MQ=255)
             gCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCgg                                                                          <  2:3915343/93‑1 (MQ=255)
              cAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGt                                                                         <  2:2790610/93‑1 (MQ=255)
              cAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGt                                                                         <  1:2927950/93‑1 (MQ=255)
                gagCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAACAACAATCTCACCAGCGCCTATCTGCGCCCGCGGATCTCCTTCGGTGa                                                                       >  1:1613641/1‑93 (MQ=255)
                gagCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGa                                                                       >  1:2439408/1‑93 (MQ=255)
                 agCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGAt                                                                      <  2:1766351/93‑1 (MQ=255)
                 agCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGAt                                                                      <  2:4733366/93‑1 (MQ=255)
                    aTTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGtt                                                                   >  1:384265/1‑93 (MQ=255)
                      tgatgaGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTgg                                                                 <  2:1297586/93‑1 (MQ=255)
                      tgatgaGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTgg                                                                 <  2:302393/93‑1 (MQ=255)
                      tgatgaGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTgg                                                                 >  1:4064007/1‑93 (MQ=255)
                        atgaGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCa                                                               >  1:3924099/1‑93 (MQ=255)
                          gaGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATg                                                             <  1:1222710/93‑1 (MQ=255)
                           aGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATgg                                                            <  1:4591399/93‑1 (MQ=255)
                             cTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGa                                                          <  2:3148420/93‑1 (MQ=255)
                              tGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAg                                                         >  2:1186111/1‑93 (MQ=255)
                              tGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAg                                                         >  1:327253/1‑93 (MQ=255)
                               gATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAga                                                        >  1:4617997/1‑92 (MQ=255)
                               gATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGt                                                        <  1:742125/93‑1 (MQ=255)
                               gATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTACCCAGCGCCTATATCCGCCCGCTGAGCTTCGTCGGTGATGTTGGCATgggcga                                                        >  1:1976400/1‑90 (MQ=255)
                                aTGGAAGCTTGTCGTTCGGTCATCCGCAAAACCAATCTCCCCAGCGCCTATATCCGCCCTCTGATCTTCGTCGGTGATGTTGGCATGGGAGTa                                                       <  2:1842968/93‑1 (MQ=255)
                                 tGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTaa                                                      <  2:757256/93‑1 (MQ=255)
                                  ggAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTaaa                                                     >  2:98070/1‑93 (MQ=255)
                                  ggAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTaaa                                                     <  1:4272104/93‑1 (MQ=255)
                                   gAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAAc                                                    <  2:3813352/93‑1 (MQ=255)
                                    aaGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAAcc                                                   >  1:367269/1‑93 (MQ=255)
                                    aaGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAAcc                                                   >  1:1361544/1‑93 (MQ=255)
                                    aaGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAAcc                                                   <  1:1364632/93‑1 (MQ=255)
                                      gCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCg                                                 >  1:2029868/1‑93 (MQ=255)
                                       cTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGc                                                >  2:4152252/1‑93 (MQ=255)
                                        ttGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGcc                                               <  1:2301316/93‑1 (MQ=255)
                                        ttGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGcc                                               >  2:168883/1‑93 (MQ=255)
                                        ttGTCGTGACGTGATCCCCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGATTAAACCCGcc                                               <  2:1290506/93‑1 (MQ=255)
                                            cgtgacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCg                                           <  2:2080379/93‑1 (MQ=255)
                                             gtgacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCgg                                          >  2:2882208/1‑93 (MQ=255)
                                              tgacgtgaTCCGCAAAAACAATCTCCCCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCggg                                         <  2:4347233/93‑1 (MQ=255)
                                              tgacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCggg                                         >  1:556559/1‑93 (MQ=255)
                                               gacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGa                                        <  1:3700285/93‑1 (MQ=255)
                                               gacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGa                                        >  2:4116384/1‑93 (MQ=255)
                                               gacgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGa                                        <  2:1592249/93‑1 (MQ=255)
                                                 cgtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCAGCCAGCGGGATa                                      >  2:1155893/1‑93 (MQ=255)
                                                  gtgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATAc                                     <  2:4077110/93‑1 (MQ=255)
                                                   tgaTCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACt                                    >  1:474571/1‑93 (MQ=255)
                                                      tCCGCAAAAACAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCaa                                 >  1:2891784/1‑93 (MQ=255)
                                                           aaaaaCAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGAc                            >  1:257165/1‑93 (MQ=255)
                                                           aaaaaCAACCCCCCCAGCGCCTTAACCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGAc                            <  1:3361195/93‑1 (MQ=255)
                                                            aaaaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACg                           >  2:4030144/1‑93 (MQ=255)
                                                             aaaCAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGt                          >  2:2714231/1‑93 (MQ=255)
                                                             aaaCAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGt                          >  2:644220/1‑93 (MQ=255)
                                                              aaCAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTg                         >  1:2354167/1‑93 (MQ=255)
                                                              aaCAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTg                         >  2:3237226/1‑93 (MQ=255)
                                                               aCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGa                        >  1:204271/1‑93 (MQ=255)
                                                                cAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGAt                       <  2:3296007/93‑1 (MQ=255)
                                                                cAATCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGAt                       >  1:125569/1‑93 (MQ=255)
                                                                  aTCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTa                     <  1:3000264/93‑1 (MQ=255)
                                                                  aTCTCACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTa                     >  1:1786414/1‑93 (MQ=255)
                                                                    ctcACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATc                   >  2:648714/1‑93 (MQ=255)
                                                                     tcACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCg                  >  1:973893/1‑93 (MQ=255)
                                                                     gcACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCg                  >  2:1932298/2‑93 (MQ=255)
                                                                      ccccAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgc                 <  1:2895439/91‑1 (MQ=255)
                                                                      cACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgc                 <  1:4365235/93‑1 (MQ=255)
                                                                      cACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgc                 >  2:2020382/1‑93 (MQ=255)
                                                                      cACCAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgc                 >  1:931039/1‑93 (MQ=255)
                                                                       cccAGCGTCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAACCCCGCCAGCGGGATACTCAACCGACGTGATTATCgct                <  1:2393748/92‑1 (MQ=255)
                                                                        ccAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgctg               >  2:2685368/1‑93 (MQ=255)
                                                                        ccAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgctg               <  1:2956712/93‑1 (MQ=255)
                                                                        ccAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgctg               >  2:3272866/1‑93 (MQ=255)
                                                                        ccAGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgctg               >  1:1341512/1‑93 (MQ=255)
                                                                          aGCGCCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCgctgct             <  1:4445432/93‑1 (MQ=255)
                                                                           gcgcCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCtt            <  1:2668581/93‑1 (MQ=255)
                                                                           gcgcCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATAGCTGCtt            >  2:2054097/1‑93 (MQ=255)
                                                                            cgcCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGGTGCttt           <  2:3380511/93‑1 (MQ=255)
                                                                             gcCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTc          >  2:324791/1‑93 (MQ=255)
                                                                             gcCTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTc          >  2:1395521/1‑93 (MQ=255)
                                                                               cTATATCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTccc        <  1:1458168/93‑1 (MQ=255)
                                                                                 atacCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGt      <  2:4342567/93‑1 (MQ=255)
                                                                                  tatCCGCCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAACGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTg     >  1:1278815/1‑93 (MQ=255)
                                                                                   agccgcccgcTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTgg    <  1:1693655/91‑1 (MQ=255)
                                                                                    tccgcccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTggg   >  1:1248338/1‑93 (MQ=255)
                                                                                    tccgcccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTggg   <  1:3895323/93‑1 (MQ=255)
                                                                                    tccgaccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTggg   <  1:3688215/93‑1 (MQ=255)
                                                                                     ccgcccgcTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTgggg  <  1:339723/93‑1 (MQ=255)
                                                                                        |                                                                                         
GCTTCCCGGTTTCGCAGAGCATTGATGAGCTGATGGAAGCTTGTCGTGACGTGATCCGCAAAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGG  >  NC_000913/3950312‑3950489

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