Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,117 C→G 5.9% H126Q (CAC→CAG cbdA → cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,1170CG 5.9% 482.8 / 18.2 135H126Q (CAC→CAGcbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (67/60);  new base G (4/4);  total (71/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 141 total aligned reads displayed.

TTGGTACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCT  >  NC_000913/1037969‑1038265
                                                                                                                                                    |                                                                                                                                                    
ttGGTACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAc                                                                                                                                                      <  1:104158/149‑1 (MQ=255)
ttGGTACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAc                                                                                                                                                      <  2:224802/149‑1 (MQ=255)
    tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                  <  3:380839/149‑1 (MQ=255)
     aCTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgct                                                                                                                                                 <  1:402158/149‑1 (MQ=255)
          cTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGa                                                                                                                                            >  2:73591/1‑149 (MQ=255)
           tGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGAc                                                                                                                                           >  4:412451/1‑149 (MQ=255)
           tGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGAc                                                                                                                                           <  2:993073/149‑1 (MQ=255)
                aTTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCACCGGCTGAATAAATACCAGCACCTGCTGGTCACGTGGc                                                                                                                                      <  4:447215/149‑1 (MQ=255)
                 ttCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                                                                     >  4:288810/1‑149 (MQ=255)
                 ttCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                                                                     >  2:128656/1‑149 (MQ=255)
                  tCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtg                                                                                                                                    >  1:1014365/1‑149 (MQ=255)
                     aTTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCGACGGCTGAATAAATACCAGCCCCTGCTGGTGACGTGGCtggtg                                                                                                                                 <  1:349067/149‑1 (MQ=255)
                           aCTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTc                                                                                                                           <  4:1050900/149‑1 (MQ=255)
                               tgtgGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGtt                                                                                                                       <  2:381163/149‑1 (MQ=255)
                                 tgGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggt                                                                                                                                  >  2:1111223/1‑136 (MQ=255)
                                 tgGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggt                                                                                                                                  <  1:1111223/136‑1 (MQ=255)
                                 tgGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCa                                                                                                                     >  4:571696/1‑149 (MQ=255)
                                    gCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                                                                  <  2:489570/149‑1 (MQ=255)
                                    gCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                                                                  <  2:471459/149‑1 (MQ=255)
                                          tttttGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATctctct                                                                                                            <  4:775877/149‑1 (MQ=255)
                                             ttGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCg                                                                                                         <  1:930472/149‑1 (MQ=255)
                                                 cgcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                     <  2:771885/149‑1 (MQ=255)
                                                  gcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgtg                                                                                                    <  1:73591/149‑1 (MQ=255)
                                                   cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTgg                                                                                                   <  3:1129849/149‑1 (MQ=255)
                                                   cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTgg                                                                                                   >  3:907985/1‑149 (MQ=255)
                                                      cGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGata                                                                                                <  1:128656/149‑1 (MQ=255)
                                                      cGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTCGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGata                                                                                                <  1:728417/149‑1 (MQ=255)
                                                         tGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTg                                                                                             >  4:92076/1‑149 (MQ=255)
                                                           gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                              >  2:902233/1‑130 (MQ=255)
                                                           gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctct                                                                                                              <  1:902233/130‑1 (MQ=255)
                                                           gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                           <  4:1102638/149‑1 (MQ=255)
                                                           gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                           <  4:697056/149‑1 (MQ=255)
                                                             gATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATg                                                                                         >  3:123167/1‑149 (MQ=255)
                                                                 gAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCaa                                                                                     >  1:277848/1‑149 (MQ=255)
                                                                   aGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAc                                                                                                                                                      <  3:676376/82‑1 (MQ=255)
                                                                      aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                                >  2:114245/1‑149 (MQ=255)
                                                                      aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                                <  3:92076/149‑1 (MQ=255)
                                                                        tAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTgg                                                                              <  3:397270/149‑1 (MQ=255)
                                                                        tAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTgg                                                                              >  1:246093/1‑149 (MQ=255)
                                                                         aaTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAAc                                                                                    <  4:15126/142‑1 (MQ=255)
                                                                           tGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                             <  4:764655/83‑2 (MQ=255)
                                                                           tGGCCTTCTTCCTCGAATCCACCTTTGGCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                             >  3:764655/1‑82 (MQ=25)
                                                                            ggCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGc                                                                          <  3:571696/149‑1 (MQ=255)
                                                                             gCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCa                                                                         <  4:123167/149‑1 (MQ=255)
                                                                              ccTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgtg                                                                                                    >  1:483560/1‑121 (MQ=255)
                                                                               cttcttCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtat                                                                                                                                            >  4:893898/1‑78 (MQ=37)
                                                                               cttcttCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtat                                                                                                                                            <  3:893898/80‑3 (MQ=37)
                                                                               cttcttCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAAt                                                                       >  2:1013223/1‑149 (MQ=255)
                                                                                ttcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATa                                                                      <  2:579749/149‑1 (MQ=255)
                                                                                 tcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATAc                                                                     <  2:498205/149‑1 (MQ=255)
                                                                                       tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGatat                                                                                               <  3:654679/117‑1 (MQ=255)
                                                                                       tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                           >  2:895872/1‑121 (MQ=255)
                                                                                       tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                           <  1:895872/121‑1 (MQ=255)
                                                                                       tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAcc                                                               >  2:351229/1‑149 (MQ=255)
                                                                                        cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCg                                                              <  2:277848/149‑1 (MQ=255)
                                                                                          aaTCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGt                                                            >  1:595055/1‑149 (MQ=255)
                                                                                              cACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCaa                                                                                     >  3:545827/1‑120 (MQ=255)
                                                                                                ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAt                                                      >  2:768346/1‑149 (MQ=255)
                                                                                                ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAt                                                      <  1:354874/149‑1 (MQ=255)
                                                                                                 cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                     >  2:1034716/1‑101 (MQ=255)
                                                                                                 cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAtt                                                     >  2:617635/1‑149 (MQ=255)
                                                                                                       cGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGata                                               <  2:294109/149‑1 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACgg                                                                                  >  4:665284/1‑112 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACgg                                                                                  <  3:665284/112‑1 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCaa                                                                        <  3:962077/122‑1 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCaa                                                                        >  4:962077/1‑122 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                             <  3:480569/149‑1 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                             >  2:668070/1‑149 (MQ=255)
                                                                                                         ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                             <  1:768346/149‑1 (MQ=255)
                                                                                                             gTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGaca                                         >  4:53775/1‑149 (MQ=255)
                                                                                                              ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGacac                                        >  3:751595/1‑149 (MQ=255)
                                                                                                              ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGacac                                        >  1:690483/1‑149 (MQ=255)
                                                                                                                cttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGa                                                 >  3:510813/1‑138 (MQ=255)
                                                                                                                cttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGa                                                 <  4:510813/138‑1 (MQ=255)
                                                                                                                       ggCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTAt                               <  4:1146660/149‑1 (MQ=255)
                                                                                                                       ggCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTAt                               >  4:425520/1‑149 (MQ=255)
                                                                                                                             cAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGAt                         <  1:1013223/149‑1 (MQ=255)
                                                                                                                              aaCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGtt                                                                                                                       <  2:99328/54‑1 (MQ=255)
                                                                                                                               aCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                       >  4:894227/1‑149 (MQ=255)
                                                                                                                               aCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                       >  2:1012041/1‑149 (MQ=255)
                                                                                                                                 ggCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTg                                                  <  2:780819/120‑1 (MQ=255)
                                                                                                                                 ggCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTg                                                  >  1:780819/1‑120 (MQ=255)
                                                                                                                                  gCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                            >  1:923957/1‑93 (MQ=255)
                                                                                                                                  gCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                            <  2:923957/93‑1 (MQ=255)
                                                                                                                                      aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCtt                >  3:1081766/1‑149 (MQ=255)
                                                                                                                                      aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCtt                >  3:145495/1‑149 (MQ=255)
                                                                                                                                        taaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCa              <  4:652369/149‑1 (MQ=255)
                                                                                                                                          gtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGgaga                          <  4:247458/132‑1 (MQ=255)
                                                                                                                                          gtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGgaga                          >  3:247458/4‑135 (MQ=255)
                                                                                                                                           tcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAtt                                                     <  1:616601/105‑1 (MQ=255)
                                                                                                                                           ataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg           >  3:18573/1‑149 (MQ=255)
                                                                                                                                           ataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCAAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg           >  1:345785/1‑149 (MQ=255)
                                                                                                                                             aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAg         >  1:1025140/1‑149 (MQ=255)
                                                                                                                                             aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAg         >  3:876147/1‑149 (MQ=255)
                                                                                                                                              ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGc        >  4:655625/1‑149 (MQ=255)
                                                                                                                                              ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGc        >  4:637201/1‑149 (MQ=255)
                                                                                                                                              ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGc        <  1:922693/149‑1 (MQ=255)
                                                                                                                                               cAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAg         >  3:122228/1‑147 (MQ=255)
                                                                                                                                                    ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCt  <  2:690483/149‑1 (MQ=255)
                                                                                                                                                    ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCt  <  3:425520/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGGTACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCT  >  NC_000913/1037969‑1038265

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: