Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,873 C→T 5.2% intergenic (‑23/‑61) dicC ← / → dicA Qin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,647,8730CT 5.2% 324.4 / 5.8 97intergenic (‑23/‑61)dicC/dicAQin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB
Reads supporting (aligned to +/- strand):  ref base C (43/48);  new base T (3/2);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCC  >  NC_000913/1647727‑1648021
                                                                                                                                                  |                                                                                                                                                    
gCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                    <  4:1158684/149‑1 (MQ=255)
gCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                    >  2:462165/1‑149 (MQ=255)
gCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                    >  2:1135147/1‑149 (MQ=255)
  aTCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaa                                                                                                                                                  <  4:1072548/149‑1 (MQ=255)
  aTCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaa                                                                                                                                                  <  1:907188/149‑1 (MQ=255)
   tCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaac                                                                                                                                                 <  2:538666/149‑1 (MQ=255)
          aCCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAag                                                                                                                                          <  1:1558574/149‑1 (MQ=255)
          aCCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTAATAACTAACAag                                                                                                                                          >  2:1558574/1‑149 (MQ=255)
              tCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagatag                                                                                                                                      >  2:723793/1‑149 (MQ=255)
                gggAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagatagat                                                                                                                                    <  1:1135147/149‑1 (MQ=255)
                   agcTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAGGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGtt                                                                                                                                 <  4:1624886/147‑1 (MQ=255)
                   aaCTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGtt                                                                                                                                 >  1:1301293/1‑149 (MQ=255)
                   aaCTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGtt                                                                                                                                 <  2:1301293/149‑1 (MQ=255)
                         aTCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaa                                                                                                                           >  2:415534/1‑149 (MQ=255)
                         aTCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaa                                                                                                                           <  4:1610867/149‑1 (MQ=255)
                           ccccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaaa                                                                                                                         <  1:1341553/149‑1 (MQ=255)
                            cccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                    <  3:405322/121‑1 (MQ=255)
                            cccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                    >  4:405322/1‑121 (MQ=255)
                              cTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACAt                                                                                                                      >  4:1184628/1‑149 (MQ=255)
                               tttCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                               <  3:396623/139‑1 (MQ=255)
                               tttCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                               >  4:396623/1‑139 (MQ=255)
                               tttCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                     <  1:462165/149‑1 (MQ=255)
                                    aGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCaa                                                                                                                <  3:1184628/149‑1 (MQ=255)
                                    aGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCaa                                                                                                                >  3:236276/1‑149 (MQ=255)
                                         taAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                           >  4:454454/1‑149 (MQ=255)
                                          aaaGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                               <  2:591240/128‑1 (MQ=255)
                                          aaaGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                               >  1:591240/1‑128 (MQ=255)
                                            aGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                       <  3:1892080/118‑1 (MQ=255)
                                            aGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                       >  4:1892080/1‑118 (MQ=255)
                                                 aGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCAt                                                                                                   <  4:236276/149‑1 (MQ=255)
                                                    cAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                       <  2:1388142/110‑1 (MQ=255)
                                                    cAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                       >  1:1388142/1‑110 (MQ=255)
                                                    cAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                        <  4:1208921/109‑1 (MQ=255)
                                                    cAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                        >  3:1208921/1‑109 (MQ=255)
                                                     aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTAtgac                                                                                                                                                 <  2:180930/99‑1 (MQ=255)
                                                     aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTAtgac                                                                                                                                                 >  1:180930/1‑99 (MQ=255)
                                                      aaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGt                                                                                                                   <  3:494210/128‑1 (MQ=255)
                                                      aaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGt                                                                                                                   >  4:494210/1‑128 (MQ=255)
                                                           aTACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                                                                         <  1:415534/149‑1 (MQ=255)
                                                             aCCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATg                                                                                       <  3:454454/149‑1 (MQ=255)
                                                               ctgctgCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTaa                                                                                            >  4:1539537/1‑142 (MQ=255)
                                                               ctgctgCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTaa                                                                                            <  3:1539537/142‑1 (MQ=255)
                                                                tgctgcTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGaa                                                                                    >  4:615871/1‑149 (MQ=255)
                                                                  ctgctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAAc                                                                                  <  3:1556909/149‑1 (MQ=255)
                                                                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                           <  1:1710928/90‑1 (MQ=255)
                                                                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                           >  2:1710928/1‑90 (MQ=255)
                                                                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGt                                                                                                                   <  1:1178430/114‑1 (MQ=255)
                                                                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGt                                                                                                                   >  2:1178430/1‑114 (MQ=255)
                                                                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaa                                                                                >  3:1406370/1‑149 (MQ=255)
                                                                     ctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaa                                                                                                                                              >  2:321585/1‑86 (MQ=255)
                                                                     ctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaa                                                                                                                                              <  1:321585/86‑1 (MQ=255)
                                                                      ttGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                           <  2:938314/120‑1 (MQ=255)
                                                                      ttGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                           >  1:938314/1‑120 (MQ=255)
                                                                        gtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaaa                                                                            <  4:1406370/149‑1 (MQ=255)
                                                                             aGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       <  1:723793/149‑1 (MQ=255)
                                                                              gttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAAc                                                                      <  3:615871/149‑1 (MQ=255)
                                                                               ttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                     <  1:1852806/101‑1 (MQ=255)
                                                                               ttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                     >  2:1852806/1‑101 (MQ=255)
                                                                                tttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTa                                                                    >  1:36350/1‑149 (MQ=255)
                                                                                    tttttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTaa                                                                                            >  3:426844/1‑121 (MQ=255)
                                                                                    tttttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTaa                                                                                            <  4:426844/121‑1 (MQ=255)
                                                                                          aaCCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCg                                                                 <  3:1757592/142‑1 (MQ=255)
                                                                                          aaCCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCg                                                                 >  4:1757592/1‑142 (MQ=255)
                                                                                            ccGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCa                                                        >  2:1256582/1‑149 (MQ=255)
                                                                                                  tACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAg                                                                                                >  4:334142/1‑103 (MQ=255)
                                                                                                  tACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAg                                                                                                <  3:334142/103‑1 (MQ=255)
                                                                                                     aaaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtc                                               >  3:784265/1‑149 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       <  4:1448435/124‑1 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       <  4:810417/124‑1 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       <  4:827794/124‑1 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       >  3:810417/1‑124 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       >  3:827794/1‑124 (MQ=255)
                                                                                                      aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       >  3:1448435/1‑124 (MQ=255)
                                                                                                             gTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATACCTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAAc                                                                                           >  3:1232380/1‑97 (MQ=255)
                                                                                                             gTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAAc                                                                                           <  4:1232380/97‑1 (MQ=255)
                                                                                                              tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                           <  1:1163998/80‑1 (MQ=255)
                                                                                                              tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                           >  2:1163998/1‑80 (MQ=255)
                                                                                                              tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaaa                                                                            >  1:542126/1‑111 (MQ=255)
                                                                                                              tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaaa                                                                            <  2:542126/111‑1 (MQ=255)
                                                                                                                 gTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       >  2:1440224/1‑113 (MQ=255)
                                                                                                                 gTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTaa                                                                       <  1:1440224/113‑1 (MQ=255)
                                                                                                                  ttttAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCt                                  >  4:732984/1‑149 (MQ=255)
                                                                                                                  ttttAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCt                                  <  1:1256582/149‑1 (MQ=255)
                                                                                                                      aaGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaaa                              <  2:36350/149‑1 (MQ=255)
                                                                                                                       aGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                     <  2:628322/61‑1 (MQ=255)
                                                                                                                       aGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                     >  1:628322/1‑61 (MQ=255)
                                                                                                                        gCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaaa                              >  3:1301403/1‑147 (MQ=255)
                                                                                                                        gCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaaa                              <  4:1301403/147‑1 (MQ=255)
                                                                                                                                     cTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCAt                                                       <  4:709089/109‑1 (MQ=255)
                                                                                                                                     cTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCAt                                                       >  3:709089/1‑109 (MQ=255)
                                                                                                                                      tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGtctc              >  2:1121402/1‑149 (MQ=255)
                                                                                                                                       ttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGGctct             <  4:784265/149‑1 (MQ=255)
                                                                                                                                        tcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCtt            >  4:1283542/3‑149 (MQ=255)
                                                                                                                                          ttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATc                                                                  <  2:1133257/93‑1 (MQ=255)
                                                                                                                                          ttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATc                                                                  >  1:1133257/1‑93 (MQ=255)
                                                                                                                                             ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTaa       >  4:323954/1‑149 (MQ=255)
                                                                                                                                                 tCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGcc   >  4:1403091/1‑149 (MQ=255)
                                                                                                                                                  cATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGccc  <  3:732984/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCC  >  NC_000913/1647727‑1648021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: