Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 232,218 A→G D98G (GAC→GGC)  yafE → putative biotin synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,2180AG98.2% 521.2 / ‑0.9 169D98G (GAC→GGC) yafEputative biotin synthesis protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (92/73);  total (93/76)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 173 total aligned reads displayed.

GTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCA  >  NC_000913/232128‑232307
                                                                                          |                                                                                         
gTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCg                                                                                         >  1:4667049/1‑93 (MQ=255)
 ttATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGt                                                                                        >  1:2734222/1‑93 (MQ=255)
  tATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTa                                                                                       <  2:928367/93‑1 (MQ=255)
    ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAAt                                                                                     >  1:808609/1‑93 (MQ=255)
         cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtctc                                                                                <  2:130569/93‑1 (MQ=255)
         cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtctc                                                                                >  1:2937433/1‑93 (MQ=255)
         cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGtctc                                                                                >  1:2937432/1‑93 (MQ=255)
          ccaccaCGGGCATGATGTGGTTGCACCACGGCGAGAAGTGAATAGGTTATTGAAACCTGGCGGAGGGCTGATTTTGATGGCGTAATGTCTCCg                                                                              <  2:1967578/93‑1 (MQ=255)
          ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTcc                                                                               >  2:1166766/1‑93 (MQ=255)
           cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCg                                                                              >  1:179003/1‑93 (MQ=255)
           cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCg                                                                              <  1:1496663/93‑1 (MQ=255)
           cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCg                                                                              <  2:2156796/93‑1 (MQ=255)
            aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGGAACCTGGCGGTAGGCTGATTGTGATTGGCGTAATGTCTCCgg                                                                             >  2:3127329/1‑93 (MQ=255)
            aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCgg                                                                             >  2:3127330/1‑93 (MQ=255)
             tCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCggg                                                                            >  2:2312380/1‑93 (MQ=255)
               aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTc                                                                          <  1:1319049/93‑1 (MQ=255)
               aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTc                                                                          >  1:2653458/1‑93 (MQ=255)
               aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTc                                                                          <  2:1101223/93‑1 (MQ=255)
                 tGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAc                                                                        >  2:3530904/1‑93 (MQ=255)
                 tGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAc                                                                        >  1:3216271/1‑93 (MQ=255)
                   gCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCAccc                                                                      <  1:3200340/93‑1 (MQ=255)
                     atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAg                                                                    <  2:1145069/93‑1 (MQ=255)
                     atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAg                                                                    >  2:3997286/1‑93 (MQ=255)
                      tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGt                                                                   <  2:1012616/93‑1 (MQ=255)
                      tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGt                                                                   <  2:173783/93‑1 (MQ=255)
                       gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                                  <  2:4526919/93‑1 (MQ=255)
                         tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcg                                                                >  2:4721640/1‑93 (MQ=255)
                         tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcg                                                                <  2:88914/93‑1 (MQ=255)
                         tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcc                                                                >  1:2521039/1‑92 (MQ=255)
                          gTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTgcgc                                                               <  1:3790134/93‑1 (MQ=255)
                            tGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGa                                                             >  2:2315250/1‑93 (MQ=255)
                             ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                                            <  1:2455576/93‑1 (MQ=255)
                             ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                                            <  1:2105387/93‑1 (MQ=255)
                                gcagcaCTGCGTGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         >  1:4189701/1‑93 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  2:641411/93‑1 (MQ=255)
                                gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                                         <  2:3388323/93‑1 (MQ=255)
                                  agcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTg                                                       <  1:1378488/93‑1 (MQ=255)
                                  agcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTg                                                       <  1:4280523/93‑1 (MQ=255)
                                   gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                                      >  1:567314/1‑93 (MQ=255)
                                   gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                                      >  1:4693984/1‑93 (MQ=255)
                                      cTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTa                                                   >  1:667466/1‑93 (MQ=255)
                                       tGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                                  >  1:1154016/1‑93 (MQ=255)
                                        gCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                                 <  1:4476582/93‑1 (MQ=255)
                                           agaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGAc                                              >  1:1106323/1‑93 (MQ=255)
                                           agaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGAc                                              <  1:2173460/93‑1 (MQ=255)
                                           agaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGAc                                              >  2:3322137/1‑93 (MQ=255)
                                            gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             <  1:4412670/93‑1 (MQ=255)
                                            gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                             <  2:472163/93‑1 (MQ=255)
                                               gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          >  2:4175741/1‑93 (MQ=255)
                                               gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  1:353327/93‑1 (MQ=255)
                                               gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                                          <  1:4078261/93‑1 (MQ=255)
                                                tGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                                         >  1:3080819/1‑93 (MQ=255)
                                                   aTAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                                      >  2:4684381/1‑93 (MQ=255)
                                                    tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGc                                     <  2:180189/93‑1 (MQ=255)
                                                     aGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                                    >  2:4597407/1‑93 (MQ=255)
                                                     aGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                                    >  1:2844593/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   >  2:2759653/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   >  2:3020307/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   >  1:1674665/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                                   >  2:1983140/1‑93 (MQ=255)
                                                      ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTACAAGCAt                                   <  2:4667983/93‑1 (MQ=255)
                                                       gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCCtt                                  >  2:1237026/1‑93 (MQ=255)
                                                       gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                                  <  1:1961820/93‑1 (MQ=255)
                                                       gATATTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAAGGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                                  <  2:3198353/93‑1 (MQ=255)
                                                        atatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                                 >  1:2009490/1‑93 (MQ=255)
                                                        atatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                                 <  1:2567665/93‑1 (MQ=255)
                                                        atatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                                 <  2:4146200/93‑1 (MQ=255)
                                                         tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                >  2:2648461/1‑93 (MQ=255)
                                                         tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                <  2:2549054/93‑1 (MQ=255)
                                                         tatTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                                >  1:4269401/1‑93 (MQ=255)
                                                          atTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                               >  1:2156867/1‑93 (MQ=255)
                                                          atTGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                               >  2:3620637/1‑93 (MQ=255)
                                                           ttGAAACCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                              >  1:1722333/1‑93 (MQ=255)
                                                              aaaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt                           <  2:1014733/93‑1 (MQ=255)
                                                               aaCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa                          <  1:373192/93‑1 (MQ=255)
                                                                aCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc                         <  2:2475677/93‑1 (MQ=255)
                                                                aCCTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc                         >  1:2439493/1‑93 (MQ=255)
                                                                  cTGGCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct                       <  2:4365849/93‑1 (MQ=255)
                                                                    ggCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct                     >  2:932003/1‑93 (MQ=255)
                                                                    ggCGGTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct                     <  1:1845846/93‑1 (MQ=255)
                                                                       ggTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCAc                  >  1:3213150/1‑93 (MQ=255)
                                                                        gTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg                 <  2:3715643/93‑1 (MQ=255)
                                                                        gTAGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg                 >  1:2409635/1‑93 (MQ=255)
                                                                          aGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAGGCGATACCTCTCACGTa               >  2:1463586/1‑93 (MQ=255)
                                                                          aGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTa               >  2:3328260/1‑93 (MQ=255)
                                                                          aGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTa               <  2:721844/93‑1 (MQ=255)
                                                                          aGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTa               >  2:2981471/1‑93 (MQ=255)
                                                                          aGGCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTa               <  1:3360743/93‑1 (MQ=255)
                                                                           ggCTGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAc              >  1:4186467/1‑93 (MQ=255)
                                                                            gCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAca             >  2:1587648/1‑92 (MQ=255)
                                                                              tGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaa           >  1:1713566/1‑93 (MQ=255)
                                                                              tGATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaa           >  2:521814/1‑93 (MQ=255)
                                                                               gATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaaa          >  1:2207163/1‑93 (MQ=255)
                                                                               gATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaaa          <  1:3257223/93‑1 (MQ=255)
                                                                               gATTGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGaaa          <  1:1416613/93‑1 (MQ=255)
                                                                                 ttGTGATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACt        >  1:772174/1‑93 (MQ=255)
                                                                                   gtgATGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTAc      >  2:3609624/1‑93 (MQ=255)
                                                                                      aTGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGcc   <  1:3841679/93‑1 (MQ=255)
                                                                                       tGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCa  <  2:2261828/93‑1 (MQ=255)
                                                                                       tGGGCGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCa  >  2:122477/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTACGAAACTACGCCA  >  NC_000913/232128‑232307

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