Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,853 T→C 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,8530TC95.8% 377.5 / 2.8 119N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (38/75);  total (41/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 119 total aligned reads displayed.

GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949
                                                                                                 |                                                                                                
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  1:222756/1‑100 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCtn                                                                                              <  1:1297522/100‑2 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                              <  1:684843/100‑1 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                              <  2:1430737/100‑1 (MQ=255)
   aaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCa                                                                                             <  2:1557760/100‑1 (MQ=255)
    aTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAt                                                                                            >  2:1178201/1‑100 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         <  1:1587587/100‑1 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         >  2:1525693/1‑100 (MQ=255)
        aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACa                                                                                        >  2:1238748/1‑100 (MQ=255)
         aaaaaTATCCAAAAATAGTATACTTTATTCCATAATAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAt                                                                                       <  1:1072877/100‑1 (MQ=255)
         aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAt                                                                                       <  1:183408/100‑1 (MQ=255)
          aaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAtt                                                                                      >  1:450968/1‑100 (MQ=255)
          aaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAtt                                                                                      <  2:1811133/100‑1 (MQ=255)
           aaaTATCCAAAGATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAtn                                                                                      <  1:895618/99‑2 (MQ=255)
           aaaTATCCAAAAATAGTATACTTTATTTCAGAAGAGTTCAATATAATGTTTGTCTTTAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtg                                                                                     <  2:1154948/100‑1 (MQ=255)
           aaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtg                                                                                     <  2:607116/100‑1 (MQ=255)
            aaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATCATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgt                                                                                    <  1:1162639/100‑1 (MQ=255)
              tatCCAAAAATAGTATACTTTATTCCAGAAGATTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                  <  2:8450/100‑1 (MQ=255)
              tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                  <  2:833075/100‑1 (MQ=255)
                tCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGc                                                                                <  2:950987/100‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCTAGAAGAGTTCAATATGATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              >  1:966568/1‑100 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              <  1:719547/100‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAAGGTAATATAAATCGCTCATTACATTGTGAGCt                                                                               >  2:61748/1‑99 (MQ=255)
                    aaaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCACTTTTTCTTAATTAAGGGTAATATAAATAGATCATTAAATTCTGAACATCa                                                                            >  2:149039/1‑100 (MQ=255)
                        tAGTATACTTTATTCTAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtt                                                                        <  2:1787948/100‑1 (MQ=255)
                          gTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCttta                                                                      <  1:418611/100‑1 (MQ=255)
                           tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                     <  2:558211/100‑1 (MQ=255)
                            ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                    <  2:479015/100‑1 (MQ=255)
                               ccgatctCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                                 <  1:1200430/94‑1 (MQ=255)
                               cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                                 >  1:721485/1‑100 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                                 <  2:1197668/99‑1 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                <  1:1262761/100‑1 (MQ=255)
                                  tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTATATTGTGAGCTTCATCTTTATTTAAttt                                                              <  1:1496672/100‑1 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                            <  1:1325458/100‑1 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  1:560607/100‑1 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  1:903951/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  1:32934/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  1:1073976/1‑100 (MQ=255)
                                       cAGAAGAGTTCAATATAATTTTTGTTTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         <  2:32934/100‑1 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTcaga                                                         >  2:1181993/1‑97 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         <  1:1759673/100‑1 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         >  2:1005427/1‑100 (MQ=255)
                                        agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                        <  1:867760/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  2:1428342/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  1:1178190/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  2:709975/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       >  2:948368/1‑100 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  1:1066100/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATAGTGAACTTCATCTTTATTTAATTTTCTGTTg                                                       >  2:123346/1‑100 (MQ=255)
                                            gagTTTAATATTATGTTTGTTTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACt                                                    <  2:1704282/100‑1 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTACAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACt                                                    >  1:970392/1‑100 (MQ=255)
                                              gTTCAATGTAATGTTTTTCTTCTATTTTTTTTTCTTTTGGGTATTTTAGATCTTTCTTTTCATTTTTTGCTTCATCTTTATTTAATTGTCTGTTGACt                                                    <  1:590734/98‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  <  2:1593094/100‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTc                                                   >  2:1532636/1‑99 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  2:1211915/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:1085710/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  2:348525/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:997973/100‑1 (MQ=255)
                                                  aaTATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGct                                              >  1:1240709/1‑100 (MQ=255)
                                                   atatAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGctc                                             <  1:1710402/100‑1 (MQ=255)
                                                   atatAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGct                                              >  1:328159/1‑99 (MQ=255)
                                                      taATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTcc                                          <  1:1668833/100‑1 (MQ=255)
                                                       aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                         <  2:1385277/100‑1 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                        <  2:450970/100‑1 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                        <  2:1271869/100‑1 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                        <  1:696614/100‑1 (MQ=255)
                                                         tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                       <  2:298195/100‑1 (MQ=255)
                                                         tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                       >  1:734602/1‑100 (MQ=255)
                                                          gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                      <  2:1455725/100‑1 (MQ=255)
                                                          gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGATCTCCGTg                                       <  1:558712/99‑1 (MQ=255)
                                                           tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                     <  2:488442/100‑1 (MQ=255)
                                                           tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                     <  1:219876/100‑1 (MQ=255)
                                                            ttGTCTTTAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATa                                    <  1:277345/100‑1 (MQ=255)
                                                             tGTCTCTAATTTTTCTTACTACAGGGTCATATAGATCGCTCACTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATat                                   <  1:1710506/100‑2 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                  <  2:719229/100‑1 (MQ=255)
                                                               tCTTCAATTTTTCTTTCTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCATCTCTCCGTAATAAc                                  <  1:879931/99‑1 (MQ=255)
                                                               tCTTCAATTTTTCTTACTTCAGGGTAATATATATCGCTCATTACATTGTGAGCTGCATCTTTATTTAATTTTTTGTTGACTCCAGCTGTCCGTGATaat                                  <  1:382063/99‑2 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                >  1:410611/1‑100 (MQ=255)
                                                                   cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             >  2:121465/1‑100 (MQ=255)
                                                                   cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  1:1525682/100‑1 (MQ=255)
                                                                       ttttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtata                         <  1:1757444/100‑1 (MQ=255)
                                                                          tCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAGTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtt                      <  1:527014/100‑1 (MQ=255)
                                                                           cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAATGGTTTTATAATTa                     >  2:1495158/1‑100 (MQ=255)
                                                                             tACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATGTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     <  1:985220/98‑1 (MQ=255)
                                                                              aCTTCAGGGTAATATAGATCGTTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                  <  2:194807/100‑1 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 <  1:894313/100‑1 (MQ=255)
                                                                                 tCAGGGTAATATAGATTGCTCCTTACATTGTGAGCTTCATCTTTATTTAATTTTGTGTTGACTCCAGCTCTCCGTGATAATGGTTTGATAATTAGATgtt               <  1:149039/100‑3 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:461157/1‑100 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  2:12430/1‑100 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:773060/1‑100 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:1559543/1‑100 (MQ=255)
                                                                                      gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  2:477504/100‑1 (MQ=255)
                                                                                         atatAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa       >  2:825811/1‑100 (MQ=255)
                                                                                         atatAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa       >  2:1613257/1‑100 (MQ=255)
                                                                                           ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  2:1692157/100‑1 (MQ=255)
                                                                                           ataGATCGCTCACTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     >  2:1195915/1‑100 (MQ=255)
                                                                                            taGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa    <  1:585606/100‑1 (MQ=255)
                                                                                            taGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa    >  2:1345820/1‑100 (MQ=255)
                                                                                             aGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGat   >  2:1249650/1‑100 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  2:749328/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949

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