Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,908 C→G 5.5% intergenic (+17/+22) 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,9080CG 5.5% 391.0 / 11.1 110intergenic (+17/+22)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (52/52);  new base G (4/2);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 110 total aligned reads displayed.

GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATC  >  NC_000913/2525761‑2526053
                                                                                                                                                    |                                                                                                                                                 
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCt                                                                                                                                                  >  3:415552/1‑149 (MQ=255)
   cAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTccc                                                                                                                                               <  3:545674/149‑1 (MQ=255)
        gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGTGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAg                                                                                                                                          <  2:1343706/149‑1 (MQ=255)
        gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAg                                                                                                                                          <  3:22247/149‑1 (MQ=255)
         cACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGt                                                                                                                                         >  1:1544488/1‑149 (MQ=255)
          aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTa                                                                                                                                        >  1:315696/1‑149 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTaa                                                                                                                                       <  3:620536/149‑2 (MQ=255)
              tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                    <  4:1259578/149‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAt                                                                                                                                   <  3:577138/149‑1 (MQ=255)
                     aTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                             >  1:616361/1‑149 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGC‑CCCGGCAg                                                                                                                                          <  2:655346/135‑1 (MQ=255)
                      tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGC‑CCCGGCAg                                                                                                                                          >  1:655346/1‑135 (MQ=255)
                         aaCATCACATATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCtata                                                                                                                         <  1:956557/149‑1 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAg                                                                                                                                          >  3:239113/1‑126 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAg                                                                                                                                          >  2:1275503/1‑126 (MQ=255)
                               acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAg                                                                                                                                          <  1:1275503/126‑1 (MQ=255)
                                       ctacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                           <  4:415552/149‑1 (MQ=255)
                                        tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGC‑CCCGGCAGt                                                                                                                                         >  1:1163519/1‑118 (MQ=255)
                                        tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGC‑CCCGGCAGt                                                                                                                                         <  2:1163519/118‑1 (MQ=255)
                                           ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                           >  3:1549775/1‑145 (MQ=255)
                                           ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                           <  4:1549775/145‑1 (MQ=255)
                                             ggATGCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCt                                                                                                     >  2:322235/1‑149 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                 <  2:1544488/149‑1 (MQ=255)
                                                   aGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                               >  1:598270/1‑149 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCAc                                                                                                                               <  1:380770/115‑1 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCAc                                                                                                                               >  2:380770/1‑115 (MQ=255)
                                                     tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                             <  1:322235/149‑1 (MQ=255)
                                                      aTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                            >  3:1205566/1‑149 (MQ=255)
                                                                      cAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                               <  2:1660348/130‑1 (MQ=255)
                                                                      cAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                               >  1:1660348/1‑130 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                             >  1:909131/1‑99 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                             <  2:909131/99‑1 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAt                                                                                       <  3:1011921/137‑1 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAt                                                                                       >  4:1011921/1‑137 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGcc                                                                           <  3:755762/149‑1 (MQ=255)
                                                                                 aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                            >  1:702196/1‑122 (MQ=255)
                                                                                 aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                            <  2:702196/122‑1 (MQ=255)
                                                                                  aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACgg                                                                                                            >  4:338591/1‑105 (MQ=255)
                                                                                  aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                <  2:315696/149‑1 (MQ=255)
                                                                                     tttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                     <  3:74702/125‑1 (MQ=255)
                                                                                     tttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                     >  4:74702/1‑125 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                                                            <  4:1435509/149‑1 (MQ=255)
                                                                                            cAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATc                                                                                                                       >  3:256131/1‑84 (MQ=255)
                                                                                            cAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATc                                                                                                                       <  4:256131/84‑1 (MQ=255)
                                                                                             aGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGt                                                          >  3:1151446/1‑144 (MQ=255)
                                                                                             aGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGt                                                          <  4:1151446/144‑1 (MQ=255)
                                                                                               gATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                >  2:242566/1‑136 (MQ=255)
                                                                                               gATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                <  1:242566/136‑1 (MQ=255)
                                                                                                   cGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCta                                                                                                                           >  1:685468/1‑73 (MQ=255)
                                                                                                   cGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCta                                                                                                                           <  2:685468/73‑1 (MQ=255)
                                                                                                    gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                              >  1:139956/1‑149 (MQ=255)
                                                                                                            aCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggcg                                      >  1:1534954/1‑149 (MQ=255)
                                                                                                               ccATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                       <  2:1311816/81‑1 (MQ=255)
                                                                                                               ccATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATccc                                                                                                       >  1:1311816/1‑81 (MQ=255)
                                                                                                                cATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                  >  4:1168738/1‑149 (MQ=255)
                                                                                                                cATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                  <  4:1205566/149‑1 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                 >  1:268994/1‑85 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTccc                                                                                                 <  2:268994/85‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                          <  1:171293/107‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                  >  2:683013/1‑131 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggcg                                      >  4:1639491/1‑143 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggcg                                      <  3:1639491/143‑1 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                >  1:1033968/1‑149 (MQ=255)
                                                                                                                  tCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCAGaaa                                                  <  1:683013/131‑1 (MQ=255)
                                                                                                                   cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                               <  3:1654775/149‑1 (MQ=255)
                                                                                                                      aTTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                              >  3:1079758/1‑67 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                <  2:1471600/144‑1 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                >  1:1471600/1‑144 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                           <  2:598270/149‑1 (MQ=255)
                                                                                                                            cTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                        >  1:1217580/1‑115 (MQ=255)
                                                                                                                            cTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                        <  2:1217580/115‑1 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAt                                                                                                                                   <  4:679301/39‑1 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAt                                                                                                                                   >  3:679301/1‑39 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                            >  1:992746/1‑126 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                            <  2:992746/126‑1 (MQ=255)
                                                                                                                             ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                     >  3:284945/1‑149 (MQ=255)
                                                                                                                              tCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                  <  3:574632/119‑1 (MQ=255)
                                                                                                                              tCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                  >  4:574632/1‑119 (MQ=255)
                                                                                                                                 gattgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGttcttc                 <  2:1534954/149‑1 (MQ=255)
                                                                                                                                  attgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                                        <  1:1044070/93‑1 (MQ=255)
                                                                                                                                  attgattACAGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCCGCGGCGGTGGAAAGTGAACGGttcttca                >  4:1627329/1‑148 (MQ=255)
                                                                                                                                   ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTg                            <  2:226726/136‑1 (MQ=255)
                                                                                                                                   ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTg                            >  1:226726/1‑136 (MQ=255)
                                                                                                                                   ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt               <  2:1033968/149‑1 (MQ=255)
                                                                                                                                   ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt               >  4:494259/1‑149 (MQ=255)
                                                                                                                                    tgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTg              >  2:1650985/1‑149 (MQ=255)
                                                                                                                                      attACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg            >  1:443608/1‑149 (MQ=255)
                                                                                                                                       ttACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                    >  3:827017/1‑92 (MQ=255)
                                                                                                                                       ttACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc           <  4:284945/149‑1 (MQ=255)
                                                                                                                                       atACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                               <  2:489340/128‑1 (MQ=255)
                                                                                                                                        tACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa          >  2:1802471/1‑149 (MQ=255)
                                                                                                                                         aCTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                           <  2:36250/51‑1 (MQ=255)
                                                                                                                                         aCTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGa                                                                                                           >  1:36250/1‑51 (MQ=255)
                                                                                                                                           tGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa       >  3:1344275/1‑149 (MQ=255)
                                                                                                                                           tGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa       >  4:62416/1‑149 (MQ=255)
                                                                                                                                             ccGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc     >  4:1787416/1‑149 (MQ=255)
                                                                                                                                               ggg‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCAt   >  4:288898/1‑149 (MQ=255)
                                                                                                                                                gg‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATc  >  3:1270062/1‑149 (MQ=255)
                                                                                                                                                 g‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa       <  4:1710162/143‑1 (MQ=255)
                                                                                                                                                 g‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa       >  3:1710162/1‑143 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATC  >  NC_000913/2525761‑2526053

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: ‑