Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,107,717 T→A 9.3% F455I (TTC→ATC) ‡ clsC → cardiolipin synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,107,7170TA 9.3% 275.9 / 14.4 96F455I (TTC→ATC) ‡clsCcardiolipin synthase C
Reads supporting (aligned to +/- strand):  ref base T (44/43);  new base A (2/7);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGGGATGCGCAAA  >  NC_000913/1107568‑1107865
                                                                                                                                                     |                                                                                                                                                    
aaaCGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGtt                                                                                                                                                      <  2:1275899/150‑1 (MQ=255)
  aggCTGGCACAGTTTATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTc                                                                                                                                                    <  2:859858/148‑1 (MQ=255)
   cGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCt                                                                                                                                                   <  2:796564/150‑1 (MQ=255)
       ggCACAGTTAATTGATAAACGCTTTTTTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGaa                                                                                                                                               >  2:247319/1‑150 (MQ=255)
            aGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                          <  1:1028934/150‑1 (MQ=255)
             gTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGt                                                                                                                                         >  1:425814/1‑150 (MQ=255)
               tAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTa                                                                                                                                       <  2:836789/150‑1 (MQ=255)
                aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                                      <  2:839956/150‑3 (MQ=255)
                aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACACGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                                      <  2:1074675/150‑3 (MQ=255)
                 aTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttctt                                                                                                                                     <  1:555425/150‑4 (MQ=255)
                    gATACACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                                  <  2:463792/150‑7 (MQ=255)
                    gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                                  <  2:468916/150‑7 (MQ=255)
                    gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                                  <  2:436751/150‑7 (MQ=255)
                    gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                                  <  1:586821/150‑7 (MQ=255)
                    gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTc                                                                                                                                  >  2:717751/1‑150 (MQ=255)
                                 ttCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGa                                                                                                                     <  1:1151870/150‑1 (MQ=255)
                                 ttCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGa                                                                                                                     <  2:974478/150‑1 (MQ=255)
                                  tCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGat                                                                                                                    >  1:1100745/1‑150 (MQ=255)
                                   cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                                   >  2:337650/1‑150 (MQ=255)
                                   cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                                   >  2:295534/1‑150 (MQ=255)
                                   cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                                   >  2:406043/1‑150 (MQ=255)
                                         cAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGccc                                                                                                             >  1:64696/1‑150 (MQ=255)
                                          aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                          >  1:15404/1‑120 (MQ=255)
                                          aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                          <  2:15404/120‑1 (MQ=255)
                                          aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGATGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCg                                                                                                            >  1:1193979/1‑150 (MQ=255)
                                             atgatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTgg                                                                                                         <  1:1526520/150‑1 (MQ=255)
                                              tgatgCGGCCTGGCAGCTCCGTCTGGACAGGGGGGGACGGATCAACTGGGTTGATCGGCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGa                                                                                                        >  2:926600/1‑150 (MQ=255)
                                               gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                                       >  2:1389521/1‑150 (MQ=255)
                                               gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                                       >  2:1347181/1‑150 (MQ=255)
                                                 tgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATg                                                                                                     <  2:1524066/150‑1 (MQ=255)
                                                   cGGCCTGGCAGCTCGGTGTGGACAGGTGGGGACGTATCAAGTGTTTTGATCGTCATGCTATGAAAGTGATTATTCTCTAAAAAGAACCCGCGACAAGTTTCTGGAAGAGGGTTATTGTCAGACTGGCGTCGATATTGCCGGTGGAATTGt                                                                                                   <  2:421598/150‑1 (MQ=255)
                                                   cGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGt                                                                                                   <  2:181662/150‑1 (MQ=255)
                                                    ggCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGtt                                                                                                  >  1:1168541/1‑150 (MQ=255)
                                                      ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                                >  2:24928/1‑150 (MQ=255)
                                                      ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                                <  1:510239/150‑1 (MQ=255)
                                                       cTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAtt                                                                                               <  1:59705/150‑1 (MQ=255)
                                                          gCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTa                                                                                            <  2:351816/150‑1 (MQ=255)
                                                           cAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaa                                                                                           >  1:674201/1‑150 (MQ=255)
                                                            aGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaaa                                                                                          <  2:425814/150‑1 (MQ=255)
                                                             gCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAt                                                                                         >  1:1161829/1‑150 (MQ=255)
                                                                ccGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAtgtg                                                                                      >  2:931305/1‑150 (MQ=255)
                                                                 cGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTgg                                                                                     <  1:406043/150‑1 (MQ=255)
                                                                           ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                           >  1:175826/1‑150 (MQ=255)
                                                                           ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                           >  1:451556/1‑150 (MQ=255)
                                                                              ggggACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTgg                                                                        >  2:1530244/1‑150 (MQ=255)
                                                                                   cGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCtt                                                                   <  1:1408705/150‑1 (MQ=255)
                                                                                   cGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCtt                                                                   >  2:312167/1‑150 (MQ=255)
                                                                                    ggATCAACTGGGTTGATCGTGATGCAAAGTATTAGATTATTCTGTAAATAGTACCGGCCACCAGTTTCTGGAAGCGGGTTATTGTCAGACTGGCGTCGATATTGCGCGTGGATTGGTTATTGTAAATGTGTCTTAACGGGCTGGTGGttt                                                                  <  2:918169/150‑1 (MQ=255)
                                                                                     gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                                 >  2:121147/1‑150 (MQ=255)
                                                                                     gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                                 >  2:510868/1‑150 (MQ=255)
                                                                                     gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                                 <  1:24928/150‑1 (MQ=255)
                                                                                        cAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                              <  2:969659/150‑1 (MQ=255)
                                                                                         aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                             >  1:254669/1‑150 (MQ=255)
                                                                                         aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                             >  1:768311/1‑150 (MQ=255)
                                                                                              ggTTGATCGTCATGCAGAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCa                                                        <  2:1091747/150‑1 (MQ=255)
                                                                                               gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAATAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGGCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                                       <  2:255058/150‑1 (MQ=255)
                                                                                               gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                                       >  1:1136980/1‑150 (MQ=255)
                                                                                                   aTCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCtt                                                   >  1:1525443/1‑150 (MQ=255)
                                                                                                      gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                                <  1:447682/150‑1 (MQ=255)
                                                                                                      gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                                <  1:442560/150‑1 (MQ=255)
                                                                                                         aTGCAAAGAAAGTGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                             <  1:1390920/150‑1 (MQ=255)
                                                                                                         aTGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGTAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTTGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                             <  1:1396557/150‑1 (MQ=255)
                                                                                                          tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                            >  2:707308/1‑150 (MQ=255)
                                                                                                          tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                            <  1:1389521/150‑1 (MQ=255)
                                                                                                           gCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCg                                           <  1:717751/150‑1 (MQ=255)
                                                                                                            cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                              >  1:889709/1‑130 (MQ=255)
                                                                                                            cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                              <  2:889709/130‑1 (MQ=255)
                                                                                                            cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                          >  1:224893/1‑150 (MQ=255)
                                                                                                            cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                          <  1:295534/150‑1 (MQ=255)
                                                                                                             aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                                         >  2:1184448/1‑150 (MQ=255)
                                                                                                             aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                                         <  2:1168541/150‑1 (MQ=255)
                                                                                                               agaaagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGttt                                       <  2:674201/150‑1 (MQ=255)
                                                                                                                  aagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGGTTAACAACCGGTTTAcc                                    >  2:1281358/1‑150 (MQ=255)
                                                                                                                  aagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTAcc                                    <  1:312167/150‑1 (MQ=255)
                                                                                                                   agatATTATTCTCAACAAAGAACCCGCCTCCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTTGCGTACCGGGCTGGTGCTTTATCGTTCTGCCGCTTAACAACCGGTTTACCa                                   <  2:813375/150‑1 (MQ=255)
                                                                                                                   agagATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCa                                   >  2:1077285/1‑150 (MQ=255)
                                                                                                                    gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                                  <  1:1132095/150‑1 (MQ=255)
                                                                                                                    gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                                  >  1:783326/1‑150 (MQ=255)
                                                                                                                     agaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCaga                                 <  1:677336/150‑1 (MQ=255)
                                                                                                                      tattattCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                                <  2:974262/149‑1 (MQ=255)
                                                                                                                      gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                                <  1:121147/150‑1 (MQ=255)
                                                                                                                      gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                                <  1:1532119/150‑1 (MQ=255)
                                                                                                                      gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                                <  2:25003/150‑1 (MQ=255)
                                                                                                                           ttCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaa                           >  1:1426347/1‑150 (MQ=255)
                                                                                                                             ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa                         >  1:791520/1‑150 (MQ=255)
                                                                                                                             ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa                         >  1:1269805/1‑150 (MQ=255)
                                                                                                                               cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc                       >  1:548383/1‑150 (MQ=255)
                                                                                                                               cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc                       >  2:1245960/1‑150 (MQ=255)
                                                                                                                                aaaaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGATGGTGCTTTATCGTTCTCACGATTAACAACCGGTTTACCAGAGAACAAAAACt                      <  1:707308/150‑1 (MQ=255)
                                                                                                                                  aaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTc                    <  1:247319/150‑1 (MQ=255)
                                                                                                                                   ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa                   >  1:1214353/5‑150 (MQ=255)
                                                                                                                                   ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa                   >  1:1374988/5‑150 (MQ=255)
                                                                                                                                      gAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTa                >  1:1445145/1‑150 (MQ=255)
                                                                                                                                        aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                                  >  1:462739/1‑130 (MQ=255)
                                                                                                                                        aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                                  <  2:462739/130‑1 (MQ=255)
                                                                                                                                         cccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACg             >  1:205512/1‑150 (MQ=255)
                                                                                                                                          ccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACgg            >  2:524583/1‑150 (MQ=255)
                                                                                                                                                    tttCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGGTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGGGATGCGCaaa  <  2:258363/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGGGATGCGCAAA  >  NC_000913/1107568‑1107865

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: ‑