Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,613 A→C 5.2% Y307D (TAC→GAC) ‡ iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,660,6130AC 5.2% 715.3 / 23.8 193Y307D (TAC→GAC) ‡iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (94/89);  new base C (3/7);  total (97/96)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 200 total aligned reads displayed.

GAGAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGAT  >  NC_000913/2660465‑2660760
                                                                                                                                                    |                                                                                                                                                   
  gaAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGt                                                                                                                                                   <  1:263046/149‑1 (MQ=255)
      cGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTa                                                                                                                                                     >  4:964757/1‑143 (MQ=255)
                tATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAg                                                                                                                                              <  2:204385/140‑1 (MQ=255)
                tATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAg                                                                                                                                              >  1:204385/1‑140 (MQ=255)
                     gCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAAt                                                                                                                                <  3:439738/149‑1 (MQ=255)
                      ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAggt                                                                                                                                            >  4:64718/1‑134 (MQ=255)
                      ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGtt                                                                                                                                                  >  1:1525832/1‑130 (MQ=255)
                      ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTga                                                                                                                                    >  3:864893/1‑144 (MQ=255)
                      ccAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATg                                                                                                                               >  2:1053288/1‑149 (MQ=255)
                         gCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTc                                                                                                                            <  4:1051423/149‑1 (MQ=255)
                           tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGtt                                                                                                                                                  <  3:1163010/125‑1 (MQ=255)
                           tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGtt                                                                                                                                                  >  4:1163010/1‑125 (MQ=255)
                           tcgtcgTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCgg                                                                                                                          <  4:23602/149‑1 (MQ=255)
                            cgtcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgc                                                                                                                         >  1:1705229/1‑149 (MQ=255)
                             gtcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAg                                                                                                                                              <  4:1856780/127‑1 (MQ=255)
                             gtcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAg                                                                                                                                              >  3:1856780/1‑127 (MQ=255)
                              tcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACTACGTTGAAg                                                                                                                                              <  3:123670/126‑1 (MQ=255)
                              tcgtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGgcgc                                                                                                                       >  2:1431978/1‑149 (MQ=255)
                                gtTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCAc                                                                                                                     <  2:1705229/149‑1 (MQ=255)
                                 ttCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCAcc                                                                                                                    >  1:1447771/1‑149 (MQ=255)
                                        ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCa                                                                                                                      <  1:1786139/140‑1 (MQ=255)
                                        ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCa                                                                                                                      >  2:1786139/1‑140 (MQ=255)
                                        ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTcc                                                                                                             >  2:83342/1‑149 (MQ=255)
                                        ccAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTcc                                                                                                             >  4:859871/1‑149 (MQ=255)
                                                          aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCAc                                                                                                      <  1:733166/138‑1 (MQ=255)
                                                          aGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCAc                                                                                                      >  2:733166/1‑138 (MQ=255)
                                                             aCGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAAct                                                                                        >  3:1713220/1‑149 (MQ=255)
                                                                 ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAg                                                                                                                                                    <  2:1010894/85‑1 (MQ=255)
                                                                 ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAg                                                                                                                                                    >  1:1010894/1‑85 (MQ=255)
                                                                 ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGt                                                                                                                                       >  1:902195/1‑98 (MQ=255)
                                                                 ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGt                                                                                                                                       <  2:902195/98‑1 (MQ=255)
                                                                 ttCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAActtctt                                                                                    <  1:134921/149‑1 (MQ=255)
                                                                  tCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTaaa                                                                                          <  1:1431348/142‑1 (MQ=255)
                                                                  tCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                   <  1:83342/149‑1 (MQ=255)
                                                                     aGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAg                                                                                                           >  4:344742/1‑122 (MQ=255)
                                                                     aGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAg                                                                                                           <  3:344742/122‑1 (MQ=255)
                                                                         ttGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                            <  1:131871/149‑1 (MQ=255)
                                                                         ttGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                            >  3:653728/1‑149 (MQ=255)
                                                                          tGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCtt                                                                           <  2:1477359/149‑1 (MQ=255)
                                                                           gCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          <  2:1121117/149‑1 (MQ=255)
                                                                           gCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          <  4:588918/149‑1 (MQ=255)
                                                                           gCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          >  1:1121117/1‑149 (MQ=255)
                                                                            cTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGa                                                                                 >  4:188767/1‑141 (MQ=255)
                                                                            cTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGa                                                                                 <  3:188767/141‑1 (MQ=255)
                                                                              gACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAg                                                                                               >  4:385615/1‑125 (MQ=255)
                                                                                            aCCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGAc                                                         <  3:1309237/149‑1 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAg                                                                                                                                                    <  3:1815840/56‑1 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAg                                                                                                                                                    >  4:1815840/1‑56 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTc                                                                                                        <  3:1069913/100‑1 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTc                                                                                                        >  4:1069913/1‑100 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCAcc                                                                                                     <  3:176964/103‑1 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCAcc                                                                                                     >  4:176964/1‑103 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGc                                                                     >  2:519924/1‑135 (MQ=255)
                                                                                              cTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGc                                                                     <  1:519924/135‑1 (MQ=255)
                                                                                                 aagaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        <  3:603136/129‑1 (MQ=255)
                                                                                                  agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCa                                                                                                            >  3:241396/1‑92 (MQ=255)
                                                                                                  agaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAg                                                                                               <  4:426312/105‑1 (MQ=255)
                                                                                                   gaaACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACg                                                  >  1:636688/1‑149 (MQ=255)
                                                                                                    aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGc                                                 >  3:67566/1‑149 (MQ=255)
                                                                                                    aaaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGc                                                 >  4:1281976/1‑149 (MQ=255)
                                                                                                     aaCTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCa                                                <  3:372338/149‑1 (MQ=255)
                                                                                                       cTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAgg                                              <  3:1539489/149‑1 (MQ=255)
                                                                                                            aGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCcgcg                                         >  4:759082/1‑149 (MQ=255)
                                                                                                              gTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                   <  1:1194922/105‑1 (MQ=255)
                                                                                                                  ttCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACg                                   <  4:67566/149‑1 (MQ=255)
                                                                                                                     aGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTc                                <  3:759082/149‑1 (MQ=255)
                                                                                                                      gcgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          >  1:284511/1‑106 (MQ=255)
                                                                                                                      gcgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTcc                               <  3:1437743/149‑1 (MQ=255)
                                                                                                                       cgcCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          >  3:392981/1‑105 (MQ=255)
                                                                                                                         ccATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        >  2:1422599/1‑105 (MQ=255)
                                                                                                                            tAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCg                         <  4:1713220/149‑1 (MQ=255)
                                                                                                                               tCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTc                      >  4:152611/1‑149 (MQ=255)
                                                                                                                                 aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTc                                                                                                 >  2:843287/1‑72 (MQ=255)
                                                                                                                                 aGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGAAGGtt                                                     >  1:1301870/1‑116 (MQ=255)
                                                                                                                                  gCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        <  2:707217/96‑1 (MQ=255)
                                                                                                                                  gCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        >  1:707217/1‑96 (MQ=255)
                                                                                                                                     cgtcagcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTc                                <  2:265895/127‑1 (MQ=255)
                                                                                                                                     aCTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtc                >  1:1087330/1‑149 (MQ=255)
                                                                                                                                       tCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCGGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGa                                     <  2:670103/126‑1 (MQ=255)
                                                                                                                                       tCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCGGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGa                                     >  1:670103/1‑126 (MQ=255)
                                                                                                                                       tCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTc                                                                <  2:1338815/99‑1 (MQ=255)
                                                                                                                                       tCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtctc              >  4:1586202/1‑149 (MQ=255)
                                                                                                                                         agcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGt                                                                                             <  1:566086/66‑1 (MQ=37)
                                                                                                                                         agcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGt                                                                                             >  2:566086/3‑68 (MQ=37)
                                                                                                                                         agcTTCAACTACGTTGAAGGTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCAtct               <  1:864602/144‑1 (MQ=255)
                                                                                                                                         aCCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                   <  1:1650521/78‑1 (MQ=255)
                                                                                                                                         aCCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACttcttc                                                                                   >  2:1650521/1‑78 (MQ=255)
                                                                                                                                           cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          <  1:1640612/85‑1 (MQ=255)
                                                                                                                                           cTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCttt                                                                          <  3:852261/85‑1 (MQ=255)
                                                                                                                                            ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCt                                                                            <  2:391216/82‑1 (MQ=255)
                                                                                                                                            ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        <  2:468247/86‑1 (MQ=255)
                                                                                                                                            ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGa                                                                        >  1:468247/1‑86 (MQ=255)
                                                                                                                                            ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt         >  2:1852131/1‑149 (MQ=255)
                                                                                                                                            ttCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTcttctt         >  4:1439324/1‑149 (MQ=255)
                                                                                                                                              cAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCtttt       <  2:228281/149‑1 (MQ=255)
                                                                                                                                                aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgc     >  1:173881/1‑149 (MQ=255)
                                                                                                                                                aCGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtg      >  4:1344125/1‑148 (MQ=255)
                                                                                                                                                 cGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcg    <  2:1567544/149‑1 (MQ=255)
                                                                                                                                                   tAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgat  <  2:636688/149‑1 (MQ=255)
                                                                                                                                                   tAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTtgcgat  <  2:1226129/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGAAACGGATAGAGCTATGTGCCAGCTCGTCGTTCAGCCCCAGCGCGCGCAGCACGTAGGACGGTTCGAGGCTTGCTGACGTACAGGCGGAACCTGAAGAAACTGCGAGGTCTTTCAGCGCCATAATCAGCGACTCACCTTCAACGTAGTTGAAGCTGACGTTGAGAATGTTCGGCGCACCGTGTTCCAGGTCACCGTTCAGGTAAACTTCTTCGATATCTTTGATGCCGTTCCACAGACGGTTACGCAGGCCGCGCAGACGTTCCATCTCGGTCGCCATCTCTTCTTTTGCGAT  >  NC_000913/2660465‑2660760

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: