Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,875,5320A→T19.5% 12.8 ‑9.8 ‑0.077H32L (CAT→CTT) yidXputative lipoprotein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (15/0):  ref base (27/35):  total (42/35)
Fisher's exact test strand distribution p-value = 5.65e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

ATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGG  >  NC_000913/3875353‑3875726
                                                                                                                                                                                   |                                                                                                                                                                                                  
aTCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTAAg                                                                                                                                                                                >  2:242788/1‑200 (MQ=255)
                   gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTc                                                                                                                                                             >  2:19992/1‑200 (MQ=255)
                      cTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTTGAAAACAAAAAGTAGTTCAGGAAAGTATCAATCAGAGTCAGa                                                                                                                                                          >  2:751256/1‑200 (MQ=255)
                        aaaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTCCCTTAAAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGaca                                                                                                                                                        >  2:683447/1‑200 (MQ=255)
                              gaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCa                                                                                                                                                  >  2:154465/1‑200 (MQ=255)
                               aggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGAAGTTCAGGAAAGTATAAATCGCAGTCAGACACTATCAc                                                                                                                                                 >  2:767911/1‑200 (MQ=255)
                                 ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTCTCCCTGACAGCAAAGAATAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCAcgg                                                                                                                                               >  2:190739/1‑198 (MQ=255)
                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTgg                                                                                                                                              >  2:25741/1‑200 (MQ=255)
                                  gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAGAAAAAATAAGATCAGGAAAGTATCAAACGCAATCAGTCAACCTCACTgg                                                                                                                                              >  2:756479/1‑200 (MQ=255)
                                   aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTAACAATCGCAGTCAGACACCAGCACTgta                                                                                                                                             >  2:745369/1‑198 (MQ=255)
                                          agTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                        <  3:276351/150‑1 (MQ=255)
                                           gTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGt                                                                                                                                                                                       >  3:149923/1‑150 (MQ=255)
                                           gTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTGTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGt                                                                                                                                                                                       <  3:600545/150‑1 (MQ=255)
                                                  ttCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCTTTACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaaaa                                                                                                                              >  2:346558/1‑200 (MQ=255)
                                                   tCAACGGCGTAGCTACGCGTTTAAAGGATTAATTATGAAGTTGAATTTTAAGGGATTTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                        <  4:418584/141‑1 (MQ=255)
                                                     aaCGGCTTAGCTAGTCGTTTAAAGGATTAATCATGAAGTTGAATTTGAAGGGATTTTTTAAGGCTGGCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATTCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  3:572268/150‑1 (MQ=255)
                                                     aaCGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  3:69893/150‑1 (MQ=255)
                                                     aaCGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  3:69901/150‑1 (MQ=255)
                                                     aaCGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  3:584220/150‑1 (MQ=255)
                                                     aaCGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGAGTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAt                                                                                                                           <  2:132786/200‑1 (MQ=255)
                                                       cGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAg                                                                                                                                                                           >  3:307015/1‑150 (MQ=255)
                                                         gCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCaaa                                                                                                                                                                                           >  1:62629/1‑132 (MQ=255)
                                                         gCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTc                                                                                                                                                                                  <  4:400887/141‑1 (MQ=255)
                                                             aGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATgg                                                                                                                   <  2:396986/200‑1 (MQ=255)
                                                             aGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGa                                                                                                                                                                              >  4:109326/1‑141 (MQ=255)
                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaaaa                                                                                                            >  2:525555/1‑200 (MQ=255)
                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATc                                                                                                                                                                       >  4:113771/1‑141 (MQ=255)
                                                                    cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGt                                                                                                                                                              >  3:581474/1‑150 (MQ=255)
                                                                     gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCATTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATTGAACAAAAg                                                                                                           >  2:759007/1‑200 (MQ=255)
                                                                       ttAAAGGATTAATCATGAAGTTTAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAGACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  1:417487/132‑1 (MQ=255)
                                                                       ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTGTTTAAGGCTGCCTGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  1:218209/132‑1 (MQ=255)
                                                                       gtAAAGGATAAATAATTAAGTTAAATTTTAAGGGATTTTTTAAGGCTGACGGTTGATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACTCCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                                         <  2:609905/199‑1 (MQ=255)
                                                                            ggATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             <  2:213292/127‑1 (MQ=255)
                                                                            ggATTAATCATCAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                                             >  3:213291/1‑127 (MQ=255)
                                                                            ggATTAATCAGGAAGTTGAATTTTAAGGGAGTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAt                                                                                                                                                                        <  1:59097/132‑1 (MQ=255)
                                                                             gATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGATATGtt                                                                                                   >  2:700275/1‑200 (MQ=255)
                                                                             gATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTTACAAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                                                                                   >  2:758675/1‑200 (MQ=255)
                                                                                 aaTCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAGAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTg                                                                                               >  2:685892/1‑200 (MQ=255)
                                                                                  aTCATGAAGTTGAATTTTAAGGGATTTTTTACGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTATAACGAAAAGGGTAGTACACGAAAGTATCAATCGCAGTAAGACACCAGCACTGGGCTCACGCAGGAAAAAGAAAGTAATGGACCAAAAGGCTCTGCTTTCGt                                                                                              >  2:726820/1‑200 (MQ=255)
                                                                                   tCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCTTTACAGAAAAGGGAAGTTCAGGAAAGTATCAATGGCAGCCAGACACCATCACTGGGCTATCTCAGGCAAAAGATAGTAATGGAACAAAAg                                                                                                           >  2:695323/1‑186 (MQ=255)
                                                                                          gTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATAAATCGCAGTCAGACACAATCACTGGGCTCTCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAAt                                                                                      >  2:73565/1‑200 (MQ=255)
                                                                                          gTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTCCTTTACAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATccc                                                                                                                                                 >  4:337946/1‑139 (MQ=255)
                                                                                           ttGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGac                                                                                                                                                         <  1:273463/132‑1 (MQ=255)
                                                                                             gAATTGTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGacac                                                                                                                                                       <  1:248296/132‑1 (MQ=255)
                                                                                               aTTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATACCGCAAAGGGTAGTTc                                                                                                                                                                                  >  4:333041/1‑103 (MQ=255)
                                                                                               aTTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTCCTTTAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGc                                                                                                                                            >  4:95426/1‑141 (MQ=255)
                                                                                               aTTTTAAGAGGTGTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTc                                                                                                                                                                                  <  1:333035/103‑1 (MQ=255)
                                                                                                       ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAt                                                                                                                           >  3:208311/1‑150 (MQ=255)
                                                                                                       ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                                         <  2:395877/200‑1 (MQ=255)
                                                                                                           gtgttAAGGCTGCCGGGTAAATCACACGGTCTCTGATGCTTTGAGTCTGTAGGTCGTAGGCACAGGTTTCCCATACCGCAAAGGGTAGTTCGGGAAAGTATCAATCGGAGTCAGACACCATAACTGGGCtgt                                                                                                                                         <  1:124193/129‑3 (MQ=255)
                                                                                                            ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATc                                                                                                                                        >  1:293465/1‑132 (MQ=255)
                                                                                                            ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTaa                                                                                                                      >  3:371391/1‑150 (MQ=255)
                                                                                                            ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCa                                                                    >  2:196077/1‑200 (MQ=255)
                                                                                                             tttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                   <  2:149923/200‑1 (MQ=255)
                                                                                                                   gCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGt                                                                                                                        <  4:333614/141‑1 (MQ=255)
                                                                                                                       ccGGTTTATACCCACTGGCGTTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAGACCGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCAGGCAATAGATAGTAATGGCACAAAAg                                                                                                           <  3:695283/150‑1 (MQ=255)
                                                                                                                          gTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGt                                                                                                                        >  3:475703/1‑134 (MQ=255)
                                                                                                                          gTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAATTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGt                                                                                                                        <  2:475730/134‑1 (MQ=255)
                                                                                                                                         cgcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAg                                                                                                           <  1:97429/132‑1 (MQ=255)
                                                                                                                                         cgcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGttt                                                                                                  <  4:270994/141‑1 (MQ=255)
                                                                                                                                                           gTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTaa                                                                                                                      >  2:735980/1‑103 (MQ=255)
                                                                                                                                                           gTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTaa                                                                                                                      <  3:735934/103‑1 (MQ=255)
                                                                                                                                                               ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGCCACCATCACTGGGATATCGCAGGCAACAGATAGTAATAGAACAAAAGGCTCTGTATGTGTAGGGGAATCGCTGGAttctctcaccac                                                                   >  3:452199/1‑140 (MQ=255)
                                                                                                                                                               ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                   >  3:317161/1‑150 (MQ=255)
                                                                                                                                                               ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCGATGCTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                   >  3:324284/1‑150 (MQ=255)
                                                                                                                                                               ctcGTATGATCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGGTTACCTTATCAc                                                                   <  3:520288/150‑1 (MQ=255)
                                                                                                                                                                 cGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTg                                                                 >  3:199306/1‑150 (MQ=255)
                                                                                                                                                                  gTATGCTCTGGTTTCCAATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                                         <  4:142667/141‑1 (MQ=255)
                                                                                                                                                                        tctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                   <  4:334546/141‑1 (MQ=255)
                                                                                                                                                                         ctGGTTTCCCATACCGCAAAGGATAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGTGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAatat       <  2:652553/200‑1 (MQ=255)
                                                                                                                                                                           ggTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTc     >  2:338429/1‑200 (MQ=255)
                                                                                                                                                                             tttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCGATAGCAGTCAGACACCATCACTGGGCTAGCGCAGGCAAAAGATAGTAATGGAACAAAATGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATg                                                     <  3:346537/150‑1 (MQ=255)
                                                                                                                                                                              ttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                    >  3:386348/1‑150 (MQ=255)
                                                                                                                                                                              ttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAgg  <  2:208312/200‑1 (MQ=255)
                                                                                                                                                                              ttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTAGCACTAATgg                                                             >  4:202277/1‑141 (MQ=255)
                                                                                                                                                                                 ccATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                                   <  1:47499/132‑1 (MQ=255)
                                                                                                                                                                                 ccATACCGCAAAGGGTAGTGCAGGAAAGTATCAATCGCAGTCAGACACCATCACTTGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACAt                                                 <  3:213276/150‑1 (MQ=255)
                                                                                                                                                                                   aTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTg                                                                 >  1:254323/1‑132 (MQ=255)
                                                                                                                                                                                   |                                                                                                                                                                                                  
ATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGG  >  NC_000913/3875353‑3875726

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.