Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,405,592 A→G 8.8% M53V (ATG→GTG)  accB → biotin carboxyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,405,5920AG 8.8% 260.7 / 12.7 91M53V (ATG→GTG) accBbiotin carboxyl carrier protein
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base G (5/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TTCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCAAAAGCGT  >  NC_000913/3405443‑3405739
                                                                                                                                                     |                                                                                                                                                   
ttCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatga                                                                                                                                                     <  1:416665/150‑1 (MQ=255)
  cGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatgatg                                                                                                                                                   >  1:1049958/1‑150 (MQ=255)
  cGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatgatg                                                                                                                                                   >  1:8672/1‑150 (MQ=255)
  cGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatgatg                                                                                                                                                   <  2:1660724/150‑1 (MQ=255)
    tAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgca                                                                                                                                                 <  1:1213034/150‑1 (MQ=255)
    tAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgca                                                                                                                                                 <  1:1360009/150‑1 (MQ=255)
     aaGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgcag                                                                                                                                                <  1:736626/150‑1 (MQ=255)
      aGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgcagc                                                                                                                                               <  1:1464517/150‑1 (MQ=255)
      aGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgcagc                                                                                                                                               >  1:1637938/1‑150 (MQ=255)
       gATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcg                                                                                                                                              <  1:225777/150‑2 (MQ=255)
         ttAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCagcagc                                                                                                                                            >  2:1427094/1‑150 (MQ=255)
          tAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcgtaa                                                                                                                                           >  1:654854/1‑146 (MQ=255)
          tAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagcc                                                                                                                                           <  1:766925/150‑1 (MQ=255)
            aaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcgtaagc                                                                                                                                         <  2:654854/150‑7 (MQ=255)
            aaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagccag                                                                                                                                         >  2:139725/1‑150 (MQ=255)
            aaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagccag                                                                                                                                         >  1:1592198/1‑150 (MQ=255)
            aaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagccag                                                                                                                                         <  1:1864015/150‑1 (MQ=255)
             aaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagccagc                                                                                                                                        >  2:1745030/1‑150 (MQ=255)
                         gCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCATTCTAACGc                                                                                                                            <  2:218296/150‑1 (MQ=255)
                          cTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCa                                                                                                                           <  2:230372/150‑1 (MQ=255)
                           tGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAg                                                                                                                          <  1:1953622/150‑1 (MQ=255)
                             gTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGcc                                                                                                                        <  1:1828891/150‑1 (MQ=255)
                               tGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGc                                                                                                                         <  1:1382543/147‑1 (MQ=255)
                               tGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGc                                                                                                                         >  2:1382543/1‑147 (MQ=255)
                               tGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGc                                                                                                                      >  1:38217/1‑150 (MQ=255)
                                     aTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGc                                                                                                                >  1:931537/1‑150 (MQ=255)
                                     aTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGc                                                                                                                >  2:28372/1‑150 (MQ=255)
                                     aTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGc                                                                                                                <  1:1372227/150‑1 (MQ=255)
                                       cAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGTCGCTCCGGCGa                                                                                                              <  1:1994443/150‑1 (MQ=255)
                                        aGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGAc                                                                                                             >  2:674728/1‑150 (MQ=255)
                                         ggCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCACCTCAATCTAACGCAGCCGCTCCCGCCAcc                                                                                                            >  2:1585162/1‑150 (MQ=255)
                                           cATCTCCGAAATGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGt                                                                                                          >  2:271804/1‑150 (MQ=255)
                                            aTCTCCGACCTGGAACTTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCGTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTACGGCGACCGtt                                                                                                         <  2:392405/150‑1 (MQ=255)
                                            aTCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGtt                                                                                                         <  2:1076810/150‑1 (MQ=255)
                                            aTCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGtt                                                                                                         <  2:385184/150‑1 (MQ=255)
                                              ctcCGAACTGGAAATTTCTGACGGCGAAGAGGCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGttcc                                                                                                       <  1:1011707/150‑1 (MQ=255)
                                              ctcCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGttcc                                                                                                       <  1:969348/150‑1 (MQ=255)
                                              ctcCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGttcc                                                                                                       <  1:456042/150‑1 (MQ=255)
                                              ctcCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCATACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGCCCGttcc                                                                                                       >  2:17468/1‑150 (MQ=255)
                                                   aaCTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCa                                                                                                  <  1:139725/150‑1 (MQ=255)
                                                      tGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATgg                                                                                               <  2:559488/150‑1 (MQ=255)
                                                       ggAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGa                                                                                              >  2:1350720/1‑150 (MQ=255)
                                                       ggAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGa                                                                                              <  2:1592198/150‑1 (MQ=255)
                                                       ggAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGa                                                                                              >  2:817631/1‑150 (MQ=255)
                                                        gAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGaa                                                                                             >  1:1323689/1‑150 (MQ=255)
                                                                tGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcagc                                                                                     <  2:275128/150‑1 (MQ=255)
                                                                       gAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCg                                                                              <  2:1533882/150‑1 (MQ=255)
                                                                       gAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCCCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCg                                                                              >  2:796178/1‑150 (MQ=255)
                                                                       gAAGAGGCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCg                                                                              <  1:1307629/150‑1 (MQ=255)
                                                                         agagTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGa                                                                            <  2:1637938/150‑1 (MQ=255)
                                                                          gagTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGaa                                                                           <  1:41831/150‑1 (MQ=255)
                                                                                gTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGt                                                                     >  1:1357467/1‑150 (MQ=255)
                                                                                 tACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTg                                                                    <  1:1427094/150‑1 (MQ=255)
                                                                                 tACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTg                                                                    <  1:255531/150‑1 (MQ=255)
                                                                                  aCGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTgg                                                                   <  1:1585162/150‑1 (MQ=255)
                                                                                   cGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGt                                                                  >  1:1320678/1‑150 (MQ=255)
                                                                                         aGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATc                                                            >  1:994278/1‑150 (MQ=255)
                                                                                         aGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATc                                                            >  2:1418236/1‑150 (MQ=255)
                                                                                           ccGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGt                                                          <  2:8672/150‑1 (MQ=255)
                                                                                           ccGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGt                                                          >  1:1085483/1‑150 (MQ=255)
                                                                                               gCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGt                                                      >  1:232428/1‑150 (MQ=255)
                                                                                                  gCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTcc                                                   <  1:1149340/150‑1 (MQ=255)
                                                                                                  gCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTcc                                                   <  1:185925/150‑1 (MQ=255)
                                                                                                    tCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTcccc                                                 <  2:1049958/150‑1 (MQ=255)
                                                                                                     ccTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGCCCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCg                                                >  2:1530109/1‑150 (MQ=255)
                                                                                                        gCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtg                                             >  1:1468230/1‑150 (MQ=255)
                                                                                                         ccGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtgg                                            <  2:1323689/150‑1 (MQ=255)
                                                                                                          cGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggt                                           >  1:1630640/1‑150 (MQ=255)
                                                                                                            cAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggttg                                         <  2:1961415/150‑1 (MQ=255)
                                                                                                              aGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGATCCCCGAtggttggt                                       >  2:1767609/1‑150 (MQ=255)
                                                                                                                tttCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTAc                                     <  1:817631/150‑1 (MQ=255)
                                                                                                                 ttCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACt                                    <  2:1978235/150‑1 (MQ=255)
                                                                                                                    ccTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTc                                 <  2:1357467/150‑1 (MQ=255)
                                                                                                                    ccTGTGATGCAACAACCTTACGCTGCACCAATGATCCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTc                                 >  2:1032672/1‑150 (MQ=255)
                                                                                                                     cTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggttggttagaaat                                >  1:1634567/1‑143 (MQ=255)
                                                                                                                     cTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCt                                >  1:1402739/1‑150 (MQ=255)
                                                                                                                     cTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCt                                >  1:610825/1‑150 (MQ=255)
                                                                                                                       gtgATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTAc                              <  2:1468230/150‑1 (MQ=255)
                                                                                                                        tgATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTAcc                             <  1:1745030/150‑1 (MQ=255)
                                                                                                                           tGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCa                          <  1:717745/150‑1 (MQ=255)
                                                                                                                              aacaacCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCCTACGTTCCCCGATGGTAGGTACTTTCTACCGCAccc                       >  2:1028758/1‑150 (MQ=255)
                                                                                                                              aacaaGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTa                                      >  2:879594/1‑135 (MQ=255)
                                                                                                                              aacaaGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTa                                      <  1:879594/135‑1 (MQ=255)
                                                                                                                                  aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAg                   >  1:1418345/1‑150 (MQ=255)
                                                                                                                                  aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAg                   >  1:1033125/1‑150 (MQ=255)
                                                                                                                                   gCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGc                  >  1:1839741/1‑150 (MQ=255)
                                                                                                                                   ctggctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCCTGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGCTGGTTGGTACTTTCGCCCGCACCCCAAGc                  >  2:1025763/7‑150 (MQ=255)
                                                                                                                                       ctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg              >  2:1915887/3‑150 (MQ=255)
                                                                                                                                       ctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg              >  2:1139143/3‑150 (MQ=255)
                                                                                                                                       aCGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg              <  2:611844/150‑1 (MQ=255)
                                                                                                                                        cGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGa             >  2:1568622/1‑150 (MQ=255)
                                                                                                                                        cGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGa             >  1:1187297/1‑150 (MQ=255)
                                                                                                                                         gCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGAc            >  1:317672/1‑150 (MQ=255)
                                                                                                                                            gCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCa         <  2:1839741/150‑1 (MQ=255)
                                                                                                                                               ccAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCaaaa      >  1:691763/1‑150 (MQ=255)
                                                                                                                                                   tgatgCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCAAAAGCGt  >  1:109263/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCAAAAGCGT  >  NC_000913/3405443‑3405739

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