Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,411,645 9722 bp→A yhdJ[yhdY] yhdJ, yhdU, envR, acrE, acrF, yhdV, yhdX, [yhdY]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3411645 3421366 9722 44 [0] [0] 43 yhdJ–[yhdY] yhdJ,yhdU,envR,acrE,acrF,yhdV,yhdX,[yhdY]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 34116440 (0.000)43 (0.730)
+A
39/492 0.3 100% intergenic (+77/‑9) fis/yhdJ DNA‑binding transcriptional dual regulator Fis/DNA adenine methyltransferase
?NC_000913 3421367 = 0 (0.000)coding (43/1104 nt) yhdY putative ABC transporter membrane subunit YhdY

ACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAA  >  NC_000913/3421128‑3421612
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
aggttaatgacatctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaagccacccgtggtaaccagacccgtgctgcgctgatgatgggcatgaaccgtggtacgttgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaaggcttgattaaaaaggcgctactcggcattgggaagcgtcttttttataggtgtcacaaaggaCATTACTTAAGc                                                                                                                                                                                                                                            <  2:42607‑M2/12‑1 (MQ=255)
    gagtgacctatatgagctggtactggctgaaggagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattaaggttagctaaatgcttgtttaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCcgcg                                                                                                                                                                                                                                        >  2:351771‑M2/236‑251 (MQ=255)
             ctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagaaccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccgtttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGgt                                                                                                                                                                                                                                 <  2:146404‑M2/23‑1 (MQ=255)
                               tgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGt                                                                                                                                                                                                             >  1:229619‑M2/209‑251 (MQ=255)
                                     agaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTcc                                                                                                                                                                                                        <  2:411659‑M2/48‑1 (MQ=255)
                                         cagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTg                                                                                                                                                                                                   >  1:340167‑M2/199‑251 (MQ=255)
                                                        atggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgct                                                                                                                                                                                    >  1:255746‑M2/184‑251 (MQ=255)
                                                                     acacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATATCCTGCTGACTATTGGCTg                                                                                                                                                                        <  2:171227‑M2/80‑1 (MQ=255)
                                                                         ccgtggtaaccagacacgtgatgcgctgatgatgggcatcaaccgtggtacgatgtgtaaaaaattgaaaaaatacggcatgaactaattcaggttagcaaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTTTGGGTGCGAAAAAATCTGGTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTgcatct                                                                                                                                                                   >  2:158217‑M2/167‑249 (MQ=255)
                                                                          cgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATg                                                                                                                                                                  >  2:76904‑M2/166‑251 (MQ=255)
                                                                                aaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGa                                                                                                                                                             <  1:338054‑M2/91‑1 (MQ=255)
                                                                                        ccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAAc                                                                                                                                                    >  1:410954‑M2/152‑251 (MQ=255)
                                                                                            gctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGa                                                                                                                                                 <  2:421105‑M2/103‑1 (MQ=255)
                                                                                               gcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATcc                                                                                                                                              <  1:375736‑M2/106‑1 (MQ=255)
                                                                                                        atgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCt                                                                                                                                    >  1:221914‑M2/136‑251 (MQ=255)
                                                                                                            gcatcaaccgtgatacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttaactaaatgcttgattaaaaaggcgctactcggcatggggaaacgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTGGATGTGGGAACTGATCCCACCGTTGCTGAAc                                                                                                                                <  1:223588‑M2/120‑1 (MQ=255)
                                                                                                              atcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaatgccggctagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCTCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAGCTg                                                                                                                              <  1:197847‑M2/122‑1 (MQ=255)
                                                                                                                        gtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTg                                                                                                                    >  1:255736‑M2/120‑251 (MQ=255)
                                                                                                                         tacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgattgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGACTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGc                                                                                                                   >  2:248403‑M2/119‑251 (MQ=255)
                                                                                                                             ctgcgtacaaaattgaaaaaatacggcatggactatttcaggttagctaaatgattgattaaaaaggcgctactcggcatggggaagcgccttttttatcggtgtcacaaaggaAATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTATCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCCTATGGTTGATGTGGGAACTGGTCCCACCGTTGCTGAACTGGGCATTTTTGCaggg                                                                                                               >  2:143582‑M2/115‑250 (MQ=255)
                                                                                                                                 gtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAAc                                                                                                           >  2:383203‑M2/111‑251 (MQ=255)
                                                                                                                                  taaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACt                                                                                                          >  1:186466‑M2/110‑251 (MQ=255)
                                                                                                                                     aaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTggg                                                                                                       >  1:93108‑M2/107‑251 (MQ=255)
                                                                                                                                           gaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCt                                                                                                 >  1:440576‑M2/101‑251 (MQ=255)
                                                                                                                                           gaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCt                                                                                                 >  1:325795‑M2/101‑251 (MQ=255)
                                                                                                                                                   acggcatgaaataattcaggttagctaaatgcttgattaaaaaggcgctactcggaatggggaagagccttttttataggtgtcacacaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGCTCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGt                                                                                         >  2:140351‑M2/93‑251 (MQ=255)
                                                                                                                                                        atgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCg                                                                                     >  2:191737‑M2/88‑250 (MQ=255)
                                                                                                                                                          gaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACt                                                                                  <  2:282796‑M2/166‑1 (MQ=255)
                                                                                                                                                           aactaagtcaggttagctaaatgcttgattaaaaaggcgctactcggcatggtgaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAAGCTAACTGGGTTGGCTCAACGCGTGCCGACTg                                                                                 <  2:45717‑M2/167‑1 (MQ=255)
                                                                                                                                                                      gttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttagaggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCg                                                                      <  2:245625‑M2/178‑1 (MQ=255)
                                                                                                                                                                        tagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGt                                                                    >  1:401070‑M2/72‑251 (MQ=255)
                                                                                                                                                                        tacctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGAGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCATGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCGGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCACCGCGTGCCGAATGCACAAAAGCCGGt                                                                    >  2:285330‑M2/72‑251 (MQ=255)
                                                                                                                                                                         agctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGt                                                                    <  2:287433‑M2/180‑1 (MQ=255)
                                                                                                                                                                              aatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGt                                                              >  1:381732‑M2/66‑251 (MQ=255)
                                                                                                                                                                                                                 aagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGAGCGATGGTGTGGGTGCGAAAAAAGCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtat                           <  2:272582‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                                   gcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgt                         <  2:325795‑M2/223‑1 (MQ=255)
                                                                                                                                                                                                                    cgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgta                        <  1:158063‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                      ccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgtatg                      <  1:159019‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                      ccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgtatg                      <  2:27702‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                          ttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCt                  >  2:408942‑M2/22‑251 (MQ=255)
                                                                                                                                                                                                                               taggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAAAGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTAcc             >  1:213441‑M2/17‑251 (MQ=255)
                                                                                                                                                                                                                                       acaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGAc     <  1:76904‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                                          aaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCaa  <  2:93108‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
ACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAA  >  NC_000913/3421128‑3421612

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: