Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,664,593 G→T 7.1% A30S (GCC→TCC)  dmsD → redox enzyme maturation protein DmsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,664,5930GT 7.1% 280.4 / 7.0 84A30S (GCC→TCC) dmsDredox enzyme maturation protein DmsD
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base T (2/4);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCG  >  NC_000913/1664454‑1664742
                                                                                                                                           |                                                                                                                                                      
tGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACctct                                                                                                                                             >  1:1100895/1‑150 (MQ=255)
   aGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGGTGCGGCTTCCGCACTCTCTgga                                                                                                                                          <  1:932615/150‑2 (MQ=255)
   aGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGGTGCGGCTTCCGCACTCTCTgga                                                                                                                                          <  1:532640/150‑2 (MQ=255)
     gTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTg                                                                                                                                        <  1:186397/150‑1 (MQ=255)
     attAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGAGCGCGTGTTGGGCGCGCTGTTTTATTAAGCTCAAGAGAGTGCGGAAGCCGCACATCTGGTTg                                                                                                                                        >  1:326068/2‑150 (MQ=255)
      ttAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGc                                                                                                                                       <  2:922810/150‑1 (MQ=255)
      ttAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGc                                                                                                                                       <  1:201285/150‑1 (MQ=255)
      ttAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGGTGCGGCTTCCGCACTCTCTggagcg                                                                                                                                       <  1:1007118/150‑5 (MQ=255)
           aaaTTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTAc                                                                                                                                  <  1:1321016/150‑1 (MQ=255)
            aaTTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACt                                                                                                                                 >  1:155271/1‑150 (MQ=255)
            aaTTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACt                                                                                                                                 >  2:362091/1‑150 (MQ=255)
             aTTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGAGCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTc                                                                                                                                <  2:53399/150‑1 (MQ=255)
                ggTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCGGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCAcc                                                                                                                             <  2:870212/150‑1 (MQ=255)
                     ggggCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGa                                                                                                                        <  1:321849/150‑1 (MQ=255)
                        gCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACgg                                                                                                                     >  2:552002/1‑150 (MQ=255)
                           aCCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTg                                                                                                                  >  1:740465/1‑150 (MQ=255)
                             ccGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTggg                                                                                                                <  1:1325997/150‑1 (MQ=255)
                                cacaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGaaa                                                                                                             >  1:744327/1‑150 (MQ=255)
                                 acaTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAAc                                                                                                            <  1:691546/150‑1 (MQ=255)
                                  caTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACt                                                                                                           >  1:580226/1‑150 (MQ=255)
                                  caTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGGTGCGGCTTCCGCACTCTCTgga                                                                                                                                          <  1:777379/119‑2 (MQ=38)
                                  caTTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGGTGCGGCTTCCGCACTCTCTgga                                                                                                                                          >  2:777379/1‑118 (MQ=38)
                                    tttCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGAGCGCTGTTTTATTACGCTACAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCa                                                                                                         >  2:1177474/1‑150 (MQ=255)
                                     ttCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAg                                                                                                        <  2:580226/150‑1 (MQ=255)
                                        aGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTgg                                                                                                     >  1:725645/1‑150 (MQ=255)
                                             gTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCctct                                                                                                >  1:952932/1‑150 (MQ=255)
                                               tAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAc                                                                                              <  1:1419276/150‑1 (MQ=255)
                                                aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCCGTGGCCTCTAcc                                                                                             >  2:1592283/1‑150 (MQ=255)
                                                aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAcc                                                                                             <  1:946402/150‑1 (MQ=255)
                                                aGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTAcc                                                                                             <  1:552002/150‑1 (MQ=255)
                                                 ggAATCACCAATTTTTCACAGCAAGATAATATTACTGTCGCGGAGCGCGTGATGGGAGAGATGTATTATAACACTCCAGAGAGTGCGGAAGCCGCACCTCTCGTTGCGGTACTCACAAGTGACGGCTGGGAAACTCAGTGGCCTCTACCa                                                                                            >  2:1452495/1‑150 (MQ=255)
                                                  gAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAg                                                                                           >  2:661166/1‑150 (MQ=255)
                                                   aaTGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGa                                                                                          <  1:269017/150‑1 (MQ=255)
                                                   aaTGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGa                                                                                          <  1:669955/150‑1 (MQ=255)
                                                     tGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAg                                                                                        <  1:87224/150‑1 (MQ=255)
                                                       aCCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCg                                                                                      <  1:601293/150‑1 (MQ=255)
                                                        cccATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGt                                                                                     >  2:1509664/1‑150 (MQ=255)
                                                             ttttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTa                                                                                >  2:654820/1‑150 (MQ=255)
                                                               ttCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAgc                                                                              >  2:1346829/1‑150 (MQ=255)
                                                                   caGCATGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGcgccg                                                                          >  1:1514202/1‑150 (MQ=255)
                                                                   caGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGcgccg                                                                          >  1:809017/1‑150 (MQ=255)
                                                                   caGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGcgccg                                                                          >  2:1435961/1‑150 (MQ=255)
                                                                       aaGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTgg                                                                      <  2:744327/150‑1 (MQ=255)
                                                                           tAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGAc                                                                  >  1:407909/1‑150 (MQ=255)
                                                                              tttttCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCACCAGCGGCGCTAATGACgct                                                               <  2:559644/150‑3 (MQ=255)
                                                                                  tCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCtttt                                                           <  1:51575/150‑1 (MQ=255)
                                                                                  tCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCACCAGCGGCGCTAATGACgcttctg                                                           <  1:31981/150‑7 (MQ=255)
                                                                                   cTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTc                                                          <  2:228492/150‑1 (MQ=255)
                                                                                    tGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCa                                                         <  2:1100895/150‑1 (MQ=255)
                                                                                    tGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCa                                                         <  2:1406336/150‑1 (MQ=255)
                                                                                      tCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGa                                                       <  2:725645/150‑1 (MQ=255)
                                                                                         cggcgCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGacac                                                    <  2:205351/150‑1 (MQ=255)
                                                                                          ggcgCGCGTGTTGGGCGCGCTGTTTTATTACGCTCGAGAGAGTGCTGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACTTCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGAcaca                                                   <  2:326068/150‑1 (MQ=255)
                                                                                                 gtgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGaa                                            <  1:362091/150‑1 (MQ=255)
                                                                                                 gtgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGaa                                            <  1:446616/150‑1 (MQ=255)
                                                                                                  tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAag                                           <  2:740465/150‑1 (MQ=255)
                                                                                                  tgtTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAag                                           <  1:654820/150‑1 (MQ=255)
                                                                                                    ttGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAagag                                         <  2:973152/150‑1 (MQ=255)
                                                                                                      gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGAc                                       >  2:663976/1‑150 (MQ=255)
                                                                                                      gggCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGAc                                       >  1:916111/1‑150 (MQ=255)
                                                                                                        gcgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTc                                     >  2:312314/1‑150 (MQ=255)
                                                                                                         cgcgcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCa                                    <  2:153338/150‑1 (MQ=255)
                                                                                                            gcTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGc                                 >  1:289151/1‑150 (MQ=255)
                                                                                                             cTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGCCTCACGcc                                >  1:1572364/1‑150 (MQ=255)
                                                                                                              tGTATTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGccc                               >  2:250498/1‑150 (MQ=255)
                                                                                                                    attaCGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCt                         >  1:750475/1‑150 (MQ=255)
                                                                                                                      taCGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTgg                       >  1:305860/1‑150 (MQ=255)
                                                                                                                       aCGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGc                      <  1:719197/150‑1 (MQ=255)
                                                                                                                       aCGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGc                      >  2:719197/1‑150 (MQ=255)
                                                                                                                         gCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGTCCTGGCAg                    >  1:839108/1‑150 (MQ=255)
                                                                                                                         gCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACAAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAg                    >  2:555946/1‑150 (MQ=255)
                                                                                                                            ccAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgt                 <  1:296064/150‑1 (MQ=255)
                                                                                                                            ccAGAGGTGCGGCTTCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtt                >  2:813558/1‑150 (MQ=255)
                                                                                                                             cAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGcc                          <  1:930004/140‑1 (MQ=255)
                                                                                                                             cAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGcc                          >  2:930004/1‑140 (MQ=255)
                                                                                                                               gagagTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttg              >  1:1338068/1‑150 (MQ=255)
                                                                                                                               gagagTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCCGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttg              >  2:879069/1‑150 (MQ=255)
                                                                                                                                 gagTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt            <  1:261934/150‑1 (MQ=255)
                                                                                                                                 gagTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCgtttgtt            <  1:1177474/150‑1 (MQ=255)
                                                                                                                                   gTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGtttgtttg          >  1:409737/1‑150 (MQ=255)
                                                                                                                                     gCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTc        <  2:155271/150‑1 (MQ=255)
                                                                                                                                      cGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCACCAGCGGCGCTAATGACgcttc                                                             <  2:540190/96‑5 (MQ=38)
                                                                                                                                      cGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCACCAGCGGCGCTAATGACgcttc                                                             >  1:540190/1‑92 (MQ=38)
                                                                                                                                         aaGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGcc    >  1:1406789/1‑150 (MQ=255)
                                                                                                                                           gCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCg  >  2:692597/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                      
TGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGTTGGGCGCGCTGTTTTATTACGCTCCAGAGAGTGCGGAAGCCGCACCTCTGGTTGCGGTACTCACCAGTGACGGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCG  >  NC_000913/1664454‑1664742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: