Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,410,558 C→G 6.1% W265C (TGG→TGC) ‡ yfbS ← putative transporter YfbS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,410,5580CG 6.1% 220.8 / 4.8 66W265C (TGG→TGC) ‡yfbSputative transporter YfbS
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAGGT  >  NC_000913/2410409‑2410707
                                                                                                                                                     |                                                                                                                                                     
tCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAAcgccgc                                                                                                                                                       >  1:1806084/1‑150 (MQ=255)
     gcgcAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCg                                                                                                                                                  <  2:1052963/150‑1 (MQ=255)
        cAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTc                                                                                                                                               <  1:283295/150‑1 (MQ=255)
        cAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTc                                                                                                                                               >  2:1168170/1‑150 (MQ=255)
             ccATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCTGGCGGCGttttcgtc                                                                                                                                          <  1:528566/150‑7 (MQ=255)
             ccATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAacac                                                                                                                                          >  1:731210/1‑150 (MQ=255)
             ccATCGGCTCCAGCAATTGCTCACTACAAAATTGCCAGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACTCTGGCGGCGttttcg                                                                                                                                            <  1:623092/148‑5 (MQ=255)
             ccATCGGCTCCAGCAATTGCTCACTACAAAATTGCCAGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCTGGCGGCGttttcg                                                                                                                                            >  2:623092/1‑144 (MQ=255)
                 cGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGAt                                                                                                                                      <  1:1919198/150‑1 (MQ=255)
                     tCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACg                                                                                                                                  <  2:994609/150‑1 (MQ=255)
                     tCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACg                                                                                                                                  <  1:613286/150‑1 (MQ=255)
                     tCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACg                                                                                                                                  <  2:633269/150‑1 (MQ=255)
                     tCCAGCAATTACTCACTACAAAATTTCCGGATATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCTCGCGAAATTCAGAATCCCCATTAACGTTCACGATAACGCGAAGAAAACGCCGCCAGCGTTCAACACCGATGACg                                                                                                                                  <  2:1327565/150‑1 (MQ=255)
                      ccAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGt                                                                                                                                 >  1:1885057/1‑150 (MQ=255)
                       cAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGtt                                                                                                                                <  1:1808254/150‑1 (MQ=255)
                        aGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTa                                                                                                                               >  1:1177004/1‑150 (MQ=255)
                        aGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTa                                                                                                                               >  1:332864/1‑150 (MQ=255)
                                    cTACAAAATTGCCGGAGATCGACGTCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGc                                                                                                                   >  1:1961477/1‑150 (MQ=255)
                                      aCAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                                                                 <  2:1518182/150‑1 (MQ=255)
                                         aaaTTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACg                                                                                                              >  2:1042448/1‑150 (MQ=255)
                                           aTTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTa                                                                                                            <  2:1640889/150‑1 (MQ=255)
                                             tGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAAt                                                                                                          >  1:860662/1‑150 (MQ=255)
                                              gCCGGAGATCGACATCAGCCGCAGACATATCAATAGGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAAtt                                                                                                         >  2:1983007/1‑150 (MQ=255)
                                              gCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAAtt                                                                                                         >  2:1794352/1‑150 (MQ=255)
                                                cGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACAcc                                                                                                                                         <  1:1331960/116‑1 (MQ=255)
                                                cGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACAcc                                                                                                                                         >  2:1331960/1‑116 (MQ=255)
                                                cGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTg                                                                                                       >  1:367225/1‑150 (MQ=255)
                                                 ggAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTga                                                                                                      >  2:1153042/1‑150 (MQ=255)
                                                    gaTCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGat                                                                                                   >  1:1174948/1‑150 (MQ=255)
                                                    gaTCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGat                                                                                                   <  1:1856421/150‑1 (MQ=255)
                                                    gaTCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGat                                                                                                   >  1:1436622/1‑150 (MQ=255)
                                                         aCATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATcc                                                                                              <  2:1367420/150‑1 (MQ=255)
                                                                 cGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTg                                                                                      <  2:860662/150‑1 (MQ=255)
                                                                  gCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGa                                                                                     >  2:1982778/1‑150 (MQ=255)
                                                                   cAGACATATCAATAAGCAAAACATCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGAc                                                                                    >  1:438762/1‑150 (MQ=255)
                                                                        atatCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTCCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATc                                                                               <  2:1883092/150‑1 (MQ=255)
                                                                         tatCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCa                                                                              <  1:1794352/150‑1 (MQ=255)
                                                                         tatCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGAAGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCa                                                                              >  1:422927/1‑150 (MQ=255)
                                                                          atCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGCTCATCCAGCCGTTGACCAATCAt                                                                             >  2:1471707/1‑150 (MQ=255)
                                                                             aaTAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTgg                                                                          <  1:1572601/150‑1 (MQ=255)
                                                                             aaTAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTgg                                                                          >  1:1445558/1‑150 (MQ=255)
                                                                              aTAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTgg                                                                          >  1:1118514/1‑149 (MQ=255)
                                                                              aTAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTgg                                                                          <  2:1118514/149‑1 (MQ=255)
                                                                              aTAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGc                                                                         <  2:367225/150‑1 (MQ=255)
                                                                                 aGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGAt                                                                      >  2:525904/1‑150 (MQ=255)
                                                                                            aCGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAAt                                                           <  2:1229922/150‑1 (MQ=255)
                                                                                                gcgcgAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGcc                                                       <  2:1406938/150‑1 (MQ=255)
                                                                                                   cgAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGa                                                    <  2:967357/150‑1 (MQ=255)
                                                                                                   cgAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGa                                                    >  1:1568258/1‑150 (MQ=255)
                                                                                                        ttCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTcg                                               <  2:871969/150‑1 (MQ=255)
                                                                                                          cAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTcgcg                                             >  1:602184/1‑150 (MQ=255)
                                                                                                          cAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTcgcg                                             <  1:1681706/150‑1 (MQ=255)
                                                                                                            gAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                                           <  2:1806084/150‑1 (MQ=255)
                                                                                                            gAAACACCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGCTCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCa                                           >  1:1212255/1‑150 (MQ=255)
                                                                                                             aaaCCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCAc                                          >  2:484734/1‑150 (MQ=255)
                                                                                                                          gTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCg                             >  1:124326/1‑150 (MQ=255)
                                                                                                                            tCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGat                           <  2:1961477/150‑1 (MQ=255)
                                                                                                                             cACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGata                          >  1:1312726/1‑150 (MQ=255)
                                                                                                                                   aaCGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACg                    <  1:1052760/150‑1 (MQ=255)
                                                                                                                                     cgcgACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGa                  <  2:1568258/150‑1 (MQ=255)
                                                                                                                                     cgcgACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGa                  >  1:87136/1‑150 (MQ=255)
                                                                                                                                      gcgACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGAt                 <  1:1982778/150‑1 (MQ=255)
                                                                                                                                       cgacgaAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATa                <  2:332864/150‑1 (MQ=255)
                                                                                                                                       cgacgaAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATa                <  1:1930961/150‑1 (MQ=255)
                                                                                                                                          gttgaaCGCTGGCGGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGa             >  2:1564176/5‑150 (MQ=255)
                                                                                                                                             aaaCGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATcg          <  1:155539/150‑1 (MQ=255)
                                                                                                                                               aCGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATcgcg        <  1:914680/150‑1 (MQ=255)
                                                                                                                                                   cgccAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAg    <  1:484734/150‑1 (MQ=255)
                                                                                                                                                   cgccAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAg    >  1:1997822/1‑150 (MQ=255)
                                                                                                                                                    gccAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAgg   <  2:731210/150‑1 (MQ=255)
                                                                                                                                                     ccAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAGGt  >  2:683759/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCGCGCAGTACCATCGGCTCCAGCAATTGCTCACTACAAAATTGCCGGAGATCGACATCAGCCGCAGACATATCAATAAGCAAAACGTCACGCGCGCGAAATTCAGAAACCCCATTAACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACGTCGCGCACGCCCGGTCAGTCGATATTCACGGATAAGATCGCGAAAGGT  >  NC_000913/2410409‑2410707

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