Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,540 A→G I350T (ATC→ACC)  mlc ← putative NAGC‑like transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,5400AG94.4% 242.2 / 4.7 89I350T (ATC→ACC) mlcputative NAGC‑like transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (39/45);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

GGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGC  >  NC_000913/1665448‑1665629
                                                                                            |                                                                                         
ggTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGa                                                                                           <  1:3282799/93‑1 (MQ=255)
  ttAGAAAACTGAGTTCTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGgta                                                                                         <  2:861302/93‑1 (MQ=255)
  ttAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGgta                                                                                         <  1:3939818/93‑1 (MQ=255)
   tAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGgtat                                                                                        >  1:2542047/1‑93 (MQ=255)
    aGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATc                                                                                       >  2:4486234/1‑93 (MQ=255)
    aGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATc                                                                                       <  1:100272/93‑1 (MQ=255)
     gAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCt                                                                                      >  2:1720806/1‑93 (MQ=255)
       aaaCTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGc                                                                                    >  2:1710508/1‑93 (MQ=255)
       aaaCTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGc                                                                                    <  1:3705068/93‑1 (MQ=255)
        aaCTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGcc                                                                                   >  1:513549/1‑93 (MQ=255)
          cTGAGTGCTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGc                                                                                 <  1:4455647/93‑1 (MQ=255)
          cTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGc                                                                                 >  1:925507/1‑93 (MQ=255)
          cTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGc                                                                                 <  1:2670031/93‑1 (MQ=255)
           tGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCt                                                                                <  2:850926/93‑1 (MQ=255)
           tGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCt                                                                                >  2:2873437/1‑93 (MQ=255)
              gTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTa                                                                             >  2:2330264/1‑93 (MQ=255)
               tACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTAc                                                                            <  2:2551751/93‑1 (MQ=255)
               tACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTAc                                                                            >  2:4087304/1‑93 (MQ=255)
                aCTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACt                                                                           <  2:3315483/93‑1 (MQ=255)
                aCTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACt                                                                           >  2:3460707/1‑93 (MQ=255)
                 ctctCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACtt                                                                          >  2:3758809/1‑93 (MQ=255)
                 ctctCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACtt                                                                          >  2:306198/1‑93 (MQ=255)
                  tctcAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTa                                                                         >  1:3543759/1‑93 (MQ=255)
                  tctcAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTa                                                                         <  1:2193733/93‑1 (MQ=255)
                    tcAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGTATATCTGCCGCTTTACTTAAc                                                                       <  2:2981759/93‑1 (MQ=255)
                    tcAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAAc                                                                       <  2:2014414/93‑1 (MQ=255)
                    tcAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGGATCTGCCGCTTTACTTAAc                                                                       >  2:2553457/1‑93 (MQ=255)
                     cAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACg                                                                      >  2:580021/1‑93 (MQ=255)
                       aCGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGt                                                                    >  1:3237428/1‑93 (MQ=255)
                       aCGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGt                                                                    >  1:4018975/1‑93 (MQ=255)
                        cGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTg                                                                   <  1:1150574/93‑1 (MQ=255)
                         gCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGa                                                                  >  1:3773574/1‑93 (MQ=255)
                          cTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGGGAg                                                                 >  2:1971996/1‑93 (MQ=255)
                           tGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGc                                                                >  1:2310045/1‑93 (MQ=255)
                            gATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGcc                                                               >  2:46465/1‑93 (MQ=255)
                            gATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGcc                                                               >  2:1522877/1‑93 (MQ=255)
                              tgtgCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGCGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCaa                                                             >  1:3724592/1‑93 (MQ=255)
                               gtgCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAAt                                                            >  2:895878/1‑93 (MQ=255)
                               gtgCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACAGGGAAGATGGTATCTGCCGCTTTACTTAACGGTGAGCCGAt                                                            >  2:3067329/1‑93 (MQ=255)
                                tgCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATc                                                           >  1:1077062/1‑93 (MQ=255)
                                tgCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATc                                                           <  1:4011080/93‑1 (MQ=255)
                                 gCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCa                                                          <  1:3941074/93‑1 (MQ=255)
                                   tGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGt                                                        >  2:4171903/1‑93 (MQ=255)
                                      cTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAtt                                                     <  2:3928150/93‑1 (MQ=255)
                                       tataCGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAttt                                                    >  1:72442/1‑93 (MQ=255)
                                        ataCACAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAtttt                                                   <  1:1760721/93‑1 (MQ=255)
                                         taCGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAttttt                                                  >  1:450590/1‑93 (MQ=255)
                                            gCAGGAAGGGCCTGCTGGCGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgt                                               <  2:1681910/93‑1 (MQ=255)
                                            gCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAtttttttt                                               >  2:1395443/1‑91 (MQ=255)
                                            gCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgt                                               <  1:891190/93‑1 (MQ=255)
                                            gCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgt                                               >  2:2796590/1‑93 (MQ=255)
                                            gCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgt                                               >  2:1840508/1‑93 (MQ=255)
                                             cAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgtg                                              <  1:114698/93‑1 (MQ=255)
                                             cAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgtg                                              <  2:4091461/93‑1 (MQ=255)
                                              aGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTgg                                             <  2:4174742/93‑1 (MQ=255)
                                               ggAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTggg                                            >  1:491258/1‑93 (MQ=255)
                                               ggAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTggg                                            <  2:2235702/93‑1 (MQ=255)
                                                gAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGt                                           <  2:3121652/93‑1 (MQ=255)
                                                    ggCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaa                                       <  1:2228646/93‑1 (MQ=255)
                                                     gCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaat                                      <  2:2454531/93‑1 (MQ=255)
                                                     gCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaat                                      <  2:849278/93‑1 (MQ=255)
                                                     gCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaat                                      <  2:3888204/93‑1 (MQ=255)
                                                        tgctgACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaataaa                                   <  2:1795947/93‑1 (MQ=255)
                                                          ctgACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAAtt                                 <  1:4154090/93‑1 (MQ=255)
                                                          ctgACGGATGCTGGCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAAtt                                 <  2:1900478/93‑1 (MQ=255)
                                                           tgACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTc                                <  1:949072/93‑1 (MQ=255)
                                                           tgACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTc                                <  1:232241/93‑1 (MQ=255)
                                                            gACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCa                               >  1:400367/1‑93 (MQ=255)
                                                            gACGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCa                               >  1:4419017/1‑93 (MQ=255)
                                                             aCGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCAc                              <  2:3427017/93‑1 (MQ=255)
                                                              cGGATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCAcc                             >  1:2810141/1‑93 (MQ=255)
                                                               ggATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCa                            <  1:924965/93‑1 (MQ=255)
                                                               ggATGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCa                            >  1:1305811/1‑93 (MQ=255)
                                                                  tGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCatga                         >  1:2400880/1‑93 (MQ=255)
                                                                  tGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCatga                         >  2:3487899/1‑93 (MQ=255)
                                                                  tGCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTAACCatga                         >  1:1210330/1‑93 (MQ=255)
                                                                   gCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCatgat                        <  2:2301186/93‑1 (MQ=255)
                                                                   gCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCatgat                        <  1:3680875/93‑1 (MQ=255)
                                                                   gCTGTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGCTAAATAAATTCACCatgat                        <  2:4578611/93‑1 (MQ=255)
                                                                      gTCTGAGATGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGc                     <  2:3723795/93‑1 (MQ=255)
                                                                           agaTGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGaa                <  2:1808298/93‑1 (MQ=255)
                                                                            gaTGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAAt               <  2:4288971/93‑1 (MQ=255)
                                                                            gaTGACAGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTTACTATGATGGCAAGAAt               >  2:3027888/1‑93 (MQ=255)
                                                                             aTGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATg              <  1:3370778/93‑1 (MQ=255)
                                                                             aTGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATg              <  1:2759232/93‑1 (MQ=255)
                                                                             aTGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATg              <  1:355187/93‑1 (MQ=255)
                                                                              tGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATgc             <  2:3999600/93‑1 (MQ=255)
                                                                              tGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATgc             <  2:2745089/93‑1 (MQ=255)
                                                                              tGACCGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATgc             <  2:115262/93‑1 (MQ=255)
                                                                                 ccGGGAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGcc          >  1:4220623/1‑93 (MQ=255)
                                                                                    ggAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGa       >  2:3839022/1‑93 (MQ=255)
                                                                                     gAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGAc      <  2:3297821/93‑1 (MQ=255)
                                                                                     gAAGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGAc      <  1:21491/93‑1 (MQ=255)
                                                                                      aaGAGGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACa     <  2:3156259/93‑1 (MQ=255)
                                                                                      aaGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACa     <  1:1949397/93‑1 (MQ=255)
                                                                                         agGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATgc  <  1:850543/93‑1 (MQ=255)
                                                                                         agGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATgc  >  2:2133753/1‑93 (MQ=255)
                                                                                         agGGTATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATgc  <  2:905849/93‑1 (MQ=255)
                                                                                            |                                                                                         
GGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGC  >  NC_000913/1665448‑1665629

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