Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009133,608,6820T→C14.4% 14.3 ‑12.4 0.098intergenic (+33/+69)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
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 (0/14):  ref base (53/30):  total (53/45)
Fisher's exact test strand distribution p-value = 4.11e-06
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCA  >  NC_000913/3608484‑3608878
                                                                                                                                                                                                            |                                                                                                                                                                                                       
tCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTc                                                                                                                                                                                                       <  4:335282/200‑1 (MQ=255)
   cGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTcccc                                                                                                                                                                                                    >  4:144786/1‑200 (MQ=255)
    gCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCt                                                                                                                                                                                                   <  4:122752/200‑1 (MQ=255)
                 aGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                      <  4:587800/200‑1 (MQ=255)
                     cATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                  >  4:224068/1‑200 (MQ=255)
                     cATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                  >  4:55795/1‑200 (MQ=255)
                        cACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                                                               >  4:1066196/1‑200 (MQ=255)
                          cTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                                             >  4:327023/1‑200 (MQ=255)
                          cTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                                             >  4:906784/1‑200 (MQ=255)
                             ttGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGTGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGGGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  4:457050/200‑1 (MQ=255)
                             ttGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  4:489610/200‑1 (MQ=255)
                             ttGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  4:526953/200‑1 (MQ=255)
                              tGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                                                         <  4:953624/200‑1 (MQ=255)
                                                 aTCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                                                    >  4:555580/1‑200 (MQ=255)
                                                  tCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTc                                                                                                                                                                                                       <  2:709253/150‑1 (MQ=255)
                                                  tCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTc                                                                                                                                                                                                       <  2:292859/150‑1 (MQ=255)
                                                  tCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGCTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGtn                                                                                                                                                                                                        <  2:461757/150‑2 (MQ=255)
                                                      ccgcGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCt                                                                                                                                                                                                   <  2:689266/147‑1 (MQ=255)
                                                           tCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                                                          <  4:93356/200‑1 (MQ=255)
                                                                  gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                                   >  4:250849/1‑200 (MQ=255)
                                                                   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                      >  2:324349/1‑150 (MQ=255)
                                                                   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                    <  4:324371/152‑1 (MQ=255)
                                                                   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                                  >  4:377841/1‑200 (MQ=255)
                                                                   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                                                  >  4:817615/1‑200 (MQ=255)
                                                                     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGcccc                                                                                                                                                                                             >  3:425236/1‑141 (MQ=255)
                                                                     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGcccc                                                                                                                                                                                             >  3:261169/1‑141 (MQ=255)
                                                                     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGcccc                                                                                                                                                                                             >  3:137047/1‑141 (MQ=255)
                                                                     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                    <  2:700668/150‑1 (MQ=255)
                                                                        ggggATGCCGGGGTTGGGGCGGGCGGTGCTGGAGGAAACGGGGGCGACGTTGCTGGTGACAGCGAAGGCTTAAAGTTTTTGGCGCAGGAGATAAGGAAAACCAATGGCAACATTGAGCGCGACCGGCCCCCCCCCCCCTCTGGGGAGGGGGTTGGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             <  4:950252/199‑1 (MQ=255)
                                                                        ggggAGGACGGGGTGGGGGCGGGCAGTGCGGGCAGATACGGGGGGGACGTTGCGGGTGACAGCGAGTGCGTTAAGTTTTTTGCGCAGGAGATAAGGCAAACCGATCGAAAAATTGGGCGGGATCGGCCCCCCCGCCCCTCGGGGGAGGGGGTTAGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             <  4:831844/199‑1 (MQ=255)
                                                                        aGGGAGGAGCGGGTGGGGGCTGGAGGTGGTGGAAGATAGGGGGCGGCGGGTGTTGGGAAAGGCAAAGGCGTTAAGTTTTTTGCGCGGGAGATAAGGAAAACCGTTCGAAATATGGGCCGCGATCGTCCCCCCCGCCCCTCGGGGGAGAGGTTTAGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             <  4:1089081/200‑1 (MQ=255)
                                                                          ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                        <  3:792087/141‑1 (MQ=255)
                                                                            aTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                      <  3:843698/141‑1 (MQ=255)
                                                                             tGACGGGGTTGGGGGGGGCGGTGCTGGCAGATAGGGGGGGGGCGGGGCTGGTAAGGGCAATGGGGTTAAGATGTTGCGGCAGGAGATAAGGCAAACGGACCGCAAAATTGGGCGCGATCGGCCCCCTCCCCCTTCGGGGGAGGGGGTGGGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                                                        <  4:392382/200‑1 (MQ=255)
                                                                               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGGTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  2:644727/150‑1 (MQ=255)
                                                                               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                                                   >  3:579902/1‑141 (MQ=255)
                                                                               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                                                   >  3:286649/1‑141 (MQ=255)
                                                                               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  2:55794/150‑1 (MQ=255)
                                                                               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  2:850307/150‑1 (MQ=255)
                                                                                 cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                                                        <  2:372757/150‑1 (MQ=255)
                                                                                  gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                                                >  3:313565/1‑141 (MQ=255)
                                                                                   ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                        >  1:546221/1‑132 (MQ=255)
                                                                                       gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                    <  1:845300/132‑1 (MQ=255)
                                                                                       gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTGGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                    <  1:786414/132‑1 (MQ=255)
                                                                                        tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  3:555183/141‑1 (MQ=255)
                                                                                                 gTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                                                          <  1:656754/132‑1 (MQ=255)
                                                                                                     gggCAGATAGGGGGCCACGGGTGCGGGGGCAAGCAAATGCTTTAAGTTTGTTGGGAGGGAAAAAGGCAAACCCGTTCGAAACTTGGGGCGCAACCGGCCCCCCCCCCCTTTGGGGGAGGGGTTTGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                <  4:482009/199‑1 (MQ=255)
                                                                                                        cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                                                                  <  1:51029/132‑1 (MQ=255)
                                                                                                        cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                             >  4:344762/1‑200 (MQ=255)
                                                                                                           aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGCTCGGTACCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAGTGGGCATCGAAGCCaa                                                                                                                                            >  2:674489/1‑150 (MQ=255)
                                                                                                              cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                                  >  3:521790/1‑141 (MQ=255)
                                                                                                               ggggaCGAGGGGGCTGGTGACAGCGATGGCGTTAAGATTGTGGGGAAGAAGAGAAGGCAAACCGACCGCAAAAGTGGGCGCGGTGGGCCCCTTCGCCCTTCTGGGGGGGGGGTGAGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                                      <  4:404825/200‑1 (MQ=255)
                                                                                                                    gggcgggTGCTGGGAACAGGAATGGGGTTAAGGTTGTTGCGCAGGAGAAAAGGCAAACCGATGGCAACATTGGGCCCGATGGGCCCCCCGGCCCCTCTGGGGGGGGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                                                    <  2:321641/145‑1 (MQ=255)
                                                                                                                       cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                              >  4:46323/1‑200 (MQ=255)
                                                                                                                          gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                                                >  1:793216/1‑132 (MQ=255)
                                                                                                                          gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                                                >  1:325833/1‑132 (MQ=255)
                                                                                                                          gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                                       >  3:548389/1‑141 (MQ=255)
                                                                                                                          gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                                                            >  4:512448/1‑200 (MQ=255)
                                                                                                                               ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                                                         >  2:926397/1‑150 (MQ=255)
                                                                                                                                    cAGCGAATGAGTTAAGATTGTGGCGCAGGAAAAAAGGAAAACAGATGGAAACATGGAGCGCAATGGGTCCCCGCGCCCCTGTGGGGAGAGGGTTGGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             <  3:325847/141‑1 (MQ=255)
                                                                                                                                      gCGAATGCGTTAGGGTTGTTGCGCAGGAGATAAGGCAAACCGATGGCAACATTGGGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                                <  4:672162/200‑1 (MQ=255)
                                                                                                                                             gCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAACAACCGCCTTTACGAAACaat                                                                                                            >  2:495712/1‑149 (MQ=255)
                                                                                                                                              cGTTACGATCGTTGCGCAGGAGATAAGGCAAACCTATCGCAACATCGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGTGGAAAAGGCGGCATCGAAGCCAACCTGCCCCTATCAAc                                                                                                                             >  1:267217/1‑132 (MQ=37)
                                                                                                                                              cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGTCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             >  1:407044/1‑132 (MQ=255)
                                                                                                                                              cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                                    >  3:96811/1‑141 (MQ=255)
                                                                                                                                              cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACACTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGTCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                             >  1:163897/1‑132 (MQ=255)
                                                                                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCGTTACGa                                                                                                                  >  3:42725/1‑141 (MQ=255)
                                                                                                                                                   aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATc                                                                                                               >  3:573350/1‑141 (MQ=255)
                                                                                                                                                   aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                                                    >  4:95940/1‑200 (MQ=255)
                                                                                                                                                    gATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                                   >  4:124319/1‑200 (MQ=255)
                                                                                                                                                         ttgCGCAGGAGAAAGGCCAAACGGATCGAACAATTAGGCGCGATGGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                <  2:817548/150‑1 (MQ=255)
                                                                                                                                                          ttcgcAGGAGAAAGGGAAAACCGCCCGAAACATTGAGCGCGGTCGTCCCCCTCGCCCCTCGGGGGAGAGGGTAGGGGTAAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                             <  4:607774/198‑1 (MQ=255)
                                                                                                                                                                 ggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                                                     >  3:654366/1‑138 (MQ=255)
                                                                                                                                                                   agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                                      >  4:940173/1‑200 (MQ=255)
                                                                                                                                                                   agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                         >  1:725390/1‑132 (MQ=255)
                                                                                                                                                                   agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                >  3:35180/1‑141 (MQ=255)
                                                                                                                                                                   agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                                     >  4:520610/1‑200 (MQ=255)
                                                                                                                                                                   agaTAAGCCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCTCCTATGGGGAGATGGCTAGGGAGAGGGGAAAAGGCAGCATCGAAGCAAATCATCCCCTATCAACCGCCTTACCGAATCAAATAACTTTAAGGCAGCCCAtac                                                                                       >  2:436474/1‑148 (MQ=255)
                                                                                                                                                                    gaaaaGGAAAACGGATCGAAACATTGAGCGCGTTCGGCCCCCCCCCCCCTCGGGGGAGAGGGTTGGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                                      <  2:555516/147‑1 (MQ=255)
                                                                                                                                                                    gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                                      >  2:392119/1‑150 (MQ=255)
                                                                                                                                                                    aaaaaGGCAAAACGGACCGAAACTTTGGGCGCATTGGTCCCCCTCCCCCCTCGGGGGGGAGGGTTAGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                <  3:621062/138‑1 (MQ=255)
                                                                                                                                                                     aTAGGGCAAACGGATCGAAACATTGGGCGAGTTCGGCCCCCCCGCCCCTGGGGGGGGGGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACg                                                                                                       <  1:573329/132‑1 (MQ=255)
                                                                                                                                                                        gagggaaagccGACCGCAACATGGGGCACGATCAACCCCCCCGCCCATCTGGGGAGAGGGTTGGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                                                    <  2:224060/141‑1 (MQ=255)
                                                                                                                                                                        acgcaaaaCCGATGGAAACATTGCGCGCGCTCGGCCCCCCCCCCCTTCTGGGGAGAGGGTTAGGTGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                                                     <  1:654344/127‑1 (MQ=255)
                                                                                                                                                                        aGGCAAACCATTCGAAACTTGAGCCGAGATCGCCCCCCCCCCCCCTTGGGGGAGAGGTTAGGGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                                                           <  3:601269/141‑1 (MQ=255)
                                                                                                                                                                         ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAACTAACGATaa                                                                                                   >  1:32303/1‑132 (MQ=255)
                                                                                                                                                                         ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                     >  4:18558/1‑178 (MQ=255)
                                                                                                                                                                             aaaCCGACCGCAACATTGGGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                <  1:286637/132‑1 (MQ=255)
                                                                                                                                                                              aaacggagggaaaCAGTGGGCGCGATGGGGCCCCTCGCCCCTAGGGGGAGAGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                                                <  1:459433/121‑1 (MQ=255)
                                                                                                                                                                                accgaggacAACATGGAGCGCGATGGGTCCCCGCGCCCCTCGGGGGGGAGGTTGAGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                                                            <  2:325188/142‑1 (MQ=255)
                                                                                                                                                                                     tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                                                                >  3:601258/1‑93 (MQ=255)
                                                                                                                                                                                     tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                          >  2:680434/1‑99 (MQ=255)
                                                                                                                                                                                     tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                                       >  2:1067316/1‑150 (MQ=255)
                                                                                                                                                                                     ncGCAACGTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCTAATCAGCCACTATCAACCGCCTTTACGAATCAAATAACGATCGGGCAGTCCATCCGTCTCCATAGTAAcaa                                                                       >  2:72/2‑148 (MQ=37)
                                                                                                                                                                                       gCAACATTGAGCGCGATCGGTCCCCTCGCGCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACGGCGGCGTCGAAGACAATCAGCCCCTACCAACCGCCCTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGAAACCAGt                                                                     >  2:941687/1‑150 (MQ=38)
                                                                                                                                                                                                 gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                                    >  3:442865/1‑141 (MQ=255)
                                                                                                                                                                                                        tCGGTCCCCCCCCCCCTCTGGGGGGAGGGTTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                                                     <  2:18557/150‑1 (MQ=255)
                                                                                                                                                                                                         cGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCa  >  4:197847/1‑200 (MQ=255)
                                                                                                                                                                                                           gTCCCCTCGCCCCTCTGGGGTGAGGGCTAGGGTGTGGGGACAAGGCGGCATCTAAGCCAATCACATCCTATCAACAGCCTTTACGAATCAAATAACGCTAAGGCCGTCAATCCGTCTCTCTAGCATCCAggt                                                                    >  1:714037/1‑130 (MQ=25)
                                                                                                                                                                                                            |                                                                                                                                                                                                       
TCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCA  >  NC_000913/3608484‑3608878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 19 ≤ 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.