Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,306 G→C 5.9% I13M (ATC→ATG yhfZ ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,3060GC 5.9% 416.1 / 14.4 118I13M (ATC→ATGyhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (46/65);  new base C (4/3);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 122 total aligned reads displayed.

TTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTG  >  NC_000913/3512159‑3512451
                                                                                                                                                   |                                                                                                                                                 
ttCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACTGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGAGTCGtg                                                                                                                                                   <  2:436940/148‑1 (MQ=255)
  cAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtgatg                                                                                                                                                >  4:771409/1‑149 (MQ=255)
   aaTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtgatga                                                                                                                                               <  3:72217/149‑1 (MQ=255)
   aaTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtgatga                                                                                                                                               >  4:673055/1‑149 (MQ=255)
     tCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgacg                                                                                                                                             >  3:398015/1‑149 (MQ=255)
          aTCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGc                                                                                                                                         >  1:452681/1‑148 (MQ=255)
           tCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                       <  3:599919/149‑3 (MQ=255)
            cGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                       <  1:900918/148‑3 (MQ=255)
             gCTCCGCTTGATTCCAGTGTTTTCAACGCGGCCTGCGTCAGGCCAACGGCTGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttc                                                                                                                                     <  4:804822/149‑1 (MQ=255)
                ccGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttctt                                                                                                                                   >  2:648534/1‑148 (MQ=255)
                  gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                <  3:325287/149‑1 (MQ=255)
                  gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                 <  2:432528/148‑1 (MQ=255)
                  gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                 <  1:222651/148‑1 (MQ=255)
                  gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                 <  2:403528/148‑1 (MQ=255)
                   cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                               <  3:252909/149‑1 (MQ=255)
                   cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                               <  3:738786/149‑1 (MQ=255)
                   cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                               <  3:661712/149‑1 (MQ=255)
                    ttGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTg                                                                                                                              >  4:509487/1‑149 (MQ=255)
                    ttGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTg                                                                                                                              >  3:577929/1‑149 (MQ=255)
                      gATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGAt                                                                                                                            >  4:872169/1‑149 (MQ=255)
                      gATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGAt                                                                                                                            >  4:47674/1‑149 (MQ=255)
                         tCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                         <  3:319235/149‑1 (MQ=255)
                         tCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                         <  4:459733/149‑1 (MQ=255)
                            aGCGTTTTCAACGCGGCCTGCGTCAGGCTAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAct                                                                                                                       >  2:349653/1‑147 (MQ=255)
                             gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGtt                                                                                                                     <  4:577929/149‑1 (MQ=255)
                              cGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTc                                                                                                                    <  3:775379/149‑1 (MQ=255)
                                ttttCAACGCTGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCg                                                                                                                   >  1:506089/1‑148 (MQ=255)
                                   tCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAGCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGc                                                                                                               <  4:675377/149‑1 (MQ=255)
                                     aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                             >  3:884084/1‑149 (MQ=255)
                                     aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                             >  3:754887/1‑149 (MQ=255)
                                     aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                              >  2:67913/1‑148 (MQ=255)
                                           gCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATaaa                                                                                                       <  4:227522/149‑1 (MQ=255)
                                               gCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtttt                                                                                                   <  3:111431/149‑1 (MQ=255)
                                                cGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtttt                                                                                                   <  2:856613/148‑1 (MQ=255)
                                                 gTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCa                                                                                                 <  3:872169/149‑1 (MQ=255)
                                                 gTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCa                                                                                                 >  3:581502/1‑149 (MQ=255)
                                                    aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                              >  3:433176/1‑149 (MQ=255)
                                                      gCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaaa                                                                                            <  4:866407/149‑1 (MQ=255)
                                                       ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                                                                                           >  3:810692/1‑149 (MQ=255)
                                                       ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                                                                                           >  4:533004/1‑149 (MQ=255)
                                                       ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAAGAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTaaa                                                                                            <  1:482096/148‑1 (MQ=255)
                                                         aaCGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                               <  3:19951/111‑1 (MQ=255)
                                                         aaCGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                               >  4:19951/1‑111 (MQ=255)
                                                           cGGATGAACGGCATTCATTTGCCAGCTCATCTATGTTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                       <  3:44324/149‑1 (MQ=255)
                                                           cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGa                                                                                        <  2:182546/148‑1 (MQ=255)
                                                           cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGa                                                                                        <  1:510162/148‑1 (MQ=255)
                                                             gATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAat                                                                                      <  1:109568/148‑1 (MQ=255)
                                                                gAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                  <  4:398015/149‑1 (MQ=255)
                                                                  aCGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtat                                                                                <  3:509487/149‑1 (MQ=255)
                                                                        ttcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                <  2:357091/95‑1 (MQ=255)
                                                                        ttcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc                                                                          <  3:899508/149‑1 (MQ=255)
                                                                         tcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc                                                                          >  1:111148/1‑148 (MQ=255)
                                                                          cattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa                                                                         >  1:484866/1‑148 (MQ=255)
                                                                          cattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg                                                                        <  3:47674/149‑1 (MQ=255)
                                                                             ttGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa                                                                     <  4:581502/149‑1 (MQ=255)
                                                                                   gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                               <  4:785868/149‑1 (MQ=255)
                                                                                   gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGa                                                                <  2:52594/148‑1 (MQ=255)
                                                                                   gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAATAAATTGAATATTTATATCTTCAGGAATTTGaa                                                               <  4:449355/149‑1 (MQ=255)
                                                                                     tCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAg                                                              <  2:506089/148‑1 (MQ=255)
                                                                                     tCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGa                                                             >  4:249282/1‑149 (MQ=255)
                                                                                      cATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGa                                                             >  1:494622/1‑148 (MQ=255)
                                                                                          tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                        >  3:856783/1‑149 (MQ=255)
                                                                                          tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                        >  4:886293/1‑149 (MQ=255)
                                                                                           aTGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGa                                                                >  2:783108/1‑140 (MQ=255)
                                                                                           aTGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGa                                                                <  1:783108/140‑1 (MQ=255)
                                                                                            tGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                       >  2:511424/1‑148 (MQ=255)
                                                                                              gTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCa                                                    >  3:274283/1‑149 (MQ=255)
                                                                                               ttttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                  >  4:260210/1‑118 (MQ=255)
                                                                                                     aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATCAACGTTCGACGCACAATAAATTATCAGTAAAGTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCga                                             >  4:775619/1‑149 (MQ=255)
                                                                                                     aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCga                                             >  3:497389/1‑149 (MQ=255)
                                                                                                     aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCga                                             <  3:40628/149‑1 (MQ=255)
                                                                                                        gATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCgaga                                           <  1:511424/148‑1 (MQ=255)
                                                                                                          ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                         <  1:880847/148‑1 (MQ=255)
                                                                                                          ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                         >  2:101560/1‑148 (MQ=255)
                                                                                                          ttACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                         >  2:45809/1‑148 (MQ=255)
                                                                                                           tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                       <  3:868853/149‑1 (MQ=255)
                                                                                                             ccGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                       >  4:709522/1‑131 (MQ=255)
                                                                                                             ccGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                       <  3:709522/131‑1 (MQ=255)
                                                                                                              cGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGa                                     >  1:909556/1‑148 (MQ=255)
                                                                                                               gCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGaaa                                   <  4:392152/149‑1 (MQ=255)
                                                                                                                cACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                       >  2:810450/1‑128 (MQ=255)
                                                                                                                cACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                       <  1:810450/128‑1 (MQ=255)
                                                                                                                     ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCg                              <  1:375088/148‑1 (MQ=255)
                                                                                                                           caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACt                       >  4:178948/1‑149 (MQ=255)
                                                                                                                             acaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                      >  2:252540/1‑148 (MQ=255)
                                                                                                                              cacacaaCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTc                     >  2:19604/6‑148 (MQ=255)
                                                                                                                                 aTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATaaa                                                                                                       >  4:478397/1‑63 (MQ=255)
                                                                                                                                 aTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTcaca                  <  1:19604/148‑1 (MQ=255)
                                                                                                                                  tAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtt                >  4:764614/1‑149 (MQ=255)
                                                                                                                                    aCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtttt              <  4:884084/149‑1 (MQ=255)
                                                                                                                                    aCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAttt               <  1:616460/148‑1 (MQ=255)
                                                                                                                                       ggCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtttttt            <  1:349653/148‑1 (MQ=255)
                                                                                                                                          cAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGa        <  4:747363/149‑1 (MQ=255)
                                                                                                                                           ggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAt                 <  4:625842/138‑1 (MQ=255)
                                                                                                                                           ggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGat       >  4:237338/2‑149 (MQ=255)
                                                                                                                                            gCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGatt      >  4:104990/1‑149 (MQ=255)
                                                                                                                                             cgtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGagta     <  4:330470/149‑1 (MQ=255)
                                                                                                                                               tcgtGATGACGACGCCTTTTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattatt   <  3:736842/149‑1 (MQ=255)
                                                                                                                                               tcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattatt   <  3:47623/149‑1 (MQ=255)
                                                                                                                                                 gtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTg  <  2:73480/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTG  >  NC_000913/3512159‑3512451

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