Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,557 2 bp→GC coding (451‑452/930 nt) 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,5570CG100.0% 1022.2 / NA 315R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (168/147);  total (168/147)
*NC_0009133,950,5580GC99.7% 1002.1 / NA 316R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (167/148);  minor base A (1/0);  total (168/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT  >  NC_000913/3950470‑3950644
                                                                                       ||                                                                                      
ttCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                    <  1:8008283/92‑1 (MQ=255)
aTCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                    >  2:6200711/1‑93 (MQ=255)
  cGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCaaa                                                                                  <  2:8873886/93‑1 (MQ=255)
   gctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                 <  1:991977/93‑1 (MQ=255)
   gctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                 <  2:3337648/93‑1 (MQ=255)
    ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                                                >  2:8552982/1‑93 (MQ=255)
      gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                              <  1:1238366/93‑1 (MQ=255)
          tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                          >  1:144617/1‑93 (MQ=255)
           cccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATccc                                                                         >  1:8854107/1‑93 (MQ=255)
           cccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATccc                                                                         <  1:697406/93‑1 (MQ=255)
           cccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATccc                                                                         <  1:5177984/93‑1 (MQ=255)
            ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATccca                                                                        >  1:2119712/1‑92 (MQ=255)
             cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                       <  2:2165937/93‑1 (MQ=255)
             cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                       <  2:9318372/93‑1 (MQ=255)
             cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                       >  2:4516701/1‑93 (MQ=255)
              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                      >  1:7124792/1‑93 (MQ=255)
              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                      <  1:5426195/93‑1 (MQ=255)
              gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                      >  2:5413398/1‑93 (MQ=255)
               tGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcg                                                                     <  1:4396111/93‑1 (MQ=255)
               tGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATCGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcg                                                                     >  2:1281803/1‑93 (MQ=255)
                ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcgg                                                                    >  2:1543748/1‑93 (MQ=255)
                ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcgg                                                                    >  1:93550/1‑93 (MQ=255)
                 gggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggc                                                                   <  1:5824259/93‑1 (MQ=255)
                  ggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                  >  1:9765484/1‑93 (MQ=255)
                       gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                             >  1:3010727/1‑93 (MQ=255)
                       gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                             <  2:4262884/93‑1 (MQ=255)
                        tATCTGGGCGCAGAGGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                            >  1:8420918/1‑93 (MQ=255)
                        tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                            >  1:6655472/1‑93 (MQ=255)
                         aTCTGGGCGCAGAATCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                           >  1:6405161/1‑93 (MQ=255)
                         aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                           <  2:6697498/93‑1 (MQ=255)
                          tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                          >  1:2361485/1‑93 (MQ=255)
                          tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                          >  2:6847527/1‑93 (MQ=255)
                          tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGGTTACTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                          >  2:4765655/1‑93 (MQ=255)
                           cTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                                         >  1:8756186/1‑93 (MQ=255)
                            tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                        <  1:8460823/93‑1 (MQ=255)
                            tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                        <  1:8182736/93‑1 (MQ=255)
                            tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                        <  2:8413629/93‑1 (MQ=255)
                            tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                        <  1:8704568/93‑1 (MQ=255)
                             gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                       >  2:7645683/1‑93 (MQ=255)
                             gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                       <  2:5792159/93‑1 (MQ=255)
                               gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                     <  2:417805/93‑1 (MQ=255)
                                  cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTa                                                  <  2:7833587/93‑1 (MQ=255)
                                  cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTa                                                  >  2:1880673/1‑93 (MQ=255)
                                    gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                                <  2:2975912/93‑1 (MQ=255)
                                    gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                                >  2:8005959/1‑93 (MQ=255)
                                    gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                                >  1:6180091/1‑93 (MQ=255)
                                    gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGGAAc                                                >  2:7085044/1‑93 (MQ=255)
                                     aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACt                                               <  2:649634/93‑1 (MQ=255)
                                       gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAc                                             <  1:1850922/93‑1 (MQ=255)
                                         gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                           <  1:9021614/93‑1 (MQ=255)
                                          cTGGAGCAGGGGATGGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACCCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                          <  2:10037788/93‑1 (MQ=255)
                                          cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                          <  1:4972386/93‑1 (MQ=255)
                                           tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                         >  1:9467804/1‑93 (MQ=255)
                                            ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                        <  2:1397319/93‑1 (MQ=255)
                                               gCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                                     <  1:559551/93‑1 (MQ=255)
                                                gcggggAGCGATGCGATGGTTTACTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                                    <  1:1116981/91‑1 (MQ=255)
                                                cAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                                    <  2:8202588/93‑1 (MQ=255)
                                                cAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                                    <  1:3824480/93‑1 (MQ=255)
                                                 aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                   >  1:9606355/1‑93 (MQ=255)
                                                  ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                  >  2:6026553/1‑93 (MQ=255)
                                                  ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                  >  2:6572184/1‑93 (MQ=255)
                                                     gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                               >  1:1688546/1‑93 (MQ=255)
                                                     gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                               <  1:4683734/93‑1 (MQ=255)
                                                      aTCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgctg                              >  1:3365254/1‑93 (MQ=255)
                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                           <  1:8548313/93‑1 (MQ=255)
                                                         gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                           >  2:5997197/1‑93 (MQ=255)
                                                          atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                          >  1:1061564/1‑93 (MQ=255)
                                                          atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                          <  1:6060704/93‑1 (MQ=255)
                                                           tgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTggg                         <  1:2708828/93‑1 (MQ=255)
                                                            gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                        <  2:8784037/93‑1 (MQ=255)
                                                            gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                        >  2:5036046/1‑93 (MQ=255)
                                                             cgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGtg                       >  1:4908142/1‑92 (MQ=255)
                                                             cgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa                       <  1:4247433/93‑1 (MQ=255)
                                                             cgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa                       >  2:325803/1‑93 (MQ=255)
                                                              gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                      <  2:4911013/93‑1 (MQ=255)
                                                               atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGc                     <  1:5371538/93‑1 (MQ=255)
                                                                tgGTTTCCTCCTGGAACCGCGCAGCACCAAGCACCAACCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg                    <  2:3214342/93‑1 (MQ=255)
                                                                tgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg                    >  2:1070497/1‑93 (MQ=255)
                                                                tgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg                    <  1:7282147/93‑1 (MQ=255)
                                                                 ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                   <  2:3436828/93‑1 (MQ=255)
                                                                 ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                   >  1:7872011/1‑93 (MQ=255)
                                                                 ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                   <  2:7397746/93‑1 (MQ=255)
                                                                  gTTTCCTCCTGGAACCTCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                  >  2:3973684/1‑93 (MQ=255)
                                                                  gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGag                  >  2:2717155/1‑92 (MQ=255)
                                                                  gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                  >  2:3591077/1‑93 (MQ=255)
                                                                  gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGCAGCGaa                  >  1:2301846/1‑93 (MQ=255)
                                                                     tcctccTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcg               >  2:9464893/1‑93 (MQ=255)
                                                                       ctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcg             >  1:8062540/1‑93 (MQ=255)
                                                                        tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc            >  1:3269406/1‑93 (MQ=255)
                                                                         cctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc           >  2:9752304/1‑93 (MQ=255)
                                                                         cctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc           >  2:1995192/1‑93 (MQ=255)
                                                                          ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg          >  1:352738/1‑93 (MQ=255)
                                                                          ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg          >  2:3904802/1‑93 (MQ=255)
                                                                          ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg          <  2:2351303/93‑1 (MQ=255)
                                                                          ctGGAACCGCCCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg          <  1:5123536/93‑1 (MQ=255)
                                                                           tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc         >  2:2551191/1‑93 (MQ=255)
                                                                           tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgacgc         >  2:1270400/1‑93 (MQ=255)
                                                                             gacccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTCCCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa       <  1:2488401/90‑1 (MQ=255)
                                                                              aaCCGCCCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc      <  1:6618820/93‑1 (MQ=255)
                                                                                 cgcgCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGt   >  1:7550731/1‑93 (MQ=255)
                                                                                       ||                                                                                      
ATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTT  >  NC_000913/3950470‑3950644

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