Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,324 C→A 5.7% H81N (CAT→AAT)  torC → trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,3240CA 5.7% 245.0 / 3.8 96H81N (CAT→AAT) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (48/42);  new base A (3/3);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAAT  >  NC_000913/1058176‑1058472
                                                                                                                                                    |                                                                                                                                                    
gATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGATTGCCATGACTGTc                                                                                                                                                      <  2:725297/149‑1 (MQ=255)
     tATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATc                                                                                                                                                 >  1:33136/1‑149 (MQ=255)
      aTTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATcc                                                                                                                                                <  2:188830/149‑1 (MQ=255)
      aTTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATcc                                                                                                                                                >  1:1057947/1‑149 (MQ=255)
      aTTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATcc                                                                                                                                                >  1:283514/1‑149 (MQ=255)
       ttGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGAATATccg                                                                                                                                               >  1:226071/1‑148 (MQ=255)
       ttGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATccc                                                                                                                                               >  1:1163178/1‑149 (MQ=255)
       ttGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTAGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATccc                                                                                                                                               >  1:965358/1‑149 (MQ=255)
         gCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGACCGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCccgc                                                                                                                                             >  1:144832/1‑149 (MQ=255)
         gCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCccgc                                                                                                                                             >  2:380245/1‑149 (MQ=255)
         gCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCccgc                                                                                                                                             >  2:10541/1‑149 (MQ=255)
          ccACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACTGTCGGTGCATTTCCAGAACGCCTACGGCGTGCGATCTGATTTCCATGACTGTAATATCccgcc                                                                                                                                            >  2:666066/1‑149 (MQ=255)
          ccACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCccgcc                                                                                                                                            >  1:953164/1‑149 (MQ=255)
             cacGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATACACAGTCCGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCAGAATGCCATGACTGTCATATCCCGCCGGa                                                                                                                                         >  1:1092527/1‑149 (MQ=255)
              acGTTGGGATCAAAGTCCCCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGGGAGCTGAATGCCATGACTGTCATATCCCGCCGGat                                                                                                                                        >  1:679666/1‑149 (MQ=255)
               cGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGata                                                                                                                                       <  2:631071/149‑1 (MQ=255)
                     gATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGAATATCCGGCGGGatatgacag                                                                                                                                 <  2:226071/149‑6 (MQ=255)
                     gATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGACACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAg                                                                                                                                 >  1:1081545/1‑149 (MQ=255)
                       tCAAAGTCACCAGCACAACCGAATTTTGTGTCAGGTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGc                                                                                                                               <  1:547835/149‑1 (MQ=255)
                         aaaGTCAACAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAAACGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGCGCTGCATGCCATGCCTGTCATATCCCGCCGGAAATTCCACGCAt                                                                                                                             >  2:660843/1‑149 (MQ=255)
                          acgTCACCAGAAAAAGCGAATTTTGTGTCAGTTGCCACAGTATGCATCCGGTGTATGAAGATTATAAACAGTCGGTGCATTTCCAGAACGCCTCCCGCGTGCGAGCGGAAAGCCATGACTGTCATATCCCGCCCGATATTCCAGGCATg                                                                                                                            >  2:188479/3‑149 (MQ=255)
                           aGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAAAATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGAATATCCGGCGGGATATGACAGTCATgg                                                                                                                           <  1:342428/149‑1 (MQ=38)
                             tCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGACGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGAATATCCGGCGGGATATGACAGTCATggca                                                                                                                         <  2:1153463/149‑3 (MQ=38)
                               aCCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGaa                                                                                                                       >  1:511663/1‑149 (MQ=255)
                                ccAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAg                                                                                                                      >  2:896237/1‑149 (MQ=255)
                                   gCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAgcgc                                                                                                                   <  1:135120/149‑1 (MQ=255)
                                   gCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGGGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAgcgc                                                                                                                   <  1:910002/149‑1 (MQ=255)
                                       aaCCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAAc                                                                                                               <  1:1001308/149‑1 (MQ=255)
                                       aaCCGAATTTTGTGTCAGTAGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACAGTCATAACCCGCCGGATATTCCAGGCATGGTGAAGAGCAAAc                                                                                                               >  2:27441/1‑149 (MQ=255)
                                        aCCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACt                                                                                                              >  2:1184247/1‑149 (MQ=255)
                                        aCCGAATTATGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACt                                                                                                              <  1:262709/149‑1 (MQ=255)
                                         ccGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTg                                                                                                             >  2:428669/1‑149 (MQ=255)
                                         ccGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTg                                                                                                             >  2:584773/1‑149 (MQ=255)
                                         ccGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTg                                                                                                             >  1:671012/1‑149 (MQ=255)
                                            aagtttGTGTCAGTTGTCACAGTATGCAACCGGTGGAGGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGaa                                                                                                          <  2:780876/146‑1 (MQ=255)
                                             aTTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAg                                                                                                         <  1:511945/149‑1 (MQ=255)
                                              ttttGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGc                                                                                                        >  1:312886/1‑149 (MQ=255)
                                              ttttGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGc                                                                                                        >  2:73241/1‑149 (MQ=255)
                                                 tgtgGCAGTTGCCACAGTAGGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAg                                                                                                     <  1:852846/149‑1 (MQ=255)
                                                   tgtCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCa                                                                                                   >  2:620002/1‑149 (MQ=255)
                                                    gtCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTACAGAACGCCTCCGGCGTGCGAGCTGAATGCAATGACTGTAATATCCCGCCTGATATTCCAGGCATTGTGAATCTCAAACTGGAAGCGAGCaa                                                                                                  >  1:710686/1‑149 (MQ=255)
                                                     tCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAAt                                                                                                 <  2:953164/149‑1 (MQ=255)
                                                     tCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAAt                                                                                                 >  1:1062736/1‑149 (MQ=255)
                                                     gcAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAAt                                                                                                 <  1:388380/148‑1 (MQ=255)
                                                       aGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACCCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGa                                                                                               <  1:660843/149‑1 (MQ=255)
                                                        gTTGCCACAGTATGCAACCGGTGTATGAAGATTATAAACAGTCGGTTCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGat                                                                                              <  1:985464/149‑1 (MQ=255)
                                                        gTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATTACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGat                                                                                              >  1:665998/1‑149 (MQ=255)
                                                        gTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGat                                                                                              <  1:445336/149‑1 (MQ=255)
                                                            ccACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGAGAGCTGAATGCCAAGACTGTCATATCCCGCCGGATATTCAAGGCATGGTGAAGCGCAAACTGGAAGCGAGAAATGAAATCt                                                                                          >  1:827635/1‑149 (MQ=255)
                                                            ccACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTACCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCAAATCCCGCCGGATATTCCAGGCATCGTGAAGCGCAAACTGGAAGCAAGCAATGATATCt                                                                                          >  1:270914/1‑149 (MQ=255)
                                                              acaGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTAc                                                                                        >  2:575497/1‑149 (MQ=255)
                                                              acaGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTAc                                                                                        <  1:719726/149‑1 (MQ=255)
                                                                aGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCa                                                                                      <  1:146524/149‑1 (MQ=255)
                                                                     gCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGATGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCt                                                                                 <  1:1170893/149‑1 (MQ=255)
                                                                         ccGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTAt                                                                             <  1:624993/149‑1 (MQ=255)
                                                                                aTGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCAc                                                                      <  2:255495/149‑1 (MQ=255)
                                                                                 tGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACt                                                                     >  1:1004564/1‑149 (MQ=255)
                                                                                 tGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACATTTATTGCTCACt                                                                     >  1:609584/1‑149 (MQ=255)
                                                                                  gaagaaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTc                                                                    <  2:33136/149‑1 (MQ=255)
                                                                                  gaagaaTATAAACAGTCGGTGCATTTCCAGAACGCCTACGGCGTGCGAGCTGAATGCCATGACTGTCATCTCCCGCCGGATATTACATGCATGGTGCAGCGCAAAATGGAAGCGAGCAATGATATCTTCAAGACCTTTATTgctcgctc                                                                    >  2:894120/1‑149 (MQ=255)
                                                                                   aagaaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTcc                                                                   >  1:107490/1‑149 (MQ=255)
                                                                                     gaaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGGGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCAt                                                                 <  1:1047162/149‑1 (MQ=255)
                                                                                      aaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCAtt                                                                >  1:893464/1‑149 (MQ=255)
                                                                                      aaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCAtt                                                                <  2:506097/149‑1 (MQ=255)
                                                                                                    gTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGaaaa                                                  >  1:430078/1‑149 (MQ=255)
                                                                                                      gCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAAt                                                >  2:835605/1‑149 (MQ=255)
                                                                                                      gCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTCATACAACTGAAAAAt                                                >  2:178939/1‑149 (MQ=255)
                                                                                                        aTTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTc                                              >  2:31756/1‑149 (MQ=255)
                                                                                                         tttCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGCTATCTACAAGCCCGTTATTGCTACCTCCATTGATTCACCTGAAAATTTCg                                             >  2:289235/1‑149 (MQ=255)
                                                                                                         tttCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCg                                             >  1:1084561/1‑149 (MQ=255)
                                                                                                          ttCCAGAAAGCCTCCTGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCATGCATGGTGAATCGCAAACTGGAAGAGAGCAATTATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGa                                            <  1:981726/149‑1 (MQ=255)
                                                                                                              aGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTTAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGcc                                        >  2:727986/1‑149 (MQ=255)
                                                                                                              aGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGcc                                        <  2:309183/149‑1 (MQ=255)
                                                                                                              aGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGcc                                        <  2:107490/149‑1 (MQ=255)
                                                                                                              aGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGcc                                        <  1:366994/149‑1 (MQ=255)
                                                                                                                     cTCCGGCGTGCTAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAcgcg                                 <  2:647151/149‑1 (MQ=255)
                                                                                                                          gCGTGCGAGCGGAAGTCCCTGACTGTCATATACCTCCGGATATTCCAGTCATGGTGAAGCGAAAACTGGAAGCGATCAATTATATCTACCAGACCTTTATTGCTCACTCAATTGATACAACTGAAAAATTCGAAGCCAAACGCGCGGac                            <  2:740149/149‑2 (MQ=255)
                                                                                                                           ggTGCGATCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGCAGCGCAAACTGGAAGCGTGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATCCACCTGAAAAATTCGAAGCCAAACGCGCGGAAc                           <  2:93405/148‑1 (MQ=255)
                                                                                                                             tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACtt                         <  2:283514/149‑1 (MQ=255)
                                                                                                                               cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGc                       <  1:762639/149‑1 (MQ=255)
                                                                                                                                gAGCTGAATGCCATGACTGTCATATCCCTCAGGATATTCCAGGCATGGTGAATCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGcc                      <  2:923673/149‑1 (MQ=255)
                                                                                                                                 aGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGAGAGCAATGATATCAACCAGACCTTTATTGCTCACTCCATTGATACATCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCg                     >  1:576240/1‑149 (MQ=255)
                                                                                                                                      aaTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGt                >  2:815369/1‑149 (MQ=255)
                                                                                                                                      aaTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGt                >  1:1157490/1‑149 (MQ=255)
                                                                                                                                      aaTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGt                <  2:1057947/149‑1 (MQ=255)
                                                                                                                                        tGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGa              >  2:753859/1‑149 (MQ=255)
                                                                                                                                        tGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGa              <  2:365466/149‑1 (MQ=255)
                                                                                                                                        tGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGa              >  1:125602/1‑149 (MQ=255)
                                                                                                                                         gCCATGACTGTCATATCCCGCCGGATAATCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGaa             <  1:916363/149‑1 (MQ=255)
                                                                                                                                          ccATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAAt            <  1:380245/149‑1 (MQ=255)
                                                                                                                                           cATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATg           >  1:89105/1‑149 (MQ=255)
                                                                                                                                           cATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATg           >  1:859844/1‑149 (MQ=255)
                                                                                                                                            aTGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAATCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATgc          <  2:572648/149‑2 (MQ=255)
                                                                                                                                              gACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGgc        <  2:529323/149‑1 (MQ=255)
                                                                                                                                              gACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGgc        <  1:933689/149‑1 (MQ=255)
                                                                                                                                                cTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAACCGCGCGGAACTTGCCGAGCGTGAATGGgcgc      <  2:671012/149‑1 (MQ=255)
                                                                                                                                                cTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGgcgc      <  2:1163178/149‑1 (MQ=255)
                                                                                                                                                    cATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGCGCAATCATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAACCATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAAt  <  2:710686/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAAT  >  NC_000913/1058176‑1058472

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