Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 20.2% A97P (GCA→CCA)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5890GC20.2% 362.3 / 88.9 153A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (57/65);  new base C (13/18);  total (70/83)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 169 total aligned reads displayed.

GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252490‑252688
                                                                                                   |                                                                                                   
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                     <  1:2269744/100‑1 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                    <  1:3621471/100‑1 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                    >  1:3438380/1‑100 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                    >  1:3333232/1‑100 (MQ=255)
       agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGtt                                                                                              <  1:2958729/100‑1 (MQ=255)
       agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATgctt                                                                                              >  2:2199092/1‑97 (MQ=255)
        gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTGTATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                             >  1:2773186/1‑100 (MQ=255)
         aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCTGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                            <  1:2770175/100‑1 (MQ=255)
           taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                          <  2:1989208/100‑1 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  2:3082791/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  2:552621/100‑2 (MQ=255)
               ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  1:839734/99‑2 (MQ=255)
                cATTTCCGCCACGGTTGCGCCAGTCCAGGATAGCGGCTGACGAAGCGCATTTGTTATATTTTCAGGGGGCGTTTGATGAGCTTGCATGGCTGCTAATTgg                                                                                     <  2:221794/100‑2 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAATCGCATTTGGTATATTGTCAGGGGGCGTTTTATGAGCAACCATGCTTGCTCAtca                                                                                      <  1:182198/99‑3 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     <  1:3530694/100‑1 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGGCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                      <  2:2100371/99‑1 (MQ=255)
                 aTTTCCGCCACAGTTGCGGCAATTCAGTAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  2:3700701/100‑5 (MQ=38)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  1:1669032/100‑5 (MQ=39)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  2:1851826/100‑5 (MQ=39)
                  tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTACTCATTGCCa                                                                                   >  1:1873921/1‑100 (MQ=255)
                  tgtCCGCCACAGTTGCGCCAATTCATCAGAACGAACGACGAAGCGCATTGAGTATATTGTCAGCGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                   <  2:1361606/98‑6 (MQ=37)
                   ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                  <  2:1527587/100‑1 (MQ=255)
                   ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTAATGAGCAACCATGCTTGCTCAtcaaacg                                                                                  <  2:2850507/100‑7 (MQ=38)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                <  2:2329715/100‑1 (MQ=255)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                <  2:658792/100‑1 (MQ=255)
                      cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                               >  2:263693/1‑100 (MQ=255)
                      cGCCACAGTTGAGCCAATTCAGCAGAATGAATGACGAAGCGCATTTGGTATATTGTTAGGGGGCGTTTGATGAGCAAGGATGGTTGCTCATTGGCATTCt                                                                               >  2:2793024/1‑100 (MQ=255)
                          acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                           <  1:3250050/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGTATGATGAAGCGCATTTGGCCTATTGTCAGGGGGCGTTTGATGAGCAAGCATCGTTGCTCATTGCCATTCTGAAAc                                                                          <  2:1311068/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                          <  1:3762712/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                          <  1:3359436/100‑1 (MQ=255)
                             tctGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCt                                                                        <  1:684986/98‑1 (MQ=255)
                             gTTGCGCCAATTCAGCAGAACGAATGACGAAGCGTATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCt                                                                        >  2:2140314/1‑100 (MQ=255)
                             gTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCt                                                                        >  1:2287075/1‑100 (MQ=255)
                                gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCtaga                                                                     >  2:684982/1‑97 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    >  1:2773110/1‑100 (MQ=255)
                                  gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAccc                                                                   >  2:522562/1‑100 (MQ=255)
                                  gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                   <  2:610448/100‑1 (MQ=255)
                                      aTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCAt                                                               <  2:2245255/100‑1 (MQ=255)
                                       ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                              <  1:704491/100‑1 (MQ=255)
                                       ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                              >  2:372459/1‑100 (MQ=255)
                                       ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                              <  1:333536/100‑1 (MQ=255)
                                        tCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaa                                                             >  2:2976757/1‑100 (MQ=255)
                                        tCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaa                                                             >  2:2535782/1‑100 (MQ=255)
                                         cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAGGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                            >  2:3804118/1‑100 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGTCGTTTGATCAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          <  1:210766/100‑1 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          >  2:2632120/1‑100 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                          <  2:3067019/100‑1 (MQ=255)
                                            cagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGATTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                         >  1:89558/1‑100 (MQ=255)
                                             agAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                        <  1:1152776/100‑1 (MQ=255)
                                              gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                       >  1:1882237/1‑100 (MQ=255)
                                               aaCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCt                                                      <  1:658394/100‑1 (MQ=255)
                                                aCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                     <  2:2034942/100‑1 (MQ=255)
                                                       aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCACGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAATTGGCTCGAGATaa                                              >  2:2349235/1‑100 (MQ=255)
                                                       aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                              <  1:148086/100‑1 (MQ=255)
                                                       aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                              >  1:3414199/1‑100 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAGTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                            <  2:2421534/100‑1 (MQ=255)
                                                            gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  2:2868291/100‑1 (MQ=255)
                                                                aTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                     >  1:837943/1‑100 (MQ=255)
                                                                aTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                     <  1:3200174/100‑1 (MQ=255)
                                                                 tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATg                                    <  1:3468494/100‑1 (MQ=255)
                                                                  ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATg                                    <  2:1184477/99‑1 (MQ=255)
                                                                  ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                   <  2:3239239/100‑1 (MQ=255)
                                                                   tGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCa                                  >  1:2196634/1‑100 (MQ=255)
                                                                    ggTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCAt                                 <  1:3028937/100‑1 (MQ=255)
                                                                      tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                               >  2:1481148/1‑100 (MQ=255)
                                                                      tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                               >  1:2127718/1‑100 (MQ=255)
                                                                        tatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTAGCTCGAGATAACAACCAAATGCATaaaa                             >  2:1943599/1‑100 (MQ=255)
                                                                         atTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAt                            <  1:627354/100‑1 (MQ=255)
                                                                          ttGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAtt                           >  1:407385/1‑100 (MQ=255)
                                                                             tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                        <  1:263693/100‑1 (MQ=255)
                                                                               aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                      >  1:2143138/1‑100 (MQ=255)
                                                                                gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  1:2952517/1‑100 (MQ=255)
                                                                                 ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                    >  2:1323001/1‑100 (MQ=255)
                                                                                  gggCGGTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                       >  2:732112/1‑96 (MQ=255)
                                                                                    gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                 <  1:2050012/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                <  2:11176/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                <  1:757244/100‑1 (MQ=255)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:3433315/7‑100 (MQ=38)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:3336257/7‑100 (MQ=38)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCCCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  1:999279/7‑100 (MQ=38)
                                                                                      gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               <  2:3403237/100‑1 (MQ=255)
                                                                                      gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               <  1:459838/100‑1 (MQ=255)
                                                                                       ggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAg              >  1:2632612/5‑100 (MQ=39)
                                                                                        ttGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  1:2888846/1‑100 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  1:866052/4‑100 (MQ=255)
                                                                                         tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:596803/1‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:299848/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:1571895/3‑100 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  1:1715727/2‑100 (MQ=255)
                                                                                             gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGt        <  1:3461447/100‑1 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt       <  2:3616271/100‑1 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  1:3102499/1‑100 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  1:2339177/1‑100 (MQ=255)
                                                                                                caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACGGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc     <  1:3445570/100‑1 (MQ=255)
                                                                                                 aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg    >  2:2018640/1‑100 (MQ=255)
                                                                                                  agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc   >  1:1777815/1‑100 (MQ=255)
                                                                                                   gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAACCTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGcn  <  2:394753/100‑2 (MQ=255)
                                                                                                   gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  <  1:1556824/100‑1 (MQ=255)
                                                                                                   gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  <  1:1037902/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252490‑252688

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