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→TAA) yaiSputative 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

CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCA  >  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)
                         aGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑aa                                                                                                             >  2:40515/1‑139 (MQ=255)
                               aTTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGa                                                                                                       <  1:1359045/139‑1 (MQ=255)
                               aTTTTATGAATTACAAAGCCCTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGa                                                                                                       <  2:1303506/139‑1 (MQ=255)
                                 tttATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACt                                                                                                     <  1:866812/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACtgtg                                                                                                  <  2:989171/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACtgtg                                                                                                  <  2:415038/139‑1 (MQ=255)
                                    aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTCAATGCAGC‑AAGTTTGACtgtg                                                                                                  <  2:198450/139‑1 (MQ=255)
                                     tGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGa                                                                                                 >  1:1048649/1‑139 (MQ=255)
                                           aCAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATAtt                                                                                           >  2:371356/1‑139 (MQ=255)
                                           aCAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGGTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATAtt                                                                                           >  1:1130056/1‑139 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACt                                                                                                     <  1:229577/125‑1 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACt                                                                                                     >  2:229577/1‑125 (MQ=255)
                                               aGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATAttcttc                                                                                       >  2:851930/1‑139 (MQ=255)
                                                 ccTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTcttctt                                                                                     >  2:82716/1‑139 (MQ=255)
                                                    tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑aa                                                                                 >  1:365241/1‑139 (MQ=255)
                                                    tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑aa                                                                                 >  2:72383/1‑139 (MQ=255)
                                                       cccAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACgg                                                                              <  1:40515/139‑1 (MQ=255)
                                                       cccAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACgg                                                                              <  1:371356/139‑1 (MQ=255)
                                                           gATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTc                                                                          >  2:1207565/1‑139 (MQ=255)
                                                            aTCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCa                                                                         >  1:935714/1‑139 (MQ=255)
                                                              cGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCaaa                                                                       <  2:213880/139‑1 (MQ=255)
                                                                 tATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                                <  1:66709/97‑1 (MQ=255)
                                                                 tATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                                >  2:66709/1‑97 (MQ=255)
                                                                  aTTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAAcc                                                                   >  2:1102803/1‑139 (MQ=255)
                                                                   ttGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCt                                                                  <  1:82716/139‑1 (MQ=255)
                                                                   ttGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCt                                                                  <  1:1102803/139‑1 (MQ=255)
                                                                    tGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTg                                                                 >  1:163096/1‑139 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCAt                                                           >  1:203230/1‑139 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCACC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAATACCTGAGGTAt                                                           <  1:851930/139‑1 (MQ=255)
                                                                          gttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCAt                                                           <  2:1048649/139‑1 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTg                                                                 <  1:1429524/132‑1 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTg                                                                 >  2:1429524/1‑132 (MQ=255)
                                                                           ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTT‑ATTCAGCAAATTTTGATGGTG‑AATTTTCTTCTTTAAACGAATTAAAAAACTTAAGGctaa                                                          >  2:1258181/1‑136 (MQ=255)
                                                                               ccGCGAAATTGTGCTACTGCACGCAGGCGATCATTACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATg                                                      <  2:402752/139‑1 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                    >  2:1341939/1‑78 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                    <  1:1341939/78‑1 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGa                                                     >  2:1302415/1‑139 (MQ=255)
                                                                                cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGa                                                     <  1:1207565/139‑1 (MQ=255)
                                                                                      aTTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCa                                               <  2:995226/139‑1 (MQ=255)
                                                                                                  cACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTCTGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGt                                   >  1:948577/1‑139 (MQ=255)
                                                                                                    cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGcc                                                >  1:774299/1‑124 (MQ=255)
                                                                                                    cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGcc                                                <  2:774299/124‑1 (MQ=255)
                                                                                                       aGGCGATCATGAGGCATATCATCGCGTCGCACTTAGCTTTTATATTTTTTTAATTCAGC‑AAGTTTGACTGTGAAATATTCTTCTCT‑TACGGATTCAAAAACCTGAGGCATAAATGAAAGTCACGTatatggagcttcat                              >  1:621656/1‑126 (MQ=255)
                                                                                                         gCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAg                            >  2:1454689/1‑139 (MQ=255)
                                                                                                                 gACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACATGAGGCATAAATGAAAGCCAGGTACTTGGGGTTACGTAGCCGAGAAt                    >  2:1349008/1‑139 (MQ=255)
                                                                                                                   cGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAttt                  >  2:1129382/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGT‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:962288/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGa                                                     <  1:1086940/100‑1 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGa                                                     >  2:1086940/1‑100 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:91023/1‑139 (MQ=255)
                                                                                                                       tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCgg              >  2:867155/1‑139 (MQ=255)
                                                                                                                        ataATCGCGACGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGa             >  1:1463768/1‑139 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGaa                                                                                                >  2:1433249/1‑56 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGaa                                                                                                <  1:1433249/56‑1 (MQ=255)
                                                                                                                         taATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTTGGGTTTCGTAGCCGAGAATTTGCGGaa            <  1:1129382/139‑1 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAtt                   >  2:1456709/1‑126 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAAtt                   <  1:1456709/126‑1 (MQ=255)
                                                                                                                               cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTc      >  1:839026/1‑139 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCgaga                      >  2:784989/1‑122 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCgaga                      <  1:784989/122‑1 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCg     >  2:974808/1‑139 (MQ=255)
                                                                                                                                gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCg     >  2:378152/1‑139 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACtgtg                                                                                                  >  1:515710/1‑46 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACtgtg                                                                                                  <  2:515710/46‑1 (MQ=255)
                                                                                                                                 tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCgg    <  2:365241/139‑1 (MQ=255)
                                                                                                                                  cgtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggc   >  1:740838/1‑139 (MQ=255)
                                                                                                                                   gtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggca  <  2:740838/139‑1 (MQ=255)
                                                                                                                                   gtTCTTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggca  >  2:664542/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGC‑AAGTTTGACTGTGAAATATTCTTCTTT‑AACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCA  >  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: ‑