Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,177,509 (C)6→7 coding (402/630 nt) nudF ← ADP‑sugar pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,5031.C92.7% 1069.4 / 69.7 273coding (408/630 nt)nudFADP‑sugar pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base . (10/10);  new base C (162/91);  total (172/101)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

ATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCAC  >  NC_000913/3177404‑3177598
                                                                                                    |                                                                                                
 tGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cc                                                                                               >  2:5901579/1‑101 (MQ=255)
 tGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cc                                                                                               >  1:5901576/1‑101 (MQ=255)
   aaCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccc                                                                                             >  2:5901583/1‑101 (MQ=255)
   aaCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccc                                                                                             >  1:5901584/1‑101 (MQ=255)
      gcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑ccccccc                                                                                          <  1:5900984/101‑2 (MQ=255)
      gcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑ccccccc                                                                                          >  2:5901589/1‑100 (MQ=255)
      gcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑ccccccc                                                                                          <  1:5901365/101‑2 (MQ=255)
        gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCcgg                                                                                        >  1:5901592/1‑101 (MQ=255)
        gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCcgg                                                                                        <  1:5901145/101‑1 (MQ=255)
        gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCcgg                                                                                        >  2:5901595/1‑101 (MQ=255)
         aaTATCTTCGTTTTCATCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccc                                                                                             <  1:5901209/95‑1 (MQ=255)
          atatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccccccg                                                                                         >  2:5901602/1‑96 (MQ=255)
          atatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCcgggc                                                                                      >  1:5901600/1‑99 (MQ=255)
              cTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGc                                                                                   >  2:5901607/1‑101 (MQ=255)
               ttCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑c                                                                                                <  1:5901126/86‑1 (MQ=255)
                tCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCa                                                                                 >  1:5901609/1‑101 (MQ=255)
                 cGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAg                                                                                <  1:5901039/101‑1 (MQ=255)
                  gTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAgg                                                                               >  2:5901612/1‑101 (MQ=255)
                   ttttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGa                                                                              >  2:5901616/1‑101 (MQ=255)
                   ttttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGa                                                                              <  1:5901257/101‑1 (MQ=255)
                    tttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaa                                                                             >  1:5901617/1‑101 (MQ=255)
                     ttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAgg                                                                               >  2:5901625/1‑98 (MQ=255)
                     ttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaaa                                                                            >  1:5901622/1‑101 (MQ=255)
                       catcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACt                                                                          >  1:5901628/1‑101 (MQ=255)
                       catcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACt                                                                          >  1:5901631/1‑101 (MQ=255)
                         tcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTa                                                                        >  2:5901635/1‑101 (MQ=255)
                           aGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAac                                                                      >  1:5901640/1‑101 (MQ=255)
                           aGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAac                                                                      >  2:5901641/1‑101 (MQ=255)
                             ccAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAacac                                                                    >  1:5901646/1‑101 (MQ=255)
                               aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCg                                                                  >  1:5901650/1‑101 (MQ=255)
                               aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCg                                                                  >  2:5901649/1‑101 (MQ=255)
                                 aCCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                                <  1:5901094/101‑1 (MQ=255)
                                  ccGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGtt                                                               <  1:5901504/101‑1 (MQ=255)
                                   cGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaaa                                                                            >  2:5901656/1‑87 (MQ=255)
                                   cGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGttt                                                              >  1:5901655/1‑101 (MQ=255)
                                     tGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                            >  1:5901661/1‑101 (MQ=255)
                                      gAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                           >  2:5901662/1‑101 (MQ=255)
                                        aTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTcc                                                         >  2:5901665/1‑101 (MQ=255)
                                          aCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGt                                                       >  1:5901669/1‑101 (MQ=255)
                                            cGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                                     >  2:5901672/1‑101 (MQ=255)
                                             gCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                                    >  1:5901675/1‑101 (MQ=255)
                                              cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                                   >  1:5901678/1‑101 (MQ=255)
                                               ttGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTTACTGGTG‑‑CCCCCCcgg                                                                                        <  2:5901085/62‑1 (MQ=255)
                                                tGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGa                                                                              >  2:5901682/1‑72 (MQ=255)
                                                 gCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTa                                                >  1:5901683/1‑101 (MQ=255)
                                                   ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                               <  1:5901322/101‑1 (MQ=255)
                                                   ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTg                                                             >  1:5901689/1‑86 (MQ=255)
                                                   ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATc                                              >  2:5901688/1‑101 (MQ=255)
                                                    gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAg                                                                                <  1:5900780/66‑1 (MQ=255)
                                                    gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                             <  1:5901235/101‑1 (MQ=255)
                                                       gTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGa                                                                              <  1:5901163/65‑1 (MQ=255)
                                                        tGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                         >  2:5901694/1‑101 (MQ=255)
                                                        tGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                         >  1:5901695/1‑101 (MQ=255)
                                                            gTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGcc                                     >  2:5901700/1‑101 (MQ=255)
                                                              ccACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctc                                   >  2:5901703/1‑101 (MQ=255)
                                                               cACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                                  >  1:5901704/1‑101 (MQ=255)
                                                                aCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                                 >  1:5901708/1‑101 (MQ=255)
                                                                    cGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCtt                                                                                     <  2:5901277/45‑1 (MQ=255)
                                                                    cGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                        <  1:5901070/74‑1 (MQ=255)
                                                                     gCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATc                            >  1:5901713/1‑101 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAgg                                                                               <  2:5901218/49‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaaa                                                                            <  2:5901321/52‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTa                                                                        <  2:5901031/56‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCg                                                                  <  2:5901309/62‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTcc                                                         <  2:5901240/71‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                      <  2:5901226/74‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                      <  2:5901174/74‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                                    <  2:5901453/76‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                                   <  2:5901254/77‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                             <  1:5901173/83‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                        <  2:5901378/88‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                        <  1:5901195/88‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                                      <  2:5900916/90‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCct                                    <  2:5901410/92‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCa                               <  2:5901559/97‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                           >  1:5901717/1‑101 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                           >  2:5901716/1‑101 (MQ=255)
                                                                      cccACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCcgg                                                                                        <  2:5901072/39‑1 (MQ=255)
                                                                      cccACCATAATTGACGAACGCCCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                      <  2:5900966/74‑1 (MQ=255)
                                                                       ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTg                                                             <  2:5901119/66‑1 (MQ=255)
                                                                       ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                                  <  2:5901165/77‑1 (MQ=255)
                                                                        caccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGtt                                                               <  2:5901101/63‑1 (MQ=255)
                                                                         accaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCtt                        >  2:5901721/1‑101 (MQ=255)
                                                                            aTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCa                               >  1:5901725/1‑91 (MQ=255)
                                                                             tAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgcg                    >  2:5901726/1‑101 (MQ=255)
                                                                              aaTTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa                   >  1:5901729/1‑101 (MQ=255)
                                                                               aTTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGAc                  >  2:5901732/1‑101 (MQ=255)
                                                                                 tGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACgg                <  1:5901269/101‑1 (MQ=255)
                                                                                  gACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGc                          >  1:5901739/1‑90 (MQ=255)
                                                                                  gACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACggg               >  2:5901737/1‑101 (MQ=255)
                                                                                      aaCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc           >  2:5901746/1‑101 (MQ=255)
                                                                                      aaCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc           >  1:5901743/1‑101 (MQ=255)
                                                                                      aaCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc           >  1:5901748/1‑101 (MQ=255)
                                                                                         gCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATc        >  1:5901753/1‑101 (MQ=255)
                                                                                         gCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATc        >  2:5901751/1‑101 (MQ=255)
                                                                                           tCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCtt      <  1:5901216/101‑1 (MQ=255)
                                                                                           tCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCtt      >  1:5901756/1‑101 (MQ=255)
                                                                                            cACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTc     >  2:5901762/1‑101 (MQ=255)
                                                                                            cACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTc     >  1:5901761/1‑101 (MQ=255)
                                                                                               tGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCcac  >  1:5901767/1‑101 (MQ=255)
                                                                                                    |                                                                                                
ATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCAC  >  NC_000913/3177404‑3177598

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