Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 17.3% 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,5890GC17.3% 530.6 / 97.3 204A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (72/95);  new base C (19/16);  total (93/111)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 226 total aligned reads displayed.

GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCA  >  NC_000913/252490‑252688
                                                                                                   |                                                                                                   
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCATTTGATGagcaag                                                                                                     >  1:1708871/1‑100 (MQ=255)
    ctGAGAATAAACGATTTCCGCAACAGCGGCGCGAGTTGAGCAGTACGATTGACGAATTGCCTTTTGTATCTTGTCAGGGGGCGTTTGATGTGCAAGCAt                                                                                                  >  1:4364201/2‑99 (MQ=255)
    aTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATg                                                                                                 >  1:758129/1‑100 (MQ=255)
    aTGAGAATAATCCATTGCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATCTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATg                                                                                                 <  1:4141896/100‑1 (MQ=255)
     tccgaATACTCCACTTCCGCCACCGTTGCCCCAAATCAGCAGAGCGACTGACGAAGCGCATTTTGTATATTGTTAGGTGGTGTTGGATGAGCGAGCAggg                                                                                                >  1:1903588/4‑100 (MQ=255)
     tGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                >  2:2380026/1‑100 (MQ=255)
       agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGtt                                                                                              <  2:1707436/100‑1 (MQ=255)
         aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGTATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                            <  1:2412117/100‑1 (MQ=255)
           taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTc                                                                                          <  1:2810014/100‑1 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                        <  2:2763263/100‑1 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                        >  1:4346797/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  2:747426/100‑1 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                        <  1:2779320/100‑1 (MQ=255)
              tgcATTTCCGCCACAGTTGCGTCAATTAAGCAGATCGGATGACGAAGCGCATTTGGTATATTGTCAGTGGGGGTTTGATGGGCAGCCATGGTTGCTCAt                                                                                        <  1:3246398/97‑1 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                       <  2:1134030/100‑2 (MQ=255)
               ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtca                                                                                      <  2:3641916/100‑3 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     >  2:1615394/1‑100 (MQ=255)
                cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                     <  2:2768192/100‑1 (MQ=255)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                    <  1:11987/100‑1 (MQ=255)
                 aTTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaa                                                                                    <  2:3813510/100‑5 (MQ=39)
                   ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                  <  2:909101/100‑1 (MQ=255)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                >  2:1575937/1‑100 (MQ=255)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTAATGAGCAACCATGCTTGCTCAtcaaacgcc                                                                                <  2:2658599/100‑9 (MQ=37)
                       gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                              <  1:2743735/100‑1 (MQ=255)
                        ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGa                                                                             >  1:2314990/1‑100 (MQ=255)
                         cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGa                                                                             <  2:3384970/99‑1 (MQ=255)
                          acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                           <  2:4366878/100‑1 (MQ=255)
                            aGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAcc                                                                         <  1:1862920/100‑1 (MQ=255)
                               tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                      <  2:1457685/100‑1 (MQ=255)
                                 ngacAATTCAGCAGAACGAATGACGAAGCGAATGTGGTATTTTGTCAGGAGGCGTATGATGAGCAAACATGATTGGTCATTGTCATTCTGAAAACTGaa                                                                     >  1:965567/4‑99 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    >  1:268179/1‑100 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    <  2:3059515/100‑1 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTAGTATATTGTCAGGGGGCGTTTGATTAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAc                                                                    <  1:454936/100‑1 (MQ=255)
                                  gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                   <  2:3890655/100‑1 (MQ=255)
                                   ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                  <  2:3977442/100‑1 (MQ=255)
                                    cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                 <  2:3583548/100‑1 (MQ=255)
                                     aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                <  1:2355488/100‑1 (MQ=255)
                                     aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                <  1:2123853/100‑1 (MQ=255)
                                         nagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                            >  1:2231394/2‑100 (MQ=255)
                                          agccgAACGAATGACAAAACGCATCTGATAAATTCGCAGGGGGCGTTTGATGAGCAACCATGGCTGCTCATTGCCATTCTGAAACCTGAACCTCAt                                                               >  1:835306/1‑96 (MQ=255)
                                               aaCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCt                                                      <  2:2927011/100‑1 (MQ=255)
                                                 cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  2:2465396/100‑1 (MQ=255)
                                                 cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                    <  1:1417948/100‑1 (MQ=255)
                                                  gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                   <  1:2502639/100‑1 (MQ=255)
                                                   aaTGACGAAGCGCATTGGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                  <  1:3339170/100‑1 (MQ=255)
                                                    aTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                                 >  2:2553082/1‑100 (MQ=255)
                                                     tGACGAAGCGCATTTGGTATATTGTCAGGGGGCTTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                <  1:3696152/100‑1 (MQ=255)
                                                     tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                <  2:3336594/100‑1 (MQ=255)
                                                      gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                               <  2:661182/100‑1 (MQ=255)
                                                       aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                              <  2:1657546/100‑1 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                             <  1:3423823/99‑1 (MQ=255)
                                                           aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaac                                          >  2:4158489/1‑100 (MQ=255)
                                                            gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  2:496772/100‑1 (MQ=255)
                                                            gcgcATTTCGTATATTGTCAGGGGGCGTTTGATCAGCAAGAATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                         <  1:3013834/100‑1 (MQ=255)
                                                              gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAGACCTGGACCTCATACACTGGCTCGAGATAACAACCaa                                       >  1:2972678/1‑100 (MQ=255)
                                                               cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                      <  1:509434/100‑1 (MQ=255)
                                                                 tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATg                                    <  2:4035794/100‑1 (MQ=255)
                                                                 tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATg                                    >  2:4092763/1‑100 (MQ=255)
                                                                  ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                   >  2:116318/1‑100 (MQ=255)
                                                                    ggTATATTGTCAGGGGTCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCAt                                 >  1:3403527/1‑100 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                >  1:1977339/1‑100 (MQ=255)
                                                                       atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                >  1:203535/1‑98 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAGCTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         <  1:3157555/100‑1 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgn                         <  2:2231403/100‑2 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         <  2:4299667/100‑1 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         <  1:649611/100‑1 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                         >  2:4017618/1‑100 (MQ=255)
                                                                             tccgggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                        >  1:1076960/4‑100 (MQ=255)
                                                                               aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                      <  2:1330756/100‑1 (MQ=255)
                                                                                gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  2:3937250/1‑100 (MQ=255)
                                                                                gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaa                     >  1:1612536/1‑100 (MQ=255)
                                                                                   ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATg                  <  2:1256436/100‑1 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:3710162/1‑100 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:3516054/1‑100 (MQ=255)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  1:1275382/7‑100 (MQ=38)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                >  2:3194189/7‑100 (MQ=38)
                                                                                     atggcaatgaGCAACCATGCTTGCTCATTGCCATTATGAAACCTGAACCTCATAAACTGGCTAGAGATAACAACCAAATGCACAAAATTGGGACAATgg                 >  1:994983/7‑99 (MQ=37)
                                                                                      gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa               >  2:1065667/1‑100 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  2:1433983/4‑100 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  2:3320505/4‑100 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa             >  1:4286046/4‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCTTAAAATTGGGAAAATGGCAGaa            >  1:3481812/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  2:1840593/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:2524931/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa            >  2:2571457/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCAAGATAACAAACAAATGCATAAAATTGGGAAAATGGCAGaa            >  1:256675/3‑100 (MQ=255)
                                                                                          gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           <  2:2270219/100‑1 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  2:3618568/2‑100 (MQ=255)
                                                                                          aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg           >  2:235513/2‑100 (MQ=255)
                                                                                           atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc          >  1:662087/1‑100 (MQ=255)
                                                                                           atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc          >  1:2067267/1‑100 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt       >  2:2873094/1‑100 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt       <  2:1301982/100‑1 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCGTAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGt        >  1:587474/1‑98 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  2:2164480/1‑100 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  2:2108494/1‑100 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  1:2849328/1‑100 (MQ=255)
                                                                                               gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt      >  1:817131/1‑100 (MQ=255)
                                                                                                   gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  <  1:4037374/100‑1 (MQ=255)
                                                                                                   gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAAATAAAAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa  >  1:3821141/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: