Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,958 G→T 5.5% intergenic (‑45/+28) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9580GT 5.5% 679.8 / 20.5 194intergenic (‑45/+28)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (94/89);  new base T (6/5);  total (100/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-03

Only 100 of 195 total aligned reads displayed.

AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCA  >  NC_000913/4470684‑4471227
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                             
       tAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                                                                                                                                                                                                        <  2:542019/275‑1 (MQ=255)
                    gTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCTCGTTTTCGGATTTACCGCGATCTCACCGGATGTGATGCTAATATTACCCAGCTCCACGCACTGTCCGTAAGTACCGTTCACTGCCGGCGCCTTTTCCGGCCATATATCTTTACTCATTATTTCACCTTTATTACCACTGTCACAACTCTTGACCATTATAGGGCACCACACtttacttccca                                                                                                                                                                                                                                                           >  4:86477/1‑267 (MQ=255)
                            cTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                                               <  1:53273/279‑1 (MQ=255)
                                cgggCCCGTGGTTTCGACGGTCTCTTTTTCTTTNACCGGTGCCTTGCGTGCCTGTTCCGNGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCNCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                                           <  3:135036/277‑1 (MQ=255)
                                       gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                                                                                                                                                                                                         <  2:617417/274‑1 (MQ=255)
                                                       ctgttACGTTGTCCAGCGAGTAACGGTGCTGTGCAGGCAGGTATGCCGGTACTTCGCCCGATTACGGCATTACCGGGATCTGACCGGGGGTGATGATCAGATTGCACAGCTCAATGCTCTGTACGTAAGGACCGATGGCTGCCGGTGCATTTTCCGGGGCGATGGTTTTGCTCATGATTTCTCCTGTATTACAGCGGTAAAAAGTCTGGCTCATTGTAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                                                                                                                                                                                                        <  4:677326/272‑1 (MQ=255)
                                                        ttttACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaa                                                                                                                                                                                                                        <  2:557524/274‑1 (MQ=255)
                                                                   ccAGCGACTGACTGGCCGTGGCGACAACTTCTGCCGGTACTTCGCCCGGTTTTGGATTTACCGGGATCTGCCCGCAGGTGATGATCATATTGCCCACATGAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGAACCACATTATGGTAAAACTCTTTTTca                                                                                                                                                                                                            <  4:236586/275‑1 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                            <  2:324415/275‑1 (MQ=255)
                                                                   ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                                                                                                                                                                                            <  2:611503/275‑1 (MQ=255)
                                                                          cTGACTGGCCTTTGCAGCGAAGGTTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                                                                                                                                                                                                     <  2:662024/275‑1 (MQ=255)
                                                                                    tgtgccggcaggTCTGCCGGTACTTCGCCCCGTTTCGGATTTACCGGGATCTGACCGGAGGTGATGACCATATTGCCGAGATCAACGCCCTGTACGTAAGGACCGATAGCTACCGGTGCAATGTCCGTCGAGATAGTTTTGCTCATGATTTCTCTTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCAAATTATGGGAAAACTCTTTTTCACAGTATGTGCATTTGAg                                                                                                                                                                                           <  4:97811/264‑1 (MQ=255)
                                                                                                tCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTATGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCCACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAATCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgc                                                                                                                                                                            >  1:21506/1‑278 (MQ=255)
                                                                                                    ccGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                                                                                                                                                                                                        <  4:334635/182‑1 (MQ=255)
                                                                                                                  gTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTGTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                                                                                                                                                             <  2:6800/275‑1 (MQ=255)
                                                                                                                                 gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGGTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                                                                                                                                                                                             >  7:732043/1‑148 (MQ=255)
                                                                                                                                 gggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                                                                                                                                                                                            >  5:718655/1‑149 (MQ=255)
                                                                                                                                  ggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                                                                                                                                                                                            >  8:401458/1‑148 (MQ=255)
                                                                                                                                    aTCTGACCGGAGGTGATGATCATATTGCCCCGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGcc                                                                                                                                                                                                                                                                          <  8:124330/148‑1 (MQ=255)
                                                                                                                                    aTCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                           >  2:52237/1‑275 (MQ=255)
                                                                                                                                        gACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                                                                                                                                                                                                     >  5:845202/1‑149 (MQ=255)
                                                                                                                                          ccGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGc                                                                                                                                     <  4:423734/275‑1 (MQ=255)
                                                                                                                                            ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTTCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                                                                                                                                                                                                  <  8:811053/148‑1 (MQ=255)
                                                                                                                                              aGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                                                                                                                                                                                                <  7:907823/148‑1 (MQ=255)
                                                                                                                                                  gatgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                                                                                                                                                                                                <  4:270448/176‑1 (MQ=255)
                                                                                                                                                   atgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                                                                                                                                                                                           <  8:706242/148‑1 (MQ=255)
                                                                                                                                                   atgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGAt                                                                                                                            >  2:358084/1‑275 (MQ=255)
                                                                                                                                                    tgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                                                                                                                                                                                                          >  8:291986/1‑148 (MQ=255)
                                                                                                                                                     gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                                                                                                                                                                                                        >  6:265908/1‑149 (MQ=255)
                                                                                                                                                        cATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGACGGTGCATTTTCCGTCGCGACAGGTTTGCTAATGATTTCTCCTGTATTACAGCGGTAACAAGTCTGGCTTATTATAGGGATCCAGACTTTTATTACCAACCGCAATTTCTTGGACAGCACCACAGTATGCGAAAACTCTTTTTCACAGTATTTGCATTTCAGCGCGATATCATTGGACAGTTTTCGCACGGCAAAGCTTGATGAAACCGGTTCCGCATGTCTAATCCatt                                                                                                                       >  4:487082/1‑273 (MQ=255)
                                                                                                                                                          tatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                                                                                                                                                                                                   >  5:670249/1‑149 (MQ=255)
                                                                                                                                                           atTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                                                                                                                                                                                                  >  5:512666/1‑149 (MQ=255)
                                                                                                                                                               cccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                                                                                                                                                                                              <  5:851101/149‑1 (MQ=255)
                                                                                                                                                               cccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaat                                                                                                                                                                                                                                               >  7:650375/1‑148 (MQ=255)
                                                                                                                                                               cccAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGt                                                                                                                 <  2:171767/274‑1 (MQ=255)
                                                                                                                                                                  aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                                                                                                                                            <  8:919814/148‑1 (MQ=255)
                                                                                                                                                                  aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                                           <  6:653416/149‑1 (MQ=255)
                                                                                                                                                                  aGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGc                                                                                                                                                                                                                                           <  6:845202/149‑1 (MQ=255)
                                                                                                                                                                        aagCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGca                                                                                                                                                                                                                                      >  8:861454/3‑148 (MQ=255)
                                                                                                                                                                         cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACAGGTTCGGCATGGCTGTTACAGTTGCTGTTCGGGTAGAcc                                                                                                      >  2:655760/1‑275 (MQ=255)
                                                                                                                                                                         cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTAGCGATAGTATTGCTCATGATTTCTCCTTTATTACAGCGGTACAACGTCTGGCACATTCTAGGGAGCCAGAATTTTATTACCATCCGGAATTACTTGGCCAGCACCACTTTATGGGAAAACTCTTTTTCAAAGTCTTTGCATTTGAGCGCGATATCATTGGCCCGTTTTCGAACGGCCAAGCTGGATGAAACCGTTTAGGCTTGGCTGCTACAGTTGCGTTTCGGGCAGAcc                                                                                                      >  4:154510/1‑275 (MQ=255)
                                                                                                                                                                             cTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCCACCGCAATTAATTGGCCCGCCCCCCAt                                                                                                                                                                                                                                >  5:983418/1‑149 (MQ=255)
                                                                                                                                                                              tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACcaca                                                                                                                                                                                                                                 >  4:753115/1‑147 (MQ=255)
                                                                                                                                                                              tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATa                                                                                                                           <  2:588451/249‑1 (MQ=255)
                                                                                                                                                                                tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCTGGCAGACCAGCACCTTGTc                                                                                           >  1:72489/1‑279 (MQ=255)
                                                                                                                                                                                 aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                             <  7:960203/148‑1 (MQ=255)
                                                                                                                                                                                  cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                                                                                                                                                                                                             <  2:753115/147‑1 (MQ=255)
                                                                                                                                                                                  cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                                                                                                                                                            <  7:401458/148‑1 (MQ=255)
                                                                                                                                                                                  cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                                                                                                                                                            <  8:524086/148‑1 (MQ=255)
                                                                                                                                                                                    tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                                                                                                                                                                                                         <  6:959877/149‑1 (MQ=255)
                                                                                                                                                                                      aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                                                                                                                                                                                                            >  4:593030/1‑144 (MQ=255)
                                                                                                                                                                                      aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                                                                                                                                                       >  5:650161/1‑149 (MQ=255)
                                                                                                                                                                                          ccGATAGTTGCCGGGGCATTTTCCGTCTCGTTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGCCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATgcg                                                                                     <  4:439206/275‑1 (MQ=255)
                                                                                                                                                                                          ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTCTTACAGCGGTAAAAAGTCTGGCGCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAATCTTTTTCACAGTATTTGAATTTGAGCGCGATATCATTTGCGCGTTTTCGCCCGGCAACGCAGGATGAAACCGGTTCGGCATGGCTGATACCGGTGNNNTTAGTGCTGCCCAGACCATTTTCTATGCGCGcc                                                                                 >  3:225235/1‑279 (MQ=255)
                                                                                                                                                                                          ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTTTTCGGGCAGACCAGCACATTGGCGTTgcc                                                                                     >  4:549011/1‑274 (MQ=255)
                                                                                                                                                                                            gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCt                                                                                   >  2:45378/1‑275 (MQ=255)
                                                                                                                                                                                             aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                                                                                                                                                                                                <  6:582590/149‑1 (MQ=255)
                                                                                                                                                                                                gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGt                                                                                                                                          <  2:97173/216‑1 (MQ=255)
                                                                                                                                                                                                   gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAACAGTCTGGCTCATTATAGGGTGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                                                                                                                                                                                           >  8:32127/1‑148 (MQ=255)
                                                                                                                                                                                                    ccGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                                                                                                                                                                                                         >  5:372553/1‑149 (MQ=255)
                                                                                                                                                                                                       gTGCATTTTCCGTCGCGATAGTTTTGCTCATGTTTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGTCTTTTATAACCAACCGCAATTAATTGGCTAGCACCACATTATGGGAAAACTCTTTTTCACAGTa                                                                                                                                                                                                       <  7:861454/148‑1 (MQ=255)
                                                                                                                                                                                                            ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAAGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGgc                                                                    >  2:120913/1‑274 (MQ=255)
                                                                                                                                                                                                             tttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                                                                                                                                                                                                 <  7:291986/148‑1 (MQ=255)
                                                                                                                                                                                                              ttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                                                                                                                                                                                               >  6:695746/1‑149 (MQ=255)
                                                                                                                                                                                                                ccGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                                                                                                                                                                                             >  5:577110/1‑149 (MQ=255)
                                                                                                                                                                                                                   tCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                                                                                                                                                                                           >  8:563516/1‑148 (MQ=255)
                                                                                                                                                                                                                    cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc                                                                                                                                                                                          >  8:185544/1‑148 (MQ=255)
                                                                                                                                                                                                                      ctttttttttGTTCATGATTTCTCTTTTTTGANGNGGGTAAAAGGTTTGGCTAATTATCGGGAGCCAGACTTTTATTACANACCGCAATTAATTGACCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGANATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccacca                                                     <  3:72489/273‑1 (MQ=255)
                                                                                                                                                                                                                       gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                                                                                                                                                       <  7:570132/148‑1 (MQ=255)
                                                                                                                                                                                                                        aTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGat                                                                                                                                                                                     >  6:447770/1‑149 (MQ=255)
                                                                                                                                                                                                                         tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACccac                                                       <  2:407328/274‑1 (MQ=255)
                                                                                                                                                                                                                              ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaa                                                                                                                                                                                                                                                <  1:117660/84‑1 (MQ=255)
                                                                                                                                                                                                                                 cataTGATTTCTCCTTTATTACAGCGATAAAAATTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATa                                              <  4:655760/272‑1 (MQ=255)
                                                                                                                                                                                                                                 cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGtt                                                                                                                >  4:402051/1‑209 (MQ=255)
                                                                                                                                                                                                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                                                                                                                                                           >  8:624062/1‑148 (MQ=255)
                                                                                                                                                                                                                                    aTGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGcgcg                                                                                                                                                                         >  6:906167/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATgc                                                                                      >  2:580648/1‑231 (MQ=255)
                                                                                                                                                                                                                                       aTTTCTCCTTTATTACAGCGGTAAAAAGGCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTGTCACAGTATTGGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGTGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGAt                                     <  1:155872/278‑1 (MQ=255)
                                                                                                                                                                                                                                         ttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTc                                                                                  <  4:580648/231‑1 (MQ=255)
                                                                                                                                                                                                                                              cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGAAGACCAGCACATTGTCGAGGCGCTGCGGCAGACTGGGGCGCGGTTTACCCACCACTTCATAGTTTTcactccgg                                 >  4:346905/1‑268 (MQ=255)
                                                                                                                                                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                                                                                                                                              >  6:1077013/1‑149 (MQ=255)
                                                                                                                                                                                                                                                ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGGGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATACCATTGGCGCGTTTTCGCACgg                                                                                                                                                              >  8:766186/1‑148 (MQ=255)
                                                                                                                                                                                                                                                ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                                                                                                                                              >  7:809488/1‑148 (MQ=255)
                                                                                                                                                                                                                                                ttattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACgg                                                                                                                                                              >  8:236937/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                                                                                                                                                                                                         >  2:399755/1‑65 (MQ=255)
                                                                                                                                                                                                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                     <  5:260213/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGa                                                                                                                                                    <  6:790111/147‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGa                                                                                                                                                    >  5:790111/1‑147 (MQ=255)
                                                                                                                                                                                                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCtt                    >  2:235815/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             aaaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGa                                                                                                                                                 >  8:15414/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                               >  5:386109/1‑82 (MQ=25)
                                                                                                                                                                                                                                                                 aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                                                                                                                                                                                                               <  6:386109/82‑1 (MQ=37)
                                                                                                                                                                                                                                                                  gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCgg                                                                                                                                           >  6:558551/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGtt                                                                                                                                         <  7:789016/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAACACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg                                                                                                                                      >  7:231242/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg                                                                                                                                 <  6:162462/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATgg                                                                                                                                 >  5:162462/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                                                                                                                                                                                              >  6:303690/1‑68 (MQ=255)
                                                                                                                                                                                                                                                                                aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttttt                                                                                                                                                                                                              <  5:303690/68‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGa                                                                                                                             <  5:906167/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                             
AACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCA  >  NC_000913/4470684‑4471227

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