Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,338 G→C 10.1% intergenic (+53/‑177) nadE → / → cho NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,823,3380GC10.1% 471.3 / 35.7 139intergenic (+53/‑177)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (59/66);  new base C (9/5);  total (68/71)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 145 total aligned reads displayed.

AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTG  >  NC_000913/1823191‑1823486
                                                                                                                                                   |                                                                                                                                                    
aaCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgc                                                                                                                                                     >  2:1124240/1‑149 (MQ=255)
     gTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGt                                                                                                                                                <  4:1254378/149‑1 (MQ=255)
       cGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGa                                                                                                                                              <  4:560702/149‑1 (MQ=255)
         ccAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                                                                            <  3:510621/149‑1 (MQ=255)
           aGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                          >  3:846400/1‑149 (MQ=255)
                aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGTGGCGTGAACGCCTTAt                                                                                                                                     >  4:1210356/1‑149 (MQ=255)
                aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                     >  4:998336/1‑149 (MQ=255)
                aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                     >  2:1396419/1‑149 (MQ=255)
                 aTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                    <  4:1266873/133‑1 (MQ=255)
                  tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTc                                                                                                                                   <  3:931727/149‑1 (MQ=255)
                  tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTc                                                                                                                                   <  3:434573/149‑1 (MQ=255)
                    gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAGGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                 >  4:1339910/1‑149 (MQ=255)
                    gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                 <  2:1111032/149‑1 (MQ=255)
                      gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                               <  3:786286/149‑1 (MQ=255)
                             tATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaaa                                                                                                                        <  3:1287867/149‑1 (MQ=255)
                               tCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACt                                                                                                                      >  3:1614927/1‑149 (MQ=255)
                                cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                     <  1:1521278/149‑1 (MQ=255)
                                 tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                    >  4:1383584/1‑149 (MQ=255)
                                 tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                    >  2:1583246/1‑149 (MQ=255)
                                 tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                    <  2:1577334/149‑1 (MQ=255)
                                  gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                   >  4:24658/1‑149 (MQ=255)
                                   aaaaCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                  >  2:724562/1‑149 (MQ=255)
                                           aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                          <  1:566895/149‑1 (MQ=255)
                                             aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                        <  2:70849/149‑1 (MQ=255)
                                             aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                        >  1:777675/1‑149 (MQ=255)
                                                    ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCa                                                                                                 <  3:1489061/149‑1 (MQ=255)
                                                         cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                            <  1:1124240/149‑1 (MQ=255)
                                                             aaTTACCGTTTTCGATGATTTCTGGAAAGAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                        >  2:195202/1‑149 (MQ=255)
                                                             aaTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                        >  3:215344/1‑149 (MQ=255)
                                                             aaTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                        <  4:1136439/149‑1 (MQ=255)
                                                              aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                       <  3:1330092/149‑1 (MQ=255)
                                                              aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                       >  1:1366623/1‑149 (MQ=255)
                                                                 aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAgg                                                                                    <  1:14502/149‑1 (MQ=255)
                                                                   cGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGcc                                                                                  >  4:504888/1‑149 (MQ=255)
                                                                    gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                 <  4:79969/149‑1 (MQ=255)
                                                                    gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                              >  3:1280425/1‑88 (MQ=255)
                                                                      tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGa                                                                               <  1:1117909/149‑1 (MQ=255)
                                                                       ttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGAc                                                                              <  4:1279346/149‑1 (MQ=255)
                                                                           atgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                 >  3:1130724/1‑94 (MQ=255)
                                                                             gatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGt                                                                        <  2:1122295/149‑1 (MQ=255)
                                                                              atTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTa                                                                       >  3:1403196/1‑149 (MQ=255)
                                                                                 tCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATa                                                                                      <  4:889227/131‑1 (MQ=255)
                                                                                  cTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                   <  4:1255972/149‑1 (MQ=255)
                                                                                     gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc                                                                >  4:202627/1‑149 (MQ=255)
                                                                                       aaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg                                                              >  3:569457/1‑149 (MQ=255)
                                                                                            taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                         >  1:705575/1‑149 (MQ=255)
                                                                                            taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                         >  2:1254538/1‑149 (MQ=255)
                                                                                             aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                        >  2:538059/1‑149 (MQ=255)
                                                                                             aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                        <  3:928718/149‑1 (MQ=255)
                                                                                                  tttGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAtt                                                                                                             >  4:1119986/1‑91 (MQ=255)
                                                                                                     gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                >  1:1515043/1‑149 (MQ=255)
                                                                                                     gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                <  1:1254538/149‑1 (MQ=255)
                                                                                                     gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                <  2:756191/149‑1 (MQ=255)
                                                                                                          tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                         <  4:27875/87‑1 (MQ=255)
                                                                                                          tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                           >  3:1483206/1‑149 (MQ=255)
                                                                                                          tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                           >  2:369375/1‑149 (MQ=255)
                                                                                                           atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATc                                          <  4:1483206/149‑1 (MQ=255)
                                                                                                            ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                         >  1:1155457/1‑149 (MQ=255)
                                                                                                            ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                         >  3:998938/1‑149 (MQ=255)
                                                                                                            ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                         >  1:25642/1‑149 (MQ=255)
                                                                                                              ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                       <  3:504888/149‑1 (MQ=255)
                                                                                                              ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                       >  2:960803/1‑149 (MQ=255)
                                                                                                               gATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAcc                                      <  1:924065/149‑1 (MQ=255)
                                                                                                                aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                     <  2:767784/149‑1 (MQ=255)
                                                                                                                aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                     <  1:960803/149‑1 (MQ=255)
                                                                                                                 ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCg                                    >  3:1121011/1‑149 (MQ=255)
                                                                                                                  tgtgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGc                                   <  4:240630/149‑1 (MQ=255)
                                                                                                                   gtgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGcc                                  <  1:740820/149‑1 (MQ=255)
                                                                                                                    tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCct                                 <  1:1189945/149‑1 (MQ=255)
                                                                                                                        aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttt                             >  4:98548/1‑149 (MQ=255)
                                                                                                                             gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                     >  4:1349108/1‑136 (MQ=255)
                                                                                                                             gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                     <  3:1349108/136‑1 (MQ=255)
                                                                                                                              agTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAAt                       <  4:846400/149‑1 (MQ=255)
                                                                                                                               gTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                            <  4:1594422/127‑1 (MQ=255)
                                                                                                                               gTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                            >  3:1594422/1‑127 (MQ=255)
                                                                                                                                  aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                   >  4:553941/1‑149 (MQ=255)
                                                                                                                                   tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                   <  1:105657/100‑1 (MQ=255)
                                                                                                                                   tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCAc                                                                   >  2:105657/1‑100 (MQ=255)
                                                                                                                                   tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCtt                  >  1:157539/1‑149 (MQ=255)
                                                                                                                                   tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCtt                  >  4:388562/1‑149 (MQ=255)
                                                                                                                                     tGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTcc                <  2:25642/149‑1 (MQ=255)
                                                                                                                                      gTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTccc               <  1:222445/149‑1 (MQ=255)
                                                                                                                                      gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGTGTCGAATAGGCCTGACa                                                                             >  2:1335443/1‑87 (MQ=38)
                                                                                                                                      gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACa                                                                             <  1:1335443/87‑1 (MQ=255)
                                                                                                                                       ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg              >  4:802439/1‑149 (MQ=255)
                                                                                                                                       ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCCACCCGCCTCTTTTTTAATCGCATACCg              >  4:1089573/1‑149 (MQ=255)
                                                                                                                                       ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg              >  2:469672/1‑149 (MQ=255)
                                                                                                                                       ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAACTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg              >  1:65362/1‑149 (MQ=255)
                                                                                                                                       ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATCGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg              >  1:1201752/1‑149 (MQ=255)
                                                                                                                                         cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc            <  1:1583246/149‑1 (MQ=255)
                                                                                                                                         cACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc            >  3:374397/1‑149 (MQ=255)
                                                                                                                                          aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt           <  2:745270/149‑1 (MQ=255)
                                                                                                                                          aCGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt           >  2:941648/1‑149 (MQ=255)
                                                                                                                                           cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                         <  4:671144/102‑1 (MQ=255)
                                                                                                                                           cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                         >  3:671144/1‑102 (MQ=255)
                                                                                                                                            gCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                         >  4:1688154/1‑101 (MQ=255)
                                                                                                                                                gATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca     <  4:58570/149‑1 (MQ=255)
                                                                                                                                                  tGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACt   <  3:553941/149‑1 (MQ=255)
                                                                                                                                                   gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                        <  2:693086/63‑1 (MQ=255)
                                                                                                                                                   gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTg  <  3:263487/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTG  >  NC_000913/1823191‑1823486

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