Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,587,131 G→C 5.2% Q467H (CAG→CAC 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,587,1310GC 5.2% 566.5 / 16.8 155Q467H (CAG→CACnarQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base G (69/78);  new base C (5/3);  total (74/81)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 159 total aligned reads displayed.

TTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGA  >  NC_000913/2586988‑2587274
                                                                                                                                               |                                                                                                                                               
 tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAt                                                                                                                                           >  2:1027504/1‑149 (MQ=255)
   actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAg                                                                                                                                                 <  2:1448769/141‑1 (MQ=255)
   actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAg                                                                                                                                                 >  1:1448769/1‑141 (MQ=255)
   actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAttt                                                                                                                                         >  1:644148/1‑149 (MQ=255)
    ctactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAGACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttg                                                                                                                                        >  3:125548/1‑149 (MQ=255)
     tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                                                                                                                                       >  2:700921/1‑149 (MQ=255)
     tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                                                                                                                                       <  2:1252401/149‑1 (MQ=255)
     tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                                                                                                                                       >  3:1151879/1‑149 (MQ=255)
      actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                                                                                                                                      <  2:303475/149‑1 (MQ=255)
      actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                                                                                                                                      <  1:798577/149‑1 (MQ=255)
       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgttg                                                                                                                                     <  1:1027720/149‑1 (MQ=255)
       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgttg                                                                                                                                     <  2:910699/149‑1 (MQ=255)
       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgttg                                                                                                                                     <  4:203248/149‑1 (MQ=255)
       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGACACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAt                                                                                                                                           <  2:560145/143‑1 (MQ=255)
        tttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                                                                                                                                       >  2:1071814/1‑146 (MQ=255)
        tttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgt                                                                                                                                       <  1:1071814/146‑1 (MQ=255)
         ttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCa                                                                                                                                   <  4:37465/149‑1 (MQ=255)
           ccGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCaaa                                                                                                                                 <  4:1069413/149‑1 (MQ=255)
           ccGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCaaa                                                                                                                                 <  2:89493/149‑1 (MQ=255)
            cGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAat                                                                                                                                <  1:839911/149‑1 (MQ=255)
             gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAg                                                                                                                                                 >  4:190464/1‑131 (MQ=255)
             gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt                                                                                                                               >  3:456829/1‑149 (MQ=255)
             gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt                                                                                                                               >  3:1465263/1‑149 (MQ=255)
              ccTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatta                                                                                                                              >  3:159764/1‑149 (MQ=255)
               cTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCaaa                                                                                                                                 <  2:598420/145‑1 (MQ=255)
                tGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                            <  3:187137/149‑1 (MQ=255)
                 gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                                                                                                                                      <  2:1354737/138‑1 (MQ=255)
                   cGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGAGAAATGCTGGATACGTTACAAAAGCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcgc                                                                                                                         <  4:125548/149‑1 (MQ=255)
                      tGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGaa                                                                                                                      <  3:342073/149‑1 (MQ=255)
                           cagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGt                                                                                                                 <  4:1117429/149‑1 (MQ=255)
                              gCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATATGTTGCaaa                                                                                                                                 <  4:1261408/130‑1 (MQ=255)
                                   tctcCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGt                                                                                                                                               <  1:667709/111‑1 (MQ=255)
                                    ctcCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGc                                                                                                        <  1:244893/149‑1 (MQ=255)
                                      ccccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGa                                                                                                      >  4:1462489/1‑149 (MQ=255)
                                       cccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCACCTGCAttt                                                                                                                                         <  3:396646/113‑1 (MQ=255)
                                       cccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCACCTGCAttt                                                                                                                                         >  4:396646/1‑113 (MQ=255)
                                          tCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAt                                                                                                          >  4:1142618/1‑141 (MQ=255)
                                          tCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  <  1:1053174/149‑1 (MQ=255)
                                          tCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  <  4:1345354/149‑1 (MQ=255)
                                                  gAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattat                                                                                                                             >  3:1126440/1‑114 (MQ=255)
                                                  gAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgcca                                                                                          >  3:31674/1‑149 (MQ=255)
                                                   aGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaa                                                                                         >  4:678703/1‑149 (MQ=255)
                                                     ggAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                            <  4:409747/112‑1 (MQ=255)
                                                     ggAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                            >  3:409747/1‑112 (MQ=255)
                                                      gAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgc                                                                                      >  4:76311/1‑149 (MQ=255)
                                                       aaaTGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgcc                                                                                     >  1:288231/1‑149 (MQ=255)
                                                        aaTGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGATGAATGCGATGAAGCacgccaacgcca                                                                                    <  2:644148/149‑1 (MQ=255)
                                                              ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGaaa                                                                              >  4:1187620/1‑149 (MQ=255)
                                                              ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGaaa                                                                              <  1:1027504/149‑1 (MQ=255)
                                                                aTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                            <  4:1127128/149‑1 (MQ=255)
                                                                aTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                            <  3:1187620/149‑1 (MQ=255)
                                                                aTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                            >  1:1179179/1‑149 (MQ=255)
                                                                  aCGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  >  2:241763/1‑125 (MQ=255)
                                                                  aCGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  <  1:241763/125‑1 (MQ=255)
                                                                    gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGcc                                                                         <  2:976306/148‑1 (MQ=255)
                                                                     ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGa                                                                                                      <  1:1465966/118‑1 (MQ=255)
                                                                     ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGc                                                                                  >  3:407780/1‑138 (MQ=255)
                                                                     ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGc                                                                                  <  4:407780/138‑1 (MQ=255)
                                                                     ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCAGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTc                                                                                                                           >  4:1241774/1‑97 (MQ=255)
                                                                       aCAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCa                                                                     <  4:31674/149‑1 (MQ=255)
                                                                       aCAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCa                                                                     <  4:456829/149‑1 (MQ=255)
                                                                         aaaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGc                                                                                                                                             >  2:552233/1‑75 (MQ=255)
                                                                         aaaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGc                                                                                                                                             <  1:552233/75‑1 (MQ=255)
                                                                         aaaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGa                                                                                                      <  1:594997/114‑1 (MQ=255)
                                                                          aaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGtt                                                                  <  2:1179179/149‑1 (MQ=255)
                                                                           aaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTc                                                                                                                           >  1:50945/1‑91 (MQ=255)
                                                                            aTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGc                                                                >  3:1297044/1‑149 (MQ=255)
                                                                            aTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGc                                                                <  4:1297044/149‑1 (MQ=255)
                                                                             tCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCg                                                               >  2:1077468/1‑149 (MQ=255)
                                                                              cAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCt                                                                                                              <  2:1196133/101‑1 (MQ=255)
                                                                              cAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCt                                                                                                              >  1:1196133/1‑101 (MQ=255)
                                                                                aaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCa                                                            >  3:881099/1‑149 (MQ=255)
                                                                                     gcgcCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcgcg                                                       >  4:517560/1‑149 (MQ=255)
                                                                                       gcCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGt                                                                                                                 <  3:638476/89‑1 (MQ=255)
                                                                                       gcCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGt                                                                                                                 >  4:638476/1‑89 (MQ=255)
                                                                                          aaaCTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAg                                                                                   <  2:995113/116‑1 (MQ=255)
                                                                                          aaaCTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAg                                                                                   >  1:995113/1‑116 (MQ=255)
                                                                                           aaCTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAc                                                 <  3:1462489/149‑1 (MQ=255)
                                                                                            aCTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACg                                                <  3:76311/149‑1 (MQ=255)
                                                                                               gACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCa                                             >  2:820738/1‑149 (MQ=255)
                                                                                                  ccTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAg                                                                                                 >  4:448723/1‑94 (MQ=255)
                                                                                                  ccTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAg                                                                                                 <  3:448723/94‑1 (MQ=255)
                                                                                                   cTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATc                                          >  1:925879/1‑148 (MQ=255)
                                                                                                                  ccAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCg                          <  3:678703/149‑1 (MQ=255)
                                                                                                                   cAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                         <  4:157116/149‑1 (MQ=255)
                                                                                                                    aaCCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTg                        <  3:517560/149‑1 (MQ=255)
                                                                                                                       ccTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATa                     >  1:1432657/1‑149 (MQ=255)
                                                                                                                       ccTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATa                     >  1:425065/1‑149 (MQ=255)
                                                                                                                           gCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACgg                 <  3:125895/149‑1 (MQ=255)
                                                                                                                              cTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                 <  1:780819/130‑1 (MQ=255)
                                                                                                                              cTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                 >  2:780819/1‑130 (MQ=255)
                                                                                                                               tGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAtt             >  2:221934/1‑149 (MQ=255)
                                                                                                                                tttgcaacAAATGCACCTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACgg                 <  1:414275/137‑1 (MQ=255)
                                                                                                                                tttgcaacAAATGCACCTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACgg                 >  2:414275/8‑144 (MQ=255)
                                                                                                                                ggATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTg            <  4:881099/149‑1 (MQ=255)
                                                                                                                                    gcgcAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTAt        <  2:493970/149‑1 (MQ=255)
                                                                                                                                      gcAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCg      <  1:820738/149‑1 (MQ=255)
                                                                                                                                      acAAATGCACCTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  >  3:766669/2‑57 (MQ=38)
                                                                                                                                      acAAATGCACCTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                  <  4:766669/56‑1 (MQ=38)
                                                                                                                                      acAAATGCACCTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCg      >  2:1172328/2‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGA  >  NC_000913/2586988‑2587274

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