Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,361 C→T 5.8% F274F (TTC→TTT) ‡ waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,3610CT 5.8% 383.8 / 6.9 122F274F (TTC→TTT) ‡waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base C (53/61);  new base T (4/3);  total (58/64)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 124 total aligned reads displayed.

CGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTG  >  NC_000913/3809212‑3809509
                                                                                                                                                     |                                                                                                                                                    
cgccgcCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTc                                                                                                                                                      <  2:111551/150‑1 (MQ=255)
 gccgcgAGTGGGGTCCACACCGCGCGGTTTTGATTGCCACCAGCACTCACGAAGGCGCTGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTcc                                                                                                                                                     <  1:1468579/150‑1 (MQ=255)
     ccAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCGTCATCCGGAACGCTTCCCGGAt                                                                                                                                                >  1:1806890/1‑150 (MQ=255)
      cAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGAt                                                                                                                                                <  2:306592/150‑1 (MQ=255)
       aGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATg                                                                                                                                               <  1:929897/150‑1 (MQ=255)
       aGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATAAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATg                                                                                                                                               >  2:1121686/1‑150 (MQ=255)
          gggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                            <  2:736955/150‑1 (MQ=255)
           ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATgacgg                                                                                                                                           <  1:783355/150‑5 (MQ=255)
           ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGAt                                                                                                                                           >  2:701102/1‑150 (MQ=255)
            gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGAtt                                                                                                                                          >  1:1291/1‑150 (MQ=255)
            gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGAtt                                                                                                                                          >  2:654176/1‑150 (MQ=255)
             caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACCTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTa                                                                                                                                         >  1:63198/1‑150 (MQ=255)
             caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATgacgggg                                                                                                                                         <  1:238604/150‑7 (MQ=255)
               ccacACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAc                                                                                                                                       >  1:1403395/1‑150 (MQ=255)
                cacaCCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAcc                                                                                                                                      <  1:75169/150‑1 (MQ=255)
                cacaCCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAcc                                                                                                                                      <  1:1126201/150‑1 (MQ=255)
                cacaCCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAAcc                                                                                                                                      <  2:1372201/150‑1 (MQ=255)
                   acCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTg                                                                                                                                   <  1:226700/150‑1 (MQ=255)
                        ccGGTATGGATTGCCACCAGCACTCACGAAGGCTAAGAGAGTGTGGTGATCGCCGCACATCAGGCATGGTTACAGCATTTCGCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAAAGCTTCCCGGATGCGATTACCCTTGTCCGc                                                                                                                              >  2:1215891/1‑150 (MQ=255)
                          ggTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCa                                                                                                                            >  1:436253/1‑150 (MQ=255)
                             aTGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATAAGGCATTGTTACAGAAATTCCAGAATTTATTGCTCATCCTTGTACCCCGTCATCCGGAACCCTTCCCGGATGAGATTAACCTTTTCCGCCAGGc                                                                                                                         >  2:346673/1‑150 (MQ=255)
                               ggATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                                       >  1:177483/1‑150 (MQ=255)
                                gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTgg                                                                                                                      >  1:1475281/1‑150 (MQ=255)
                                gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTgg                                                                                                                      >  1:666647/1‑150 (MQ=255)
                                gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTgg                                                                                                                      <  2:1510991/150‑1 (MQ=255)
                                  ttGCCACCATCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGAc                                                                                                                    <  1:552228/150‑1 (MQ=255)
                                  ttGCCACCATCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGAc                                                                                                                    <  1:1557141/150‑1 (MQ=255)
                                  ttGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGAc                                                                                                                    <  2:143073/150‑1 (MQ=255)
                                       accaGCACTCACGAATGCGAAGAGAGTTTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCGCTAATTTATTGGTCATCCTTGTACCCCTTCATCCGTAACGCTTCCCGGATGCGTTTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               <  1:346698/150‑1 (MQ=255)
                                       accaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               <  1:620055/150‑1 (MQ=255)
                                       accaGCACTCACGAAGGAGAAGAGAGTTTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               <  1:245801/150‑1 (MQ=255)
                                        ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTCCCCCGTCATCCGTAACGCTTGCCGGATGAGATTAACCTTGTCCGCCAGGCTGTACTAAGCt                                                                                                              >  2:1115720/1‑150 (MQ=255)
                                        ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                              >  1:1709901/1‑150 (MQ=255)
                                           gCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCtata                                                                                                           <  1:1296477/150‑1 (MQ=255)
                                            cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCTGGCTGGACTAAGCTatat                                                                                                          <  1:1799649/150‑1 (MQ=255)
                                            cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTatat                                                                                                          <  1:861121/150‑1 (MQ=255)
                                             aCTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATc                                                                                                         >  2:1728765/1‑150 (MQ=255)
                                              cTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATca                                                                                                        >  2:893228/1‑150 (MQ=255)
                                               tCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcac                                                                                                       <  1:1346727/150‑1 (MQ=255)
                                                cACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcaca                                                                                                      <  1:1213990/150‑1 (MQ=255)
                                                  cGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCGTTATTTATTGATCATCCTTGTACCCCGTCATCCGGGCCGCGTCCCGGATGCGATTGTCCTTGTCCGCGAGGCTGGACTCAGCGATATcgctcg                                                                                                    >  2:1112121/1‑145 (MQ=255)
                                                   gAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGc                                                                                                   <  2:609574/150‑1 (MQ=255)
                                                   gAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCATCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGc                                                                                                   <  1:1663796/150‑1 (MQ=255)
                                                      ggCGAAGAGAGTGTGGTGATCCCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGctct                                                                                                >  1:824122/1‑150 (MQ=255)
                                                       gCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCTCACGCTCtt                                                                                               >  1:276673/1‑150 (MQ=255)
                                                        cGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTc                                                                                              >  1:1489856/1‑150 (MQ=255)
                                                           agagagTGTGGTGATCGCCGCACATCATGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCGGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAgg                                                                                           >  2:1622130/1‑150 (MQ=255)
                                                            gagagTGTGGTGATCTCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGTCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                          <  2:1355424/150‑1 (MQ=255)
                                                            gagagTGTGGTGATCGCCGCACAGCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                          <  2:1319196/150‑1 (MQ=255)
                                                             agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAgggg                                                                                         <  1:730989/150‑1 (MQ=255)
                                                                gtgtGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAg                                                                                      >  1:1448933/1‑150 (MQ=255)
                                                                gtgtGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAg                                                                                      <  2:645722/150‑1 (MQ=255)
                                                                   tggtgATCGCCGCACATCAGGCATTGTTACAGCAGTTCCCGAATTTATTGCTCATCCTGGTACCCCGTCCTCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTTCGCTATATCACACGCTCTTCAGGGGAAGTcc                                                                                   >  2:361334/1‑150 (MQ=255)
                                                                   tggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGACCGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTcc                                                                                   <  1:141244/150‑1 (MQ=255)
                                                                    ggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTccc                                                                                  >  2:1781726/1‑150 (MQ=255)
                                                                    ggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTccc                                                                                  >  2:535155/1‑150 (MQ=255)
                                                                     gtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTcccc                                                                                 <  1:1025978/150‑1 (MQ=255)
                                                                      tgATCGCCGCACATCAGTCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCTGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTTGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCt                                                                                <  1:211334/150‑1 (MQ=255)
                                                                       gATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                               <  1:1042702/150‑1 (MQ=255)
                                                                       gATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCGTTATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                               >  2:950650/1‑150 (MQ=255)
                                                                         tCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccc                                                                             >  2:920979/1‑149 (MQ=255)
                                                                         tCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTcca                                                                             <  2:1258885/150‑1 (MQ=255)
                                                                           gccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccacc                                                                           >  1:185399/1‑150 (MQ=255)
                                                                           gccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccacc                                                                           <  1:257566/150‑1 (MQ=255)
                                                                               cacaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCa                                                                       <  1:535155/150‑1 (MQ=255)
                                                                                     cggCTTTGTTACAGCAATTCCCGAATTTATTGCTGATCCTGGGACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAgg                                                                 <  2:1697398/149‑1 (MQ=255)
                                                                                     aGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAgg                                                                 <  1:1540897/150‑1 (MQ=255)
                                                                                     aGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGt                                                                <  2:1806890/150‑1 (MQ=255)
                                                                                      ggCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAggg                                                                >  2:1463350/1‑149 (MQ=255)
                                                                                      ggCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGt                                                                <  2:267033/150‑1 (MQ=255)
                                                                                         aTTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTtgt                                                             <  1:1027033/150‑1 (MQ=255)
                                                                                          ttGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTtgtg                                                            <  1:135068/150‑1 (MQ=255)
                                                                                             ttACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGtt                                                         <  2:430227/150‑1 (MQ=255)
                                                                                              tACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTg                                                        >  1:982654/1‑150 (MQ=255)
                                                                                                cAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGc                                                      <  1:1079210/150‑1 (MQ=255)
                                                                                                cAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCAGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTACCCTACACCAGCACGCAGGTTGTGGTTGGc                                                      <  1:199320/150‑1 (MQ=255)
                                                                                                  gCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGa                                                    >  2:828376/1‑150 (MQ=255)
                                                                                                   cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGAt                                                   <  2:1381419/150‑1 (MQ=255)
                                                                                                   cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGAt                                                   <  1:471816/150‑1 (MQ=255)
                                                                                                         ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTTTCCGCCATCATGGACTAAGCTATATCACACGCTCTTAAGGGGAAGTCCACTCCACCAGCACGAAGGATGTGTGTGGCGATACGagg                                             >  1:367169/1‑148 (MQ=255)
                                                                                                         ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATg                                             <  1:654176/150‑1 (MQ=255)
                                                                                                                tATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGt                                      >  2:629118/1‑150 (MQ=255)
                                                                                                                tATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGt                                      <  1:850008/150‑1 (MQ=255)
                                                                                                                 aTTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGtt                                     <  2:446591/150‑1 (MQ=255)
                                                                                                                 aTTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGtt                                     <  1:375226/150‑1 (MQ=255)
                                                                                                                    gCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATTCGATGGGCGAgatgat                                  >  2:1809022/1‑150 (MQ=255)
                                                                                                                     cTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCGGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATg                                 <  2:1680310/150‑1 (MQ=255)
                                                                                                                           cTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCTATACGATGGGCGAGTTGATGtttatg                           >  1:1256092/1‑146 (MQ=255)
                                                                                                                             ggTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCAGACGCACTACAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGGTGGCGATACGATGGGCGAGTTGATGTTACTGTa                         >  2:1297322/1‑150 (MQ=255)
                                                                                                                               tACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATg                       >  2:1433858/1‑150 (MQ=255)
                                                                                                                                   atcgcATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCAt                   >  2:638694/5‑150 (MQ=255)
                                                                                                                                   atcgcATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTCGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCAt                   >  1:1684772/5‑150 (MQ=255)
                                                                                                                                      gcATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGc                >  1:1589174/2‑150 (MQ=255)
                                                                                                                                          ccGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAt            <  1:893228/150‑1 (MQ=255)
                                                                                                                                           cGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAtc           >  2:968793/1‑150 (MQ=255)
                                                                                                                                              aaCGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCCGCCCGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTTCTGTATGGCATTGCCGATCTCg        >  2:1786865/1‑150 (MQ=255)
                                                                                                                                                cGCTTCCCGGATGCGATTAACCTTGTCCGCCAGTCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGcc      <  2:1403395/150‑1 (MQ=255)
                                                                                                                                                cGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGcc      <  2:550211/150‑1 (MQ=255)
                                                                                                                                                  cTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCtt    >  2:1295292/1‑150 (MQ=255)
                                                                                                                                                    tCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTttg  >  1:1613393/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTG  >  NC_000913/3809212‑3809509

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