Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,116,854 A→C 5.3% K16N (AAA→AAC yceA → putative rhodanese‑related sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,116,8540AC 5.3% 321.1 / 6.9 89K16N (AAA→AACyceAputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (50/34);  new base C (3/2);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCAT  >  NC_000913/1116707‑1117000
                                                                                                                                                   |                                                                                                                                                  
aaTGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCaaaa                                                                                                                                                   <  3:151881/149‑1 (MQ=255)
     gcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTg                                                                                                                                              <  4:118252/149‑1 (MQ=255)
       aaaTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTgg                                                                                                                                             <  2:330398/148‑1 (MQ=255)
         aTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCtt                                                                                                                                          <  4:620491/149‑2 (MQ=255)
                 tAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcga                                                                                                                                   >  1:161836/1‑148 (MQ=255)
                        tgtgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccgcgaa                                                                                                                           <  3:81543/149‑1 (MQ=255)
                           ggATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaaccgcgaacc                                                                                                                         >  2:64029/1‑148 (MQ=255)
                            gATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAacca                                                                                                                        <  1:909208/148‑1 (MQ=255)
                             aTGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAaccac                                                                                                                       >  1:186858/1‑148 (MQ=255)
                                ttAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCAt                                                                                                                    >  1:571518/1‑148 (MQ=255)
                                  aaTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTc                                                                                                                 >  4:803356/1‑149 (MQ=255)
                                    ttATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGt                                                                                                               >  3:85660/1‑149 (MQ=255)
                                     tATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGt                                                                                                               >  2:242550/1‑148 (MQ=255)
                                          aTAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCtctta                                                                                                                                       <  4:481602/119‑5 (MQ=255)
                                          aTAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACa                                                                                                         >  3:831994/1‑149 (MQ=255)
                                          aTAATTGCTAGATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCtttta                                                                                                                                       >  3:481602/1‑115 (MQ=255)
                                             aTTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGt                                                                                                      <  4:39018/149‑1 (MQ=255)
                                               tGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTa                                                                                                     >  1:752227/1‑148 (MQ=255)
                                                  tataTCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                  >  2:535734/1‑148 (MQ=255)
                                                  tataTCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                  <  2:204998/148‑1 (MQ=255)
                                                  tataTCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTc                                                                                                 <  3:803356/149‑1 (MQ=255)
                                                   atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTc                                                                                                 >  1:463733/1‑148 (MQ=255)
                                                      tCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCcaca                                                                                             >  3:767997/1‑149 (MQ=255)
                                                       cATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCcaca                                                                                             >  1:274171/1‑148 (MQ=255)
                                                           ccGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGc                                                                                         >  1:69745/1‑148 (MQ=255)
                                                             gcgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCg                                                                                       >  2:679821/1‑148 (MQ=255)
                                                             gcgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGa                                                                                      >  3:768206/1‑149 (MQ=255)
                                                             gcgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGa                                                                                      >  4:685370/1‑149 (MQ=255)
                                                              cgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGa                                                                                      >  2:583187/1‑148 (MQ=255)
                                                              cgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGa                                                                                      >  2:712883/1‑148 (MQ=255)
                                                              cgGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGa                                                                                      >  1:706464/1‑148 (MQ=255)
                                                               ggATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGAt                                                                                     >  2:322142/1‑148 (MQ=255)
                                                                    tttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                <  2:512652/148‑1 (MQ=255)
                                                                      tACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGc                                                                             >  3:275579/1‑149 (MQ=255)
                                                                         tttCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGAc                                                                          >  4:527837/1‑149 (MQ=255)
                                                                           tCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGAcc                                                                         >  2:418586/1‑148 (MQ=255)
                                                                            cccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                  <  1:142031/122‑1 (MQ=255)
                                                                            cccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTAttt                                                                                                  >  2:142031/1‑122 (MQ=255)
                                                                             ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                >  4:539640/1‑139 (MQ=255)
                                                                             ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                <  3:539640/139‑1 (MQ=255)
                                                                             ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGt                                                                      >  4:683677/1‑149 (MQ=255)
                                                                               aTCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTg                                                                     >  1:572178/1‑148 (MQ=255)
                                                                                  tcGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGc                                                                 >  4:412298/1‑149 (MQ=255)
                                                                                    gCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCt                                                                <  1:679821/148‑1 (MQ=255)
                                                                                    gCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCt                                                                <  2:212503/148‑1 (MQ=255)
                                                                                         aaCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtata                                                           >  2:497274/1‑148 (MQ=255)
                                                                                         aaCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtata                                                           >  2:518681/1‑148 (MQ=255)
                                                                                         aaCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTtata                                                           <  1:535734/148‑1 (MQ=255)
                                                                                          aCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATc                                                         <  3:374408/149‑1 (MQ=255)
                                                                                           ccGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCa                                                        >  4:158926/1‑149 (MQ=255)
                                                                                             gTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAg                                                       <  2:706464/148‑1 (MQ=255)
                                                                                              tACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGc                                                      <  1:242550/148‑1 (MQ=255)
                                                                                                     tGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTAc                                               <  1:322142/148‑1 (MQ=255)
                                                                                                      gCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcg                                             >  4:807101/1‑149 (MQ=255)
                                                                                                         aGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACcgcg                                           >  2:552964/1‑148 (MQ=255)
                                                                                                           tgtTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTg                                        >  3:66772/1‑149 (MQ=255)
                                                                                                             ttACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGaa                                      >  3:13737/1‑149 (MQ=255)
                                                                                                             ttACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGa                                       <  2:571518/148‑1 (MQ=255)
                                                                                                              tACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGaa                                      <  1:609966/148‑1 (MQ=255)
                                                                                                              tACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAAt                                     >  4:330338/1‑149 (MQ=255)
                                                                                                               acacAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                <  4:550833/105‑1 (MQ=255)
                                                                                                               acacAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTaa                                                                                >  3:550833/1‑105 (MQ=255)
                                                                                                                 acaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                                  >  4:280512/1‑149 (MQ=255)
                                                                                                                  caACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                                  <  2:463733/148‑1 (MQ=255)
                                                                                                                  caACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                                  >  1:524598/1‑148 (MQ=255)
                                                                                                                    aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGttttt                               <  4:541608/149‑1 (MQ=255)
                                                                                                                    aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtttt                                >  2:23298/1‑148 (MQ=255)
                                                                                                                    aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtttt                                <  1:23298/148‑1 (MQ=255)
                                                                                                                    aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtttt                                <  1:712883/148‑1 (MQ=255)
                                                                                                                     ccGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                              <  3:280512/149‑1 (MQ=255)
                                                                                                                      cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTgg                             <  4:275579/149‑1 (MQ=255)
                                                                                                                      cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                              >  1:801778/1‑148 (MQ=255)
                                                                                                                              cAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                               <  2:234418/91‑1 (MQ=255)
                                                                                                                              cAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                               >  1:234418/1‑91 (MQ=255)
                                                                                                                               aaCGACGCGCTAAAAGCCAAAATGTTGGGTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTa                    >  3:777945/1‑149 (MQ=255)
                                                                                                                               aaCGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgtgt                     >  2:287143/1‑148 (MQ=255)
                                                                                                                               aaCGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgtgt                     <  2:574725/148‑1 (MQ=255)
                                                                                                                                acgacgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTa                    <  2:274171/148‑1 (MQ=255)
                                                                                                                                acgacgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTAt                   <  4:13737/149‑1 (MQ=255)
                                                                                                                                     gcgcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTgg               >  2:278948/1‑148 (MQ=255)
                                                                                                                                       gcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                               >  3:135024/1‑82 (MQ=255)
                                                                                                                                       gcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                               <  4:135024/82‑1 (MQ=255)
                                                                                                                                         tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCAt          >  4:108300/1‑149 (MQ=255)
                                                                                                                                          ctcaGCCAACATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTTTCTGACGCAt          >  2:439635/4‑148 (MQ=255)
                                                                                                                                           tcaGCCAACATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGTATGTTTTTGGGCGAGTGTATCTGGCGCATGa        >  3:434089/3‑149 (MQ=255)
                                                                                                                                             aGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAgg      <  4:768206/149‑1 (MQ=255)
                                                                                                                                             aGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAg       <  1:855776/148‑1 (MQ=255)
                                                                                                                                              gCCAAAATGTTGGCTGAGAGCGAACCGCGGACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAgg      >  1:671991/1‑148 (MQ=255)
                                                                                                                                              gCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCtt                                                               >  1:687403/1‑91 (MQ=255)
                                                                                                                                              gCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCtt                                                               <  2:687403/91‑1 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                            <  2:449218/123‑1 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                            >  1:449218/1‑123 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCa   >  1:418423/1‑148 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCa   >  2:780790/1‑148 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCAt  >  4:79222/1‑149 (MQ=255)
                                                                                                                                                 aaaaTGTTGGCTGAGAGCGAACCGAGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGGGACCCGTGACGCTTTATATCAGCTGGTGACCGCGCTGAATGTTTTTCGGCGAGTGTATCTGGCGCATGAGGGCAt  >  3:765304/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATGCGCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCAT  >  NC_000913/1116707‑1117000

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: