Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 329,174 A→G Y163C (TAC→TGC)  betT → high‑affinity choline transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,1740AG99.3% 446.6 / ‑4.7 140Y163C (TAC→TGC) betThigh‑affinity choline transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (74/65);  total (74/66)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

Only 100 of 146 total aligned reads displayed.

GGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329085‑329263
                                                                                         |                                                                                         
ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        >  2:2836113/1‑93 (MQ=255)
ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        <  2:4038713/93‑1 (MQ=255)
ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        <  1:1448235/93‑1 (MQ=255)
ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        <  2:1187661/93‑1 (MQ=255)
 gCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCttt                                                                                       <  1:2449725/93‑1 (MQ=255)
   cgcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAg                                                                                     >  1:2901119/1‑93 (MQ=255)
   cgcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAg                                                                                     <  2:2336076/93‑1 (MQ=255)
    gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  1:1629751/1‑93 (MQ=255)
    gcgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGAATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    >  1:1629750/1‑93 (MQ=255)
    gagggTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGc                                                                                    <  2:1729169/89‑1 (MQ=255)
     cgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCt                                                                                   <  1:4324629/93‑1 (MQ=255)
      gcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCtt                                                                                  >  2:4140147/1‑92 (MQ=255)
      gcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTa                                                                                  <  1:1373422/93‑1 (MQ=255)
       cgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGGTTTAGCTAt                                                                                 >  2:4076756/1‑93 (MQ=255)
        gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATc                                                                                <  2:2567198/93‑1 (MQ=255)
          cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  1:4258409/93‑1 (MQ=255)
          cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  2:2512290/93‑1 (MQ=255)
          cAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGt                                                                              <  1:1207223/93‑1 (MQ=255)
           aGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGtt                                                                             >  1:3518665/1‑93 (MQ=255)
            ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            <  2:675586/93‑1 (MQ=255)
               gATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATaa                                                                         >  1:13239/1‑93 (MQ=255)
                aTGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAt                                                                        <  1:3283026/93‑1 (MQ=255)
                aTGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAAt                                                                        <  2:4524537/93‑1 (MQ=255)
                 tGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAtt                                                                       >  2:802397/1‑93 (MQ=255)
                   gTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTg                                                                     >  2:1440702/1‑93 (MQ=255)
                      tGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCg                                                                  >  2:1131400/1‑93 (MQ=255)
                       ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 >  1:193231/1‑93 (MQ=255)
                       ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  1:1816157/93‑1 (MQ=255)
                        gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTTTAATTTTCCGCt                                                                >  1:155175/1‑93 (MQ=255)
                        gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCt                                                                <  1:3006290/93‑1 (MQ=255)
                          cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTACCTATCGTTATAATTTGCCGCTCa                                                              >  2:2579048/1‑93 (MQ=255)
                           gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAc                                                             >  2:3827195/1‑93 (MQ=255)
                           gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCAc                                                             >  1:2971600/1‑93 (MQ=255)
                             tGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCa                                                           >  1:3546249/1‑93 (MQ=255)
                              gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCAt                                                          <  2:4581325/93‑1 (MQ=255)
                              gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGAGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCCCCAt                                                          >  2:3369881/1‑93 (MQ=255)
                               tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         >  1:2511973/1‑93 (MQ=255)
                                ttCCCTACGGCTTAACCGGCTGGTTGATGTATGCCCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATcc                                                        <  2:3683563/93‑1 (MQ=255)
                                ttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATcc                                                        >  1:3197992/1‑93 (MQ=255)
                                 tCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCTGCTCACCATCCg                                                       <  1:3781556/93‑1 (MQ=255)
                                   ccTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     <  2:1078125/92‑1 (MQ=255)
                                   aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     >  1:4731299/1‑93 (MQ=255)
                                   aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     <  2:3453725/93‑1 (MQ=255)
                                    cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  1:2010871/93‑1 (MQ=255)
                                    cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  1:3081317/93‑1 (MQ=255)
                                    cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  1:3317280/93‑1 (MQ=255)
                                     tACGGCTTAACGGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCg                                                   >  2:2749528/1‑93 (MQ=255)
                                       cGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgc                                                 >  1:2095161/1‑93 (MQ=255)
                                          cTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCt                                              >  2:248963/1‑93 (MQ=255)
                                           ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             >  1:778226/1‑93 (MQ=255)
                                           ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             >  2:2819388/1‑93 (MQ=255)
                                           ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             <  2:3730711/93‑1 (MQ=255)
                                             aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTa                                           >  1:2717883/1‑93 (MQ=255)
                                              aCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAc                                          >  1:3813882/1‑93 (MQ=255)
                                               ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  1:3814924/93‑1 (MQ=255)
                                                cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAccc                                        <  2:2956951/93‑1 (MQ=255)
                                                  gCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGa                                      <  2:2182417/93‑1 (MQ=255)
                                                    tGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATc                                    >  2:3552172/1‑93 (MQ=255)
                                                    tGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATc                                    >  1:426857/1‑93 (MQ=255)
                                                     ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   <  2:539593/93‑1 (MQ=255)
                                                       tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 >  2:4428729/1‑93 (MQ=255)
                                                       tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 <  1:1878516/93‑1 (MQ=255)
                                                       tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 <  1:1106549/93‑1 (MQ=255)
                                                       tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 <  2:2183367/93‑1 (MQ=255)
                                                         gATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCgg                               <  1:3771938/93‑1 (MQ=255)
                                                          aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGt                              <  1:4386281/93‑1 (MQ=255)
                                                          aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGt                              >  2:3920902/1‑93 (MQ=255)
                                                            gTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTaa                            >  2:4279434/1‑93 (MQ=255)
                                                              ttGCGCTGATGGGCATTGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  2:326184/92‑1 (MQ=255)
                                                              aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  1:537889/93‑1 (MQ=255)
                                                              aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  2:326185/93‑1 (MQ=255)
                                                                  gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAt                      >  1:3297648/1‑93 (MQ=255)
                                                                  gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAt                      <  2:3367942/93‑1 (MQ=255)
                                                                   cTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAtt                     >  1:3436368/1‑93 (MQ=255)
                                                                     gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   <  1:2068372/93‑1 (MQ=255)
                                                                      aTGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAAc                  <  2:3015471/93‑1 (MQ=255)
                                                                       tGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACg                 >  1:1813501/1‑93 (MQ=255)
                                                                        gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  1:4025445/93‑1 (MQ=255)
                                                                        gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  1:3834090/93‑1 (MQ=255)
                                                                        gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  2:1987454/93‑1 (MQ=255)
                                                                        gggAATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                <  2:1261312/93‑1 (MQ=255)
                                                                          gCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGc              <  2:1649097/93‑1 (MQ=255)
                                                                           cATGGCGCTCGGATGCTTCAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGcc             >  1:800638/1‑93 (MQ=255)
                                                                            aTGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCg            >  1:1035784/1‑93 (MQ=255)
                                                                             tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           <  1:628570/93‑1 (MQ=255)
                                                                             tGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGa           >  2:3759406/1‑93 (MQ=255)
                                                                              ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          <  2:2046732/93‑1 (MQ=255)
                                                                              ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          <  2:2046731/93‑1 (MQ=255)
                                                                              ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          <  2:3617790/93‑1 (MQ=255)
                                                                                 gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  2:4366213/1‑93 (MQ=255)
                                                                                 gcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAgg       >  1:3560700/1‑93 (MQ=255)
                                                                                   tgggATACTTTAGTTATCGTTATAATTTGCCTCTCACCATCCGCTCGGCTTTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     <  2:3055592/91‑1 (MQ=255)
                                                                                   tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     <  1:817747/93‑1 (MQ=255)
                                                                                   tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     >  1:4238976/1‑93 (MQ=255)
                                                                                   tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     >  2:2857328/1‑93 (MQ=255)
                                                                                    cGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa    >  2:1855400/1‑93 (MQ=255)
                                                                                     ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCAc   >  2:2760229/1‑93 (MQ=255)
                                                                                      gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  1:4646689/1‑93 (MQ=255)
                                                                                      gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  >  2:2921136/1‑93 (MQ=255)
                                                                                      gATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGGAAACGGATTAACGGGCCGATAGGTCACt  >  1:511956/1‑93 (MQ=255)
                                                                                         |                                                                                         
GGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329085‑329263

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