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.

GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCA  >  NC_000913/2525761‑2526054
                                                                                                                                                   |                                                                                                                                                   
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑Gc                                                                                                                                                    >  4:423678/1‑148 (MQ=255)
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCAATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCt                                                                                                                                                   >  2:879615/1‑149 (MQ=255)
  gcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTcc                                                                                                                                                 >  1:276733/1‑149 (MQ=255)
      ttGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTTCTGCCGGG‑GCTCCCGGc                                                                                                                                             <  2:776622/149‑1 (MQ=255)
         cACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTCCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGt                                                                                                                                          <  2:852461/149‑1 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGGTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTa                                                                                                                                         <  2:325522/149‑1 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTa                                                                                                                                         >  2:61849/1‑149 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTa                                                                                                                                         <  2:922468/149‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTcc                                                                                                                                                 >  3:380701/1‑140 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTcc                                                                                                                                                 <  4:380701/140‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTa                                                                                                                                         <  4:81374/148‑1 (MQ=255)
            tgtgGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATg                                                                                                                                       >  1:684334/1‑149 (MQ=255)
             gtggtgGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATg                                                                                                                                       >  4:1032737/1‑148 (MQ=255)
              tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATACTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                     <  1:922938/149‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                     <  4:856946/148‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                     <  4:128975/148‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAt                                                                                                                                    >  1:966240/1‑149 (MQ=255)
                gtggATTGCAACATCACACATTACTACTGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTGCTGCCGGG‑GCTCCCGGCAGTATGAAt                                                                                                                                    <  4:683075/148‑1 (MQ=255)
                gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAtt                                                                                                                                   <  1:884819/149‑1 (MQ=255)
                 tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTc                                                                                                                                  <  1:946694/149‑1 (MQ=255)
                    aaTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCAc                                                                                                                                <  3:345581/148‑1 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActct                                                                                                                             <  1:1067495/149‑1 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActct                                                                                                                             <  1:879615/149‑1 (MQ=255)
                         aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCtat                                                                                                                           <  4:444405/148‑1 (MQ=255)
                          aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCtata                                                                                                                          >  4:1042295/1‑148 (MQ=255)
                          aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCtata                                                                                                                          >  3:713853/1‑148 (MQ=255)
                             tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTaatca                                                                                                                                     <  3:477989/134‑5 (MQ=255)
                             tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTaatca                                                                                                                                     >  4:477989/1‑130 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      <  3:218667/148‑1 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      <  3:176014/148‑1 (MQ=255)
                              cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                      >  4:447929/1‑148 (MQ=255)
                               acacgTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCTGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                     <  4:1002286/148‑1 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTccc                                                                                                                                                >  3:425171/1‑121 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTccc                                                                                                                                                <  4:425171/121‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTTCTCCATCAAATTTAACTTCTGATTGATTACTGCCTGG‑GCTCCCGGCAGTATGa                                                                                                                                      <  4:732190/126‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAGCTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGa                                                                                                                                      >  3:732190/1‑126 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                <  3:770720/148‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGACGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                >  3:570460/1‑148 (MQ=255)
                                           ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                        >  1:846352/1‑149 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                            <  4:854252/144‑1 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                            >  3:854252/1‑144 (MQ=255)
                                              gATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCt                                                                                                      <  3:479626/148‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                  <  2:24799/149‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAg                                                                                                                                           <  3:472666/108‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAg                                                                                                                                           >  4:472666/1‑108 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAg                                                                                                                                           <  2:170558/104‑1 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAg                                                                                                                                           >  1:170558/1‑104 (MQ=255)
                                                      aTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                              >  3:475752/1‑148 (MQ=255)
                                                       tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGAGAATCTACTCCATGAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                             >  4:1008024/1‑148 (MQ=255)
                                                             tGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                      <  1:57185/149‑1 (MQ=255)
                                                              gCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                      >  3:137468/1‑148 (MQ=255)
                                                                 gATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcg                                                                                   >  3:132400/1‑148 (MQ=255)
                                                                         aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                          <  1:27505/149‑1 (MQ=255)
                                                                          cTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                           <  4:18755/147‑1 (MQ=255)
                                                                              gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtc                                                                      <  3:1042295/148‑1 (MQ=255)
                                                                               cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                     >  4:207036/1‑148 (MQ=255)
                                                                                gAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  1:654971/1‑149 (MQ=255)
                                                                                 aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTCACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  3:219119/1‑148 (MQ=255)
                                                                                   gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                 <  4:713853/148‑1 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              >  4:131873/1‑148 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              <  3:1008024/148‑1 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                              <  4:435527/148‑1 (MQ=255)
                                                                                            cAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTccc                                                        <  3:1032737/148‑1 (MQ=255)
                                                                                                  gCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAg                                                  >  4:249769/1‑148 (MQ=255)
                                                                                                    gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                               >  1:419965/1‑149 (MQ=255)
                                                                                                    gATAATCTACCCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                               >  1:422781/1‑149 (MQ=255)
                                                                                                      tAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                             <  1:10065/149‑1 (MQ=255)
                                                                                                        aTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAg                                            <  4:475752/148‑1 (MQ=255)
                                                                                                            aCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggc                                        >  4:70714/1‑148 (MQ=255)
                                                                                                              tCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGcggcgg                                      <  3:645173/148‑1 (MQ=255)
                                                                                                              tCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGcggcgg                                      <  3:76980/148‑1 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTatat                                                                                                                         >  3:715537/1‑62 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTatat                                                                                                                         <  4:715537/62‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                <  2:756657/70‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACaa                                                                                                                >  1:756657/1‑70 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  2:144605/149‑1 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  2:46788/149‑1 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                 >  4:841337/1‑148 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAATATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCGTAAAGCTTCAGCGGCGGTGGaa                                 >  3:648109/1‑148 (MQ=255)
                                                                                                                    aaaTTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                <  4:530116/148‑1 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                   >  3:264600/1‑110 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                           <  2:532909/149‑1 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                           <  1:665697/149‑1 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                            >  4:130228/1‑148 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                          <  1:467387/113‑1 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                          >  2:467387/1‑113 (MQ=255)
                                                                                                                              tCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTGTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                      >  3:872756/1‑148 (MQ=255)
                                                                                                                                   ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt                >  1:1066909/1‑149 (MQ=255)
                                                                                                                                     gattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgc              <  2:1066909/149‑1 (MQ=255)
                                                                                                                                      attACTGCCGGG‑GCTCCCGGCAGTATGAATTCGCTCTATATCTGATACAACGGACCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGCCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg             <  2:276733/149‑1 (MQ=255)
                                                                                                                                      attACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg             >  1:491192/1‑149 (MQ=255)
                                                                                                                                        tACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc            >  4:158661/1‑148 (MQ=255)
                                                                                                                                        tACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc            >  3:841119/1‑148 (MQ=255)
                                                                                                                                        tACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                        <  4:487452/136‑1 (MQ=255)
                                                                                                                                        tACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                        >  3:487452/1‑136 (MQ=255)
                                                                                                                                         aCTGCCGGGAGC‑CCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAc          >  1:856326/1‑149 (MQ=255)
                                                                                                                                           tGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATTAAGCTTCa                                             <  2:1020216/112‑1 (MQ=255)
                                                                                                                                           tGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa        >  1:609952/1‑149 (MQ=255)
                                                                                                                                                gg‑gCTCCCGGCAGTATAAATTCACTCTATATCTGATACAAGGTATCCCATTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGATGGTAAGTGAACGGTTCTTCCTTCGCACCAGCAt    >  4:653650/1‑148 (MQ=255)
                                                                                                                                                   gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAACGTGAACGGTTCTTCCTGCGCACCAGCATCa  <  4:648109/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCA  >  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: ‑