Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,109,928 (A)7→6 100% intergenic (+105/‑2) sbp → / → cdh sulfate transporter subunit/CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,109,9220A.59.6% 50.1 / 199.6 104intergenic (+99/‑8)sbp/cdhsulfate transporter subunit/CDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand):  ref base A (16/26);  new base . (34/28);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

Only 100 of 114 total aligned reads displayed.

ATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGAAGA  >  NC_000913/4109823‑4110018
                                                                                                   |                                                                                                
aTTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGa                                                                                                  >  2:589722/1‑100 (MQ=255)
    ccccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaa                                                                                              >  2:1345326/1‑100 (MQ=255)
    ccccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaa                                                                                              >  1:2139323/1‑100 (MQ=255)
    ccccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaa                                                                                              >  1:2786782/1‑100 (MQ=255)
     cccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCATTACCCTTGCATCTCTTTGAGGTACAGGGaaaaa                                                                                              >  1:2343882/1‑99 (MQ=255)
     ccaaGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGagaaaa                                                                                             >  1:2299007/4‑100 (MQ=255)
      ccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaa                                                                                              >  2:1014146/1‑98 (MQ=255)
       cAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaaag                                                                                            >  2:602562/1‑98 (MQ=255)
         ggATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGAt                                                                                         <  2:1504447/100‑1 (MQ=255)
          gATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGa                                                                                       <  2:2893456/100‑1 (MQ=255)
          gATAATTATCAAACCCGGTGGTTTCTCACGACCAGGTTTTTTATTTGTCACGTTTTAGCGATACCCTTGCGTCTCTTTGAGGCACAGGG‑GAATAAGATGa                                                                                       >  1:2420153/1‑100 (MQ=255)
           aTAATTATCAAACCCGGTGGTTTCTCGGGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGa                                                                                       <  2:1182670/100‑1 (MQ=255)
             aaTTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGaaaa                                                                                    >  1:1873389/1‑100 (MQ=255)
               ataaaaaacccggtcgcgagaaaccaccgggtttTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGaaaaaa                                                                                  >  1:1739833‑M2/35‑100 (MQ=255)
               agaaaaaacccggtcgcgagaaaccaccgggtttTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGaaaaaa                                                                                  >  2:1962816‑M2/35‑100 (MQ=255)
                 aaaaaacccggtcgcgagaaaccaccgggtttTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAg                                                                                 >  1:779032‑M2/33‑100 (MQ=255)
                  tCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCg                                                                               >  2:498002/1‑100 (MQ=255)
                  aaaaacccggtcgcgagaaaccaccgggtttTTTATTTGTCACGATTTTGCGTTACCCATGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGc                                                                                >  1:1475044‑M2/32‑100 (MQ=255)
                   cAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCgg                                                                              <  2:2240349/100‑1 (MQ=255)
                    aaaCACGATGGTTTCGCGCGACCGGGTTTTTTATTTGGCACGATTTTTCGTTACACTTGCATCTCTTTGAGGTACAGGGAAAAAAA‑AAAAAAAAAGCggg                                                                             <  1:1706164/100‑1 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             <  1:62582‑M1/100‑1 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             <  1:2131327‑M1/100‑1 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             >  2:2710500‑M1/1‑100 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             <  1:1847739‑M1/100‑1 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             >  2:612235‑M1/1‑100 (MQ=255)
                    AAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCggg                                                                             >  2:62583‑M1/1‑100 (MQ=255)
                     AACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCGGGt                                                                            <  2:820424‑M1/100‑1 (MQ=255)
                       CCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCGGGtct                                                                          <  2:2193093‑M1/100‑1 (MQ=255)
                         CGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCGGGtcttc                                                                        <  2:84324‑M1/100‑1 (MQ=255)
                         CGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAA‑GATGAAAAAAGCGGGtcttc                                                                        <  2:2556445‑M1/100‑1 (MQ=255)
                                tgctcGCGACCGGTTTTTTTATTTGTCACGATGTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCttttttt                                                                 <  1:2306391/98‑1 (MQ=255)
                                  ctcGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCGTGCATTTCTTTGAGATACAAAGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTg                                                                <  1:1080524/100‑1 (MQ=255)
                                   tcGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGt                                                              <  2:2524277/100‑1 (MQ=255)
                                    cgcgACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtg                                                             <  1:863773/100‑1 (MQ=255)
                                    cgcgACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtg                                                             <  2:767567/100‑1 (MQ=255)
                                    cgcgACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGt                                                              >  1:931352/1‑99 (MQ=255)
                                     gcgACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtg                                                             >  1:1267896/1‑99 (MQ=255)
                                          cgggtttTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATa                                                       >  2:2107338‑M2/8‑100 (MQ=255)
                                          cgggtttTTTATTTGTAACGATTTTGCGTAACCCTTGCACCTCTTTGAGTTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTg                                                                <  1:1265785‑M2/85‑1 (MQ=255)
                                           gggtttTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAg                                                      <  1:1511778‑M2/94‑1 (MQ=255)
                                            ggtttTTTATTTGTGACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTA‑AAAAAAAAAAAGATAAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGt                                                     <  1:126740‑M2/95‑1 (MQ=255)
                                             gtttTTTATTTGTCACGATTTTGCGTAACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATagn                                                     <  1:1043932‑M2/96‑2 (MQ=255)
                                                tTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGAGGGTCTTCTTTTTTTGGGGCTGATCGt                                                     >  2:2617368‑M2/2‑96 (MQ=255)
                                                 tttATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCg                                                >  2:1754708/1‑100 (MQ=255)
                                                 tttATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCg                                                >  2:556674/1‑100 (MQ=255)
                                                  ttATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGc                                               <  2:213877/100‑1 (MQ=255)
                                                   tATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGc                                               <  1:112318/100‑1 (MQ=255)
                                                   tATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGcc                                              >  1:373763/1‑100 (MQ=255)
                                                      ttGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGtt                                           <  2:932534/100‑1 (MQ=255)
                                                       tGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTtg                                          <  1:1612158/100‑1 (MQ=255)
                                                       tGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTtg                                          <  1:23103/100‑1 (MQ=255)
                                                         tCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTAGAGGG‑AAATAAGATGAAAAAAGCGGGGCTTCTTTTTTTGGTGATGATAGTTATCGCCGTtgtg                                        >  1:2318279/1‑100 (MQ=255)
                                                         tCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTtgtg                                        >  2:1630070/1‑100 (MQ=255)
                                                          cACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTTTTTATTTTAGTTATCGCCGTtgtg                                        >  1:545423/1‑100 (MQ=255)
                                                          cACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTgg                                       >  1:1154621/1‑100 (MQ=255)
                                                          cACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTgg                                       >  2:242016/1‑100 (MQ=255)
                                                          cACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTgg                                       >  2:1081779/1‑100 (MQ=255)
                                                           aCGATTTTGCGTTACGCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTgg                                       <  1:86524/100‑1 (MQ=255)
                                                           aCGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGc                                      >  1:795359/1‑100 (MQ=255)
                                                            cGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCt                                     >  2:1755875/1‑100 (MQ=255)
                                                              aTTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTg                                    >  1:416053/1‑100 (MQ=255)
                                                                 ttGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCgg                                <  2:1081168/100‑1 (MQ=255)
                                                                 ttGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCgg                                <  2:332410/100‑1 (MQ=255)
                                                                 ttGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCgg                                <  1:2107331/100‑1 (MQ=255)
                                                                  tGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCgg                                <  2:100919/100‑1 (MQ=255)
                                                                  tGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCgg                                >  1:1986247/1‑100 (MQ=255)
                                                                  tGCGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGt                               <  2:1506350/100‑1 (MQ=255)
                                                                    cGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTa                              >  1:261311/1‑100 (MQ=255)
                                                                    cGTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTAt                             >  2:288522/1‑100 (MQ=255)
                                                                     gTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTAt                             >  2:836788/1‑100 (MQ=255)
                                                                     gTTACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTAtt                            <  1:23131/100‑1 (MQ=255)
                                                                      ttACCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTg                           >  2:1682779/1‑100 (MQ=255)
                                                                       tACCCGTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTg                           <  2:356737/100‑1 (MQ=255)
                                                                        aCCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTgg                          <  1:165509/100‑1 (MQ=255)
                                                                        aCCCTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGt                         >  1:1517138/1‑100 (MQ=255)
                                                                         cccTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGtt                        >  2:2253023/1‑100 (MQ=255)
                                                                         cccTTGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGtt                        <  2:2101822/100‑1 (MQ=255)
                                                                            ttGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACt                     <  1:1753600/100‑1 (MQ=255)
                                                                             tGCATCTCTTTGAGGTACAGGGAAAAAAA‑AAGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTg                    <  1:172371/100‑1 (MQ=255)
                                                                             tGCATCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTg                    >  2:172372/1‑100 (MQ=255)
                                                                              gCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATTGTTATCGCCGTTGTGGTTGCCGGTATTGGTTACt                     >  1:608848/1‑99 (MQ=255)
                                                                              gCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGTCTGCCGGTATTGGTTACTg                    >  2:2690698/1‑100 (MQ=255)
                                                                                aTCTCTTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGaa                 <  2:2624394/100‑1 (MQ=255)
                                                                                  ctctTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAAt               <  1:2199748/100‑1 (MQ=255)
                                                                                   tctTTGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAAt               <  1:2253013/99‑1 (MQ=255)
                                                                                     tttGAGGTACAGTGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAAt               <  2:2719029/98‑1 (MQ=255)
                                                                                      ttGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTaa            >  1:2115542/1‑100 (MQ=255)
                                                                                      ttGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAAc           <  1:1433176/100‑1 (MQ=255)
                                                                                      ttGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAAc           >  1:2340417/1‑100 (MQ=255)
                                                                                       tGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAAcc          >  1:1510613/1‑100 (MQ=255)
                                                                                       tGAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTAGCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAAcc          <  2:387811/100‑1 (MQ=255)
                                                                                        gAGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCg         >  2:229532/1‑100 (MQ=255)
                                                                                         aGGTACAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCgg        <  2:493690/100‑1 (MQ=255)
                                                                                          ggTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCgg        >  1:2815977/1‑100 (MQ=255)
                                                                                             aCAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAGATTAACCGGTGa     <  2:1913702/100‑1 (MQ=255)
                                                                                             aCAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCACCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGaa    >  1:1926067/1‑100 (MQ=255)
                                                                                             aCAGGG‑AAAAAAGATGAAAAAAGCGGGTCTACTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTCaa    >  1:1065456/1‑100 (MQ=255)
                                                                                              cAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGAag   >  2:8101/1‑100 (MQ=255)
                                                                                              cAGGG‑AAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGAag   >  2:806320/1‑100 (MQ=255)
                                                                                                gggAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGAaga  <  2:2815987/100‑1 (MQ=255)
                                                                                                   |                                                                                                
ATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGATGATAGTTATCGCCGTTGTGGCTGCCGGTATTGGTTACTGGAAATTAACCGGTGAAGA  >  NC_000913/4109823‑4110018

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