Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,586,946 T→A 5.5% S406T (TCC→ACC)  narQ → sensory histidine kinase in two‑component regulatory system with NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,586,9460TA 5.5% 291.6 / 2.6 87S406T (TCC→ACC) narQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base A (2/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGATGGAAGAACGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCG  >  NC_000913/2586801‑2587090
                                                                                                                                                 |                                                                                                                                                
gttgATGGAAGAACGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTc                                                                                                                                                 <  2:266341/147‑1 (MQ=255)
gttgATGGAAGAACGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTc                                                                                                                                                 <  7:16127/147‑1 (MQ=255)
 ttgATGGAAGAACGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTcc                                                                                                                                                <  6:153943/147‑1 (MQ=255)
       gaagaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGt                                                                                                                                         <  3:391598/148‑1 (MQ=255)
        aagaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGtt                                                                                                                                        >  8:400212/1‑148 (MQ=255)
         agaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTcc                                                                                                                                                <  8:402803/139‑1 (MQ=255)
         agaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTcc                                                                                                                                                >  7:402803/1‑139 (MQ=255)
         agaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGtt                                                                                                                                        >  7:415584/1‑147 (MQ=255)
         agaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGtt                                                                                                                                        >  1:69555/1‑147 (MQ=255)
          gaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGa                                                                                                                                      <  8:415584/148‑1 (MQ=255)
          gaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGa                                                                                                                                      <  2:69555/148‑1 (MQ=255)
          gaaCGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGa                                                                                                                                      <  1:9297/148‑1 (MQ=255)
             cGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAatg                                                                                                                                   <  7:400212/148‑1 (MQ=255)
             cGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAat                                                                                                                                    >  4:183336/1‑147 (MQ=255)
               tGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAatgat                                                                                                                                 <  6:167901/148‑1 (MQ=255)
                 cGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAatgatg                                                                                                                                <  7:416001/147‑1 (MQ=255)
                    ccATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTa                                                                                                                            <  5:77445/148‑1 (MQ=255)
                      aTCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTa                                                                                                                            <  6:172072/146‑1 (MQ=255)
                      aTCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTa                                                                                                                            >  5:172072/1‑146 (MQ=255)
                       tCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCg                                                                                                                         <  3:45713/148‑1 (MQ=255)
                          cccGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCa                                                                                                                      <  4:400733/148‑1 (MQ=255)
                            cgcgAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAg                                                                                                                     >  2:246355/1‑147 (MQ=255)
                            cgcgAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAg                                                                                                                     <  1:246355/147‑1 (MQ=255)
                               gAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAc                                                                                                                 >  3:61619/1‑148 (MQ=255)
                                     cACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGc                                                                                                           >  7:336911/1‑148 (MQ=255)
                                      aCGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGa                                                                                                             <  2:152237/145‑1 (MQ=255)
                                      aCGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGa                                                                                                             >  1:152237/1‑145 (MQ=255)
                                      aCGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCt                                                                                                          >  6:67887/1‑148 (MQ=255)
                                       cGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTg                                                                                                         >  2:202724/1‑148 (MQ=255)
                                       cGACTCGCTGGCTCAGGTACTTTCATACTTACGTATCCAGTTGACGTTACTGACGCGGTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTg                                                                                                         <  6:367178/148‑1 (MQ=255)
                                          cTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTg                                                                                                      <  6:220955/148‑1 (MQ=255)
                                              cTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGc                                                                                                                       <  7:346789/127‑1 (MQ=255)
                                              cTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGc                                                                                                                       >  8:346789/1‑127 (MQ=255)
                                                     aGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCg                                                                                           >  3:91372/1‑148 (MQ=255)
                                                       gTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGcc                                                                                         >  8:29561/1‑148 (MQ=255)
                                                        tACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAcg                                                                                                                <  7:379422/124‑1 (MQ=255)
                                                        tACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAcg                                                                                                                >  8:379422/1‑124 (MQ=255)
                                                             tCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGAc                                                                                     <  7:74017/146‑1 (MQ=255)
                                                                 acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCt                                                                                  <  5:67887/145‑1 (MQ=255)
                                                                 acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCACGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTTACGCTg                                                                                 >  1:264272/1‑146 (MQ=255)
                                                                 actaacGAATCCACTTGACGCTACTCAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTCGACTACTTTCCGCCTGACGCTg                                                                                 <  4:61619/146‑1 (MQ=255)
                                                                  cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCTGGGAAAAATCggccatgatactttgtgcggtggcgttatcttc                                                                                                             <  2:60713‑M2/117‑34 (MQ=255)
                                                                  cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCTGGGAAAAATCggccatgatactttgtgcggtggcgttatcttc                                                                                                             >  1:60713‑M2/1‑84 (MQ=255)
                                                                  cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcag                                                                              >  5:344724/1‑148 (MQ=255)
                                                                    tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCTGGGAAAAATCggccatgatactttgtgcggtggcgttatcttccggtat                                                                                                       >  8:188795‑M2/1‑82 (MQ=255)
                                                                    tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCTGGGAAAAATCggccatgatactttgtgcggtggcgttatcttccggtat                                                                                                       <  7:188795‑M2/121‑40 (MQ=255)
                                                                     acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcagca                                                                            <  3:183336/147‑1 (MQ=255)
                                                                      cGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAgg                                                                          >  4:82910/1‑148 (MQ=255)
                                                                           ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAcgcg                                                                                                              <  3:135984/107‑1 (MQ=255)
                                                                           ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAcgcg                                                                                                              >  4:135984/1‑107 (MQ=255)
                                                                           ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAt                                                                     >  5:384876/1‑148 (MQ=255)
                                                                           ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAt                                                                     >  2:373880/1‑148 (MQ=255)
                                                                           ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAt                                                                     >  7:77241/1‑148 (MQ=255)
                                                                            cAGTGGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtc                                                                    <  1:246140/148‑1 (MQ=255)
                                                                               ttGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcc                                                                 >  5:155191/1‑148 (MQ=255)
                                                                                ttaCGTTACTGAAGAGTTCGATACCGGAAGTTAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcc                                                                 <  6:344724/145‑1 (MQ=255)
                                                                                      tACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAgg                                                                          <  8:256562/132‑1 (MQ=255)
                                                                                      tACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAgg                                                                          >  7:256562/1‑132 (MQ=255)
                                                                                      tACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCg                                                           >  3:383081/1‑147 (MQ=255)
                                                                                      tACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCg                                                           >  8:347152/1‑147 (MQ=255)
                                                                                       aCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtc                                                                    <  4:93998/137‑1 (MQ=255)
                                                                                       aCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtc                                                                    >  3:93998/1‑137 (MQ=255)
                                                                                            aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGt                                                                                                        >  8:148131/1‑96 (MQ=255)
                                                                                            aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGt                                                                                                        <  7:148131/96‑1 (MQ=255)
                                                                                                     tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCTTTGAGGGAAATGCt                                           <  4:383081/148‑1 (MQ=255)
                                                                                                     tACCGGAAGATAACGCCACCGCACAAAGTATCATGGACGATTTTTCCCAGTCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCt                                           <  6:296304/148‑1 (MQ=255)
                                                                                                        cGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGTATCTCCCCTCCGCATTGAGGGAAATGCttg                                         <  2:264272/147‑3 (MQ=255)
                                                                                                         ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACttt                                                                                              <  4:240531/93‑1 (MQ=255)
                                                                                                         ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACttt                                                                                              >  3:240531/1‑93 (MQ=255)
                                                                                                         ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTc                                                                                             >  3:395697/1‑94 (MQ=255)
                                                                                                         ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGa                                        >  5:407665/1‑147 (MQ=255)
                                                                                                           aaGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATAc                                     <  8:336911/148‑1 (MQ=255)
                                                                                                               tAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTa                                 >  4:7674/1‑148 (MQ=255)
                                                                                                               tAACGACACCGCAAAAAGTATCATGGGGGTTTTTACCCAGGCGTTTAATGATGCTTATGGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATgg                                            <  5:126603/137‑2 (MQ=255)
                                                                                                                   gCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaa                             <  4:91372/148‑1 (MQ=255)
                                                                                                                        cGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaa                        <  7:347152/148‑1 (MQ=255)
                                                                                                                         gCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaaa                       <  6:155191/148‑1 (MQ=255)
                                                                                                                         gCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaaa                       >  2:196380/1‑148 (MQ=255)
                                                                                                                           acaAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAAc                      >  3:150959/1‑147 (MQ=255)
                                                                                                                            caAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGAAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaa                             <  6:192772/139‑1 (MQ=255)
                                                                                                                            caAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGAAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaa                             >  5:192772/1‑139 (MQ=255)
                                                                                                                                   tCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCaa             <  3:7674/148‑1 (MQ=255)
                                                                                                                                   tCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCaa             <  3:260217/148‑1 (MQ=255)
                                                                                                                                     aTGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAc           >  8:366177/1‑148 (MQ=255)
                                                                                                                                        gCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTg         <  3:82910/147‑1 (MQ=255)
                                                                                                                                         ccGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaaa                       <  6:47061/132‑1 (MQ=255)
                                                                                                                                         ccGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCaaa                       >  5:47061/1‑132 (MQ=255)
                                                                                                                                            cgttctCCCAGGCGTTGAATGTTGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCTTTGAGGGAAATGCTGGATACGTTACAAATTCAAACCAGCGCCAAACTGACCCt    <  3:5316/143‑1 (MQ=255)
                                                                                                                                               tttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCg  <  1:202724/147‑1 (MQ=255)
                                                                                                                                               tttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCg  <  1:196380/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTGATGGAAGAACGTGCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCG  >  NC_000913/2586801‑2587090

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