Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,853 T→C 63.5% N129D (AAT→GAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,8530TC63.5% 134.3 / 141.3 145N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/31);  new base C (38/54);  total (60/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 151 total aligned reads displayed.

TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATATCG  >  NC_000913/1197754‑1197952
                                                                                                   |                                                                                                   
   aaagTTCAAAAAATATCCAAAAATTGTATAATTTATTTCAAAAGAGTGTAATAAATTGTTTGCCTTCAATTTTTTTTACTTCAGGATATTATAGAtttat                                                                                                  >  1:2023437/1‑97 (MQ=255)
    aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                                 >  2:2724357/1‑100 (MQ=255)
    aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                                 <  2:103421/100‑1 (MQ=255)
        tCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCACTCATTa                                                                                             >  2:4183592/1‑100 (MQ=255)
          aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACa                                                                                           >  2:3525777/1‑100 (MQ=255)
           aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAt                                                                                          >  1:3633193/1‑100 (MQ=255)
           aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTCCAGGGTAATATAGATCGCTCATTACAt                                                                                          >  1:3633072/1‑100 (MQ=255)
             aaaTATACAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtg                                                                                        <  1:2499286/100‑1 (MQ=255)
                tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGa                                                                                     <  1:355880/100‑1 (MQ=255)
                tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                     <  2:2013304/100‑1 (MQ=255)
                 atCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAg                                                                                    <  1:3476255/100‑1 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCt                                                                                  <  2:1193454/100‑1 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCt                                                                                  >  1:561485/1‑100 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCt                                                                                  <  1:2592069/100‑1 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCt                                                                                  <  2:1871547/100‑1 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCtt                                                                                 <  1:1291921/100‑1 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                                 >  2:616634/1‑100 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                                 >  2:1230009/1‑100 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                                 >  2:3877341/1‑100 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTGCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCtt                                                                                 >  1:2825835/1‑100 (MQ=255)
                      aaaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCa                                                                               >  1:1293061/1‑100 (MQ=255)
                          tAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtt                                                                           >  1:3886527/1‑100 (MQ=255)
                             ttttttttATTCCAGTCGAGTTCAATATAATGTTTGTCTTAAATTTTTATTACTTCGGGGTAATATATATCGCTCATTACATTGTGATCTTCATCgttat                                                                        <  1:4126148/95‑1 (MQ=255)
                             tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                        <  2:837755/100‑1 (MQ=255)
                              ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttatt                                                                       >  1:29367/1‑100 (MQ=255)
                                aCTTTATTCCAGAAGAGTTCAATCTAATGTTTGTATTCACTTTTTCTTACTTCTGGGTAATATACATCGCGCCTTTCATTGTGAGCTTCATCtttattta                                                                     <  1:1465464/100‑1 (MQ=255)
                                 cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATGGATCGATCATTACATTGTGAGCTTCATCTTTATTTaa                                                                    >  1:2982563/1‑100 (MQ=255)
                                  tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                   >  2:2942860/1‑100 (MQ=255)
                                  tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                   <  2:2120548/100‑1 (MQ=255)
                                    tattttagaagaGTTCAATATAATGTTTGTCTTTAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                                 <  1:491683/94‑1 (MQ=255)
                                    tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                                 >  2:4210997/1‑100 (MQ=255)
                                    tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                                 >  1:2214630/1‑100 (MQ=255)
                                    tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                                 <  1:3769232/100‑1 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                             >  2:1445531/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                             <  1:1461469/99‑1 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                            >  2:3756065/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTCTCTGt                                                            >  2:3686905/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCTCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                            <  1:1881799/100‑1 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                            >  1:1796359/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                            <  2:3287372/100‑1 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                          <  1:2947036/100‑1 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                           >  1:1802444/1‑99 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                          <  1:532514/100‑1 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                          <  1:106362/100‑1 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                          >  1:4334967/1‑100 (MQ=255)
                                                gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                     <  2:89905/100‑1 (MQ=255)
                                                  tCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                                   >  1:4362846/1‑100 (MQ=255)
                                                  tCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                                   <  1:1016875/100‑1 (MQ=255)
                                                  tCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                                   <  2:288685/100‑1 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                                  >  2:4302013/1‑100 (MQ=255)
                                                       ataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCtctc                                              >  2:3450235/1‑100 (MQ=255)
                                                         aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                            >  1:1120326/1‑100 (MQ=255)
                                                         aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                            <  1:4060973/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                           <  2:4249772/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                           <  2:1435138/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                           <  2:1100721/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                           <  2:3860163/100‑1 (MQ=255)
                                                           tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                          <  1:4201767/100‑1 (MQ=255)
                                                            gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                         >  2:183006/1‑100 (MQ=255)
                                                            gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                         <  2:2291229/100‑1 (MQ=255)
                                                             tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                        <  1:4167087/100‑1 (MQ=255)
                                                             tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                        <  2:4334999/100‑1 (MQ=255)
                                                             tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                        >  2:3114781/1‑100 (MQ=255)
                                                               tGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATaa                                      <  1:3525753/100‑1 (MQ=255)
                                                                gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                     <  2:3006537/100‑1 (MQ=255)
                                                                gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                     >  2:1170976/1‑100 (MQ=255)
                                                                gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                     <  1:3041631/100‑1 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                   >  2:337979/1‑100 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                   >  1:2120411/1‑100 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTATTGACTCCAGCTCTCCGTGATAActg                                   >  2:1267604/1‑98 (MQ=255)
                                                                   ttCAATTTTTCTTACTTCAGGGTAATATATATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCAGCTCTCCGTGATACCGGtt                                 <  1:4186492/100‑1 (MQ=255)
                                                                   ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                                  >  2:3197354/1‑100 (MQ=255)
                                                                      aaTTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGCCTCCAGCTCTCCGTGATAACGGtttt                               >  1:1400927/1‑100 (MQ=255)
                                                                       aTTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTta                              <  1:995873/100‑1 (MQ=255)
                                                                       aTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTta                              <  2:2617357/100‑1 (MQ=255)
                                                                        tttttCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtat                             <  2:4252227/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATATATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                        <  1:3460955/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                        <  2:589589/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                        <  1:2808821/100‑1 (MQ=255)
                                                                                aCTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                     <  1:3313376/100‑1 (MQ=255)
                                                                                aCTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                     <  1:3949216/100‑1 (MQ=255)
                                                                                 cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                    >  1:2633280/1‑100 (MQ=255)
                                                                                  ttCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                   >  1:2369697/1‑100 (MQ=255)
                                                                                  tgcAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                   >  2:2995532/3‑100 (MQ=255)
                                                                                   tCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt                  <  2:2274297/100‑1 (MQ=255)
                                                                                   tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt                  <  2:1277319/100‑1 (MQ=255)
                                                                                     aGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGCCTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa                <  2:913744/100‑1 (MQ=255)
                                                                                       ggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATc              <  2:1082590/100‑1 (MQ=255)
                                                                                        gTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcn             <  2:1976855/100‑2 (MQ=255)
                                                                                        gTAATATAGATTGCTCATTACATTGTCAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc             <  1:1088331/100‑1 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc             <  1:4301981/100‑1 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc             <  2:1178622/100‑1 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc             <  2:1968192/100‑1 (MQ=255)
                                                                                             ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg        <  1:4251307/100‑1 (MQ=255)
                                                                                             ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg        <  2:1037751/100‑1 (MQ=255)
                                                                                             ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg        >  1:1376457/1‑100 (MQ=255)
                                                                                             ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg        <  2:2716144/100‑1 (MQ=255)
                                                                                                  ttGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATATc   <  1:3231704/100‑1 (MQ=255)
                                                                                                   tGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATATCg  <  1:709411/100‑1 (MQ=255)
                                                                                                   tGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATATCg  <  2:2598799/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATATCG  >  NC_000913/1197754‑1197952

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