Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 384,059 C→T *186* (TAG→TAA yaiS ← putative PIG‑L family deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913384,0590CT97.3% 214.2 / ‑2.0 73*186* (TAG→TAAyaiSputative PIG‑L family deacetylase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (43/28);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCA  >  NC_000913/383924‑384193
                                                                                                                                       |                                                                                                                                        
ccAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAg                                                                                                                                       >  2:1359045/1‑139 (MQ=255)
  aGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACACAAATCGCTATTAACCTGTTGTCTGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGAATATAATCGCGTCTTTCTTAGCt                                                                                                                                     >  1:402752/1‑139 (MQ=255)
    gTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCttt                                                                                                                                   >  1:995226/1‑139 (MQ=255)
       ttCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTAt                                                                                                                                >  2:866812/1‑139 (MQ=255)
         cATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGt                                                                                                                              >  1:989171/1‑139 (MQ=255)
                         aGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCaa                                                                                                             >  2:40515/1‑139 (MQ=255)
                               aTTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGa                                                                                                       <  1:1359045/139‑1 (MQ=255)
                               aTTTTATGAATTACAAAGCCCTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGa                                                                                                       <  2:1303506/139‑1 (MQ=255)
                                 tttATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                     <  1:866812/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                                  <  2:989171/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                                  <  2:415038/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTCAATGCAGCAAGTTTGACtgtg                                                                                                  <  2:198450/139‑1 (MQ=255)
                                     tGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGa                                                                                                 >  1:1048649/1‑139 (MQ=255)
                                           aCAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAtt                                                                                           >  2:371356/1‑139 (MQ=255)
                                           aCAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGGTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAtt                                                                                           >  1:1130056/1‑139 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                     <  1:229577/125‑1 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                     >  2:229577/1‑125 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttcttc                                                                                       >  2:851930/1‑139 (MQ=255)
                                                 ccTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTcttctt                                                                                     >  2:82716/1‑139 (MQ=255)
                                                    tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTaa                                                                                 >  1:365241/1‑139 (MQ=255)
                                                    tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTaa                                                                                 >  2:72383/1‑139 (MQ=255)
                                                       cccAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACgg                                                                              <  1:40515/139‑1 (MQ=255)
                                                       cccAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACgg                                                                              <  1:371356/139‑1 (MQ=255)
                                                           gATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTc                                                                          >  2:1207565/1‑139 (MQ=255)
                                                            aTCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCa                                                                         >  1:935714/1‑139 (MQ=255)
                                                              cGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCaaa                                                                       <  2:213880/139‑1 (MQ=255)
                                                                 tATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                                <  1:66709/97‑1 (MQ=255)
                                                                 tATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                                >  2:66709/1‑97 (MQ=255)
                                                                  aTTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAAcc                                                                   >  2:1102803/1‑139 (MQ=255)
                                                                   ttGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCt                                                                  <  1:82716/139‑1 (MQ=255)
                                                                   ttGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCt                                                                  <  1:1102803/139‑1 (MQ=255)
                                                                    tGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTg                                                                 >  1:163096/1‑139 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCAt                                                           >  1:203230/1‑139 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCACCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAATACCTGAGGTAt                                                           <  1:851930/139‑1 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCAt                                                           <  2:1048649/139‑1 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTg                                                                 <  1:1429524/132‑1 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTg                                                                 >  2:1429524/1‑132 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTATTCAGCAAATTTTGATGGTGAATTTTCTTCTTTAAACGAATTAAAAAACTTAAGGctaa                                                          >  2:1258181/1‑136 (MQ=255)
                                                                               ccGCGAAATTGTGCTACTGCACGCAGGCGATCATTACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATg                                                      <  2:402752/139‑1 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                    >  2:1341939/1‑78 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                    <  1:1341939/78‑1 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGa                                                     >  2:1302415/1‑139 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGa                                                     <  1:1207565/139‑1 (MQ=255)
                                                                                      aTTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCa                                               <  2:995226/139‑1 (MQ=255)
                                                                                                  cACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTCTGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGt                                   >  1:948577/1‑139 (MQ=255)
                                                                                                    cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGcc                                                >  1:774299/1‑124 (MQ=255)
                                                                                                    cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGcc                                                <  2:774299/124‑1 (MQ=255)
                                                                                                       aGGCGATCATGAGGCATATCATCGCGTCGCACTTAGCTTTTATATTTTTTTAATTCAGCAAGTTTGACTGTGAAATATTCTTCTCTTACGGATTCAAAAACCTGAGGCATAAATGAAAGTCACGTatatggagcttcat                              >  1:621656/1‑126 (MQ=255)
                                                                                                         gCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAg                            >  2:1454689/1‑139 (MQ=255)
                                                                                                                 gACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACATGAGGCATAAATGAAAGCCAGGTACTTGGGGTTACGTAGCCGAGAAt                    >  2:1349008/1‑139 (MQ=255)
                                                                                                                   cGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAttt                  >  2:1129382/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGTAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:962288/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGa                                                     <  1:1086940/100‑1 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGa                                                     >  2:1086940/1‑100 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:91023/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:867155/1‑139 (MQ=255)
                                                                                                                        ataATCGCGACGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGa             >  1:1463768/1‑139 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGaa                                                                                                >  2:1433249/1‑56 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGaa                                                                                                <  1:1433249/56‑1 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTTGGGTTTCGTAGCCGAGAATTTGCGGaa            <  1:1129382/139‑1 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAtt                   >  2:1456709/1‑126 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAtt                   <  1:1456709/126‑1 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTc      >  1:839026/1‑139 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCgaga                      >  2:784989/1‑122 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCgaga                      <  1:784989/122‑1 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCg     >  2:974808/1‑139 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCg     >  2:378152/1‑139 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                                  >  1:515710/1‑46 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                                  <  2:515710/46‑1 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCgg    <  2:365241/139‑1 (MQ=255)
                                                                                                                                  cgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggc   >  1:740838/1‑139 (MQ=255)
                                                                                                                                   gtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggca  <  2:740838/139‑1 (MQ=255)
                                                                                                                                   gtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggca  >  2:664542/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCA  >  NC_000913/383924‑384193

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