Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,907 G→A 9.3% intergenic (+16/+23) yfeN → / ← yfeR putative outer membrane protein/transcriptional regulator of yefH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,9070GA 9.3% 316.4 / 17.2 107intergenic (+16/+23)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base G (48/49);  new base A (4/6);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 109 total aligned reads displayed.

GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCA  >  NC_000913/2525761‑2526054
                                                                                                                                                   |                                                                                                                                                   
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGc                                                                                                                                                    >  4:423678/1‑148 (MQ=255)
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCAATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCt                                                                                                                                                   >  2:879615/1‑149 (MQ=255)
  gcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTcc                                                                                                                                                 >  1:276733/1‑149 (MQ=255)
      ttGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTTCTGCCGGGGCTCCCGGc                                                                                                                                             <  2:776622/149‑1 (MQ=255)
         cACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTCCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGt                                                                                                                                          <  2:852461/149‑1 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGGTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTa                                                                                                                                         <  2:325522/149‑1 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTa                                                                                                                                         >  2:61849/1‑149 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTa                                                                                                                                         <  2:922468/149‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTcc                                                                                                                                                 >  3:380701/1‑140 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTcc                                                                                                                                                 <  4:380701/140‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTa                                                                                                                                         <  4:81374/148‑1 (MQ=255)
            tgtgGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATg                                                                                                                                       >  1:684334/1‑149 (MQ=255)
             gtggtgGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATg                                                                                                                                       >  4:1032737/1‑148 (MQ=255)
              tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATACTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                     <  1:922938/149‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                     <  4:856946/148‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                     <  4:128975/148‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAt                                                                                                                                    >  1:966240/1‑149 (MQ=255)
                gtggATTGCAACATCACACATTACTACTGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTGCTGCCGGGGCTCCCGGCAGTATGAAt                                                                                                                                    <  4:683075/148‑1 (MQ=255)
                gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAtt                                                                                                                                   <  1:884819/149‑1 (MQ=255)
                 tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTc                                                                                                                                  <  1:946694/149‑1 (MQ=255)
                    aaTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAc                                                                                                                                <  3:345581/148‑1 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCActct                                                                                                                             <  1:1067495/149‑1 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCActct                                                                                                                             <  1:879615/149‑1 (MQ=255)
                         aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtat                                                                                                                           <  4:444405/148‑1 (MQ=255)
                          aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                          >  4:1042295/1‑148 (MQ=255)
                          aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                          >  3:713853/1‑148 (MQ=255)
                             tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaatca                                                                                                                                     <  3:477989/134‑5 (MQ=255)
                             tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaatca                                                                                                                                     >  4:477989/1‑130 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      <  3:218667/148‑1 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      <  3:176014/148‑1 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      >  4:447929/1‑148 (MQ=255)
                               acacgTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCTGGGCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                     <  4:1002286/148‑1 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTccc                                                                                                                                                >  3:425171/1‑121 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTccc                                                                                                                                                <  4:425171/121‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTTCTCCATCAAATTTAACTTCTGATTGATTACTGCCTGGGCTCCCGGCAGTATGa                                                                                                                                      <  4:732190/126‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAGCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGa                                                                                                                                      >  3:732190/1‑126 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                <  3:770720/148‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGACGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                >  3:570460/1‑148 (MQ=255)
                                           ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                        >  1:846352/1‑149 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                            <  4:854252/144‑1 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                            >  3:854252/1‑144 (MQ=255)
                                              gATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCt                                                                                                      <  3:479626/148‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                  <  2:24799/149‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAg                                                                                                                                           <  3:472666/108‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAg                                                                                                                                           >  4:472666/1‑108 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                           <  2:170558/104‑1 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                           >  1:170558/1‑104 (MQ=255)
                                                      aTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                              >  3:475752/1‑148 (MQ=255)
                                                       tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGAGAATCTACTCCATGAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                             >  4:1008024/1‑148 (MQ=255)
                                                             tGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                      <  1:57185/149‑1 (MQ=255)
                                                              gCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                      >  3:137468/1‑148 (MQ=255)
                                                                 gATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcg                                                                                   >  3:132400/1‑148 (MQ=255)
                                                                         aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                          <  1:27505/149‑1 (MQ=255)
                                                                          cTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                           <  4:18755/147‑1 (MQ=255)
                                                                              gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtc                                                                      <  3:1042295/148‑1 (MQ=255)
                                                                               cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                     >  4:207036/1‑148 (MQ=255)
                                                                                gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  1:654971/1‑149 (MQ=255)
                                                                                 aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTCACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  3:219119/1‑148 (MQ=255)
                                                                                   gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                 <  4:713853/148‑1 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              >  4:131873/1‑148 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              <  3:1008024/148‑1 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              <  4:435527/148‑1 (MQ=255)
                                                                                            cAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTccc                                                        <  3:1032737/148‑1 (MQ=255)
                                                                                                  gCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAg                                                  >  4:249769/1‑148 (MQ=255)
                                                                                                    gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                               >  1:419965/1‑149 (MQ=255)
                                                                                                    gATAATCTACCCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                               >  1:422781/1‑149 (MQ=255)
                                                                                                      tAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                             <  1:10065/149‑1 (MQ=255)
                                                                                                        aTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAg                                            <  4:475752/148‑1 (MQ=255)
                                                                                                            aCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggc                                        >  4:70714/1‑148 (MQ=255)
                                                                                                              tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGcggcgg                                      <  3:645173/148‑1 (MQ=255)
                                                                                                              tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGcggcgg                                      <  3:76980/148‑1 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTatat                                                                                                                         >  3:715537/1‑62 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTatat                                                                                                                         <  4:715537/62‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                <  2:756657/70‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                >  1:756657/1‑70 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  2:144605/149‑1 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  2:46788/149‑1 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                 >  4:841337/1‑148 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAATATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCGTAAAGCTTCAGCGGCGGTGGaa                                 >  3:648109/1‑148 (MQ=255)
                                                                                                                    aaaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  4:530116/148‑1 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  3:264600/1‑110 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                           <  2:532909/149‑1 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                           <  1:665697/149‑1 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                            >  4:130228/1‑148 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                          <  1:467387/113‑1 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                          >  2:467387/1‑113 (MQ=255)
                                                                                                                              tCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTGTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                      >  3:872756/1‑148 (MQ=255)
                                                                                                                                   ttgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt                >  1:1066909/1‑149 (MQ=255)
                                                                                                                                     gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgc              <  2:1066909/149‑1 (MQ=255)
                                                                                                                                      attACTGCCGGGGCTCCCGGCAGTATGAATTCGCTCTATATCTGATACAACGGACCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGCCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg             <  2:276733/149‑1 (MQ=255)
                                                                                                                                      attACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg             >  1:491192/1‑149 (MQ=255)
                                                                                                                                        tACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc            >  4:158661/1‑148 (MQ=255)
                                                                                                                                        tACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc            >  3:841119/1‑148 (MQ=255)
                                                                                                                                        tACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                        <  4:487452/136‑1 (MQ=255)
                                                                                                                                        tACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                        >  3:487452/1‑136 (MQ=255)
                                                                                                                                         aCTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAc          >  1:856326/1‑149 (MQ=255)
                                                                                                                                           tGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATTAAGCTTCa                                             <  2:1020216/112‑1 (MQ=255)
                                                                                                                                           tGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa        >  1:609952/1‑149 (MQ=255)
                                                                                                                                                gggCTCCCGGCAGTATAAATTCACTCTATATCTGATACAAGGTATCCCATTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGATGGTAAGTGAACGGTTCTTCCTTCGCACCAGCAt    >  4:653650/1‑148 (MQ=255)
                                                                                                                                                   gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAACGTGAACGGTTCTTCCTGCGCACCAGCATCa  <  4:648109/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCA  >  NC_000913/2525761‑2526054

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